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ABSTRACT 

Nucleate  pool  boiling  heat  transfer  is  an  integral  part 
of  any  vapor-compression  refrigeration  cycle.  With  a  view  to 
improving  overall  cycle  efficiency,  the  heat  transfer 
performance  in  the  evaporator  can  be  improved  by  using 
enhanced  boiling  surfaces.  This  thesis  looks  at  the  pool 
boiling  characteristics  of  R-114  (presently  used  in  large 
shipboard  AC  systems)  from  10  enhanced  single  copper  tubes  and 
compares  performance  with  a  smooth  copper  tube.  Since  small 
amounts  of  oil  escape  into  the  refrigerant  as  it  passes 
through  the  compressor  of  a  refrigeration  system,  tests  have 
also  been  conducted  with  up  to  10%  (by  weight)  of  a  miscible 
oil  to  see  what  effect  this  may  have  on  overall  evaporator 
performance . 
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NOMENCLATURE 


A  area 

Ab  tube  outside  surface  area  of  active  boiling 

section 

Ac  cross-sectional  area  of  tube 

Cp  specific  heat 

D  diameter 

Di  tube  inside  diameter 

Do  tube  outside  diameter 

Dl  diameter  at  the  position  of  the  thermocouple 

D2  outer  diameter  of  the  boiling  tube 

g  gravitational  acceleration 

h  heat-transfer  coefficient 

hfg  latent  heat  of  vaporization 

I  current 

Is  current  reading  by  AC  Current  Sensor 

k  thermal  conductivity  of  liquid 

kc  thermal  conductivity  of  copper 

L  active  boiling  tube  length 

Lu  non-boiling  length  of  the  test  tube 

Nu  Nusselt  number 

p  tube  outside  wall  perimeter 

Pr  Prandtl  Number 

Q  heat-transfer  rate  from  boiling  surface 

Qf  heat-transfer  rate  through  one  nonboiling  end 

Qh  heat-transfer  rate  from  cartridge  heater 

q"  heat  flux 

Ra  Rayleigh  number 

T  temperature 

Tavg  average  wall  temperature  at  the  thermocouple 

location 

Tc  critical  temperature 

Tf  film  temperature 

xiv 


Tn  temperature  of  the  thermocouple  location 

Tsat  saturation  temperature 

Two  outer  wall  temperature  of  the  boiling  test  tube 

V  voltage  across  the  cartridge  heater 

Vs  voltage  reading  by  AC-DC  true  RMS  converter 

a  thermal  diffusivity 

P  volumetric  thermal  expansion  coefficient 

5  uncertainty  in  measurement  and  calibration 

6  superheat 

[i  dynamic  viscosity 

V  kinematic  viscosity 
p  density 
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I .  INTRODUCTION 

A.  BACKGROUND 

The  United  States  Navy  currently  uses  R-114  in 
approximately  950  centrifugal  chilled-water  air-conditioning 
plants.  Ongoing  efforts  to  reduce  the  cost  and  increase  the 
efficiency  of  these  plants  have  been  directed  towards  the  use 
of  enhanced  heat-transfer  tube  surfaces  and  alternative 
refrigerants.  The  choice  of  alternative  refrigerant  must  be 
"ozone  friendly."  A  possible  "drop  in"  replacement  for  R-114 
has  been  identified:  R-124.  Studies  of  the  boiling  performance 
of  R-124  would  be  complemented  well  by  a  baseline  of  data  on 
R-114. 

B.  ADVANTAGES  OF  USING  R-114 

The  use  of  R-114  in  Naval  applications  has  been  driven  by 
concerns  for  both  cost  and  safety.  R-114  operates  under 
moderate  pressures  and  has  one  of  the  lowest  toxicity  ratings 
among  all  refrigerants.  System  components  are  significantly 
smaller  and  more  lightweight  which  is  desirable  for  shipboard 
use  and  maintenance.  While  R-ll  is  widely  used  in  industry, 
its  operating  requirements  make  it  unsuitable  onboard  ship.  R- 
11  operates  under  vacuum  and  has  a  tendency  to  absorb  ambient 
moisture.  This  acidic  combination  of  R-ll  and  moisture  attacks 
tubes  and  causes  corrosion.  R-114  is  considerably  stable  to 


heat  and  may  be  used  with  all  standard  ferrous  and  non-ferrous 
metals  except  zinc,  magnesium  and  aluminum  alloys  containing 
appreciable  amounts  of  zinc  and  magnesium.  The  thermophysical 
properties  of  R-114  are  provided  in  Appendix  A. 

C.   THESIS  OBJECTIVES 

Safety  requirements  for  future  studies  of  R-124 
necessitated  the  disassembly  and  relocation  of  the  single-tube 
apparatus  described  in  this  thesis.  Past  studies  using  the 
apparatus  have  varied  considerably  in  operating  procedure, 
thus  it  became  desirable  to  establish  a  consistent  operating 
procedure  in  the  acquisition  of  all  data  on  R-114.  With  the 
foregoing  discussion  in  mind,  the  objectives  of  this  thesis 
were  to: 

1.  Relocate  and  calibrate  the  single  tube  apparatus. 

2 .  Operate  the  apparatus  to  produce  repeatable  data  and 
confirm  earlier  results. 

3.  Obtain  boiling  data  for  R-114  and  R-114/oil  mixtures  on 
ten  tube  surfaces  to  be  used  as  a  database  for  future 
studies  on  enhanced  surfaces  and  alternative  "ozone 
safe"  refrigerants. 


II.  REVIEW  OF  REFRIGERANT  AND  REFRIGERANT/OIL  POOL 

BOILING  BEHAVIOR 


A.   NUCLEATE  POOL  BOILING  PERFORMANCE  OF  PURE  R-114 

The  heat-transfer  performance  of  pure  R-114  on  a  smooth 
tube  in  the  natural  convection  region  of  its  boiling  curve  can 
be  reasonably  predicted  using  well  known  correlations, 
Churchill  and  Chu  [Ref.  1]  and  more  recently  Churchill  and 
Usagi  [Ref.  2]  .  A  lack  of  complete  understanding  of  the 
nucleate  boiling  mechanism  makes  it  more  difficult  the  to 
predict  the  performance  in  the  nucleate  boiling  region.  Well 
known  correlations,  such  as  that  by  Rohsenow  [Ref.  1],  and 
more  recently,  Stephan-Adelsalam  [Ref.  3]  are  approximate  and 
should  be  compared  with  experimental  data  whenever  possible. 

Many  pool-boiling  experiments  have  been  carried  out  on 
finned  tubes  (typically  19  to  40  fpi)  in  pure  refrigerants. 
For  various  refrigerants,  tube  material,  fin  density  and 
system  pressure,  typical  enhancements  in  the  heat  transfer 
coefficient  (over  smooth  tubes  values)  of  up  to  3  have  been 
reported,  the  refrigerants  with  better  wetting  characteristics 
giving  the  greater  enhancement.  Interestingly  enough,  the 
boiling  performance  of  finned  tubes  becomes  very  similar  to 
that  of  smooth  tubes,  if  the  heat  flux  is  recalculated  on  the 
basis  of  total  outside  wetted  surface  area  (rather  than  root 
diameter,  which  is  commonly  used) . 


Re-entrant  cavity  surfaces  differ  from  finned  surfaces  in 
that  they  provide  a  higher  density  of  stable  nucleation  sites 
with  minimum  wall  superheat  (Yilmaz  and  Westwater  [Ref.  4] 
Carnavos  [Ref.  5]  ,  Marto  and  Lepere  [Ref.  6]  and  Wanniarachchi 
et  al.  [Refs.  7,  8]).  Significant  improvements  in  heat 
transfer  over  both  finned  and  smooth  tubes  have  been  shown  for 
re-entrant  cavity  surfaces. 

Boiling  enhancements  were  as  large  as  14.6  and  6.4  for 
High  Flux  and  Turbo-B  tubes,  respectively,  compared  to  smooth 
tube  data. 

Webb  [Ref.  9]  stated  that  the  key  to  the  high  performance 
of  the  re-entrant  surfaces  can  be  attributed  to  three  factors: 
(1)  re-entrant  cavity  within  a  critical  size  range,  (2) 
interconnected  cavities,  and  (3)  nucleation  sites  of  a  re- 
entrant shape.  If  the  cavities  are  interconnected,  adjacent 
cavities  can  activate  each  other.  Re-entrant  cavities  provide 
a  stable  vapor  trap,  which  can  remain  active  at  low  values  of 
superheat . 

During  nucleate  boiling  of  R-114  at  one  atmosphere, 
Wanniarachchi  et  al.  [Refs.  7  and  8]  reported  data  for  four 
enhanced  tubes:  a  porous  coated  High  Flux,  Thermoexcel-E, 
Thermoexcel-HE  and  a  26  fpi  GEWA-T.  Enhancements  of  9.1,  8.2, 
6,8  and  4.4  respectively  were  obtained  at  a  heat  flux  of  30 
kW/m2.  The  enhancement  provided  by  the  High  Flux  tube  reduced 
to  6  at  higher  fluxes  of  about  100  kW/m2,  but  it  still 
outperformed  the  other  three  tubes.  At  heat  fluxes  below  18 


kW/m2,  both  the  Thermoexcel  tubes  performed  the  best.  However, 
they  noted  that  the  uncertainty  in  the  data  increased  as  both 
superheat  and  heat  flux  decreased. 

Whereas  the  majority  of  pool-boiling  experiments  use 
electrically  heat  tubes,  McManus  et  al .  [Ref.  10]  used  hot 
water  to  heat  both  a  High  Flux  and  Turbo-B  tube  in  a  pool  of 
R-114.  The  range  of  heat  flux  was  very  limited.  However,  at  a 
heat  flux  of  40  Kw/m2,  enhancements  of  14.6  and  6.4 
respectively,  were  obtained  over  smooth  tube  values. 

B.   NUCLEATE  POOL  BOILING  PERFORMANCE  OF 
REFRIGERANT-OIL  MIXTURES 

Most  refrigeration  systems  use  hermetically  sealed 
compressors  which  allow  small  quantities  of  oil  to  escape  into 
the  working  fluid.  This  oil  tends  to  collect  in  the  flooded 
evaporator,  giving  up  to  10%  oil  in  some  cases.  The  nucleate 
pool  boiling  mechanism  is  further  complicated  when  a  less 
volatile  component,  such  as  miscible  refrigerant  oil,  is  mixed 
with  the  refrigerant.  Design  and  modelling  rely  heavily  on 
accurate  experimental  data  covering  a  wide  range  of  operating 
conditions . 

Compared  with  pure  refrigerants,  there  is  little  data  for 
boiling  of  refrigerant /oil  mixtures  from  finned  surfaces. 
Chaddock  [Ref.  11]  gives  a  thorough  summary  of  such  mixtures 
through  1975.  Sauer  et  al.  [Ref.  12]  report  data  from  19  fpi 
copper  tube  using  R-ll  with  two  different  oils.  For  oil 


concentrations  less  than  3%,  their  results  indicated  no 
difference  in  boiling  heat  transfer  performance  from  that  of 
the  pure  refrigerant.  For  oil  concentrations  greater  than  5%, 
significant  degradation  in  performance  was  shown.  They 
concluded  that  finned  tubes  performed  best  at  low  wall 
superheats  and  that  this  performance  was  not  impaired  by  up  to 
3%  oil  addition. 

Murphy  [Ref.  13]  carried  out  experiments  on  a  26  and  40 
fpi  tube  using  R-114/oil  mixtures.  Similar  to  Sauer  [Ref.  11], 
no  difference  in  performance  was  observed  for  oil 
concentrations  up  to  3%.  Indeed  in  some  cases,  at  a  higher 
heat  flux  of  about  100  kW/m2,  a  3%  oil  concentration  showed  a 
small  increase  in  performance  over  the  pure  refrigerant  case. 
The  mechanism  for  this  improvement  is  not  well  understood,  but 
has  been  found  by  a  number  of  investigators  for  both  flow  and 
pool  boiling.  At  a  10%  oil  concentration,  both  tubes  exhibited 
a  30%  degradation  in  performance  versus  the  pure  refrigerant 
case,  which  agreed  with  Sauer' s  findings. 

Wanniarachi  et  al.  [Refs.  7,  8]  reported  the  performance 
of  High  Flux,  Thermoexcel-E,  Thermoexcel-HE  and  GEWA-T  tubes 
in  R-114/oil  mixtures  at  one  atmosphere.  Oil  concentrations  up 
to  10%  (by  weight)  were  studied.  A  marked  increase  in  the  wall 
superheat  values  required  for  the  incipience  of  nucleate 
boiling  was  observed  in  the  presence  of  oil.  Enhancements 
mentioned  previously  for  pure  R-114  were  reduced  at  a  3%  oil 
concentration,  to  7.5,  5.5,  5.5  and  3.8  respectively.  With  10% 


oil  concentrations,  further  reductions  to  5.5,  3.8,  3,8  and 
2.9  were  observed.  At  the  highest  heat  fluxes  (>  45  kW/m2)  and 
for  oil  concentrations  greater  than  6%,  the  performance  of  the 
High  Flux  tube  degraded  at  a  rate  that  left  the  GEWA-T  tube  as 
the  best  performer. 


III.  DESCRIPTION  OF  EXPERIMENTAL  APPARATUS 

A.   OVERVIEW  OF  THE  SYSTEM 

Wanniarachchi  et  al.  [Ref.  7],  Murphy  [Ref.  13]  and 
Karasabun  [Ref.  14]  provide  complete  details  of  the  original 
apparatus  used  in  previous  studies.  Disassembly  of  the 
original  apparatus  to  facilitate  relocation  and  subsequent 
modifications  to  meet  present  site  requirements  justify  a  full 
description  of  the  apparatus.  The  original  R-12  cooling 
subsystem  was  supplanted  by  a  more  efficient  R-502  system.  The 
auxiliary  heaters  and  plastic  shield  were  removed,  however 
power  extensions  were  kept  available  for  their  use  in  future 
studies.  The  Variac  power  supply  to  the  tube  heater  was 
calibrated  and  emplaced  into  a  panel  enclosure.  See  Appendix 
D,  Apparatus  Calibration  for  more  details.  A  more  modern 
Hewlet-Packard  Computer/Data  Acquisition  system  was  installed 
and  the  data  reduction  program  DRP8  was  updated  to  account  for 
variac  calibration  and  program  language  changes. 

An  overall  schematic  representation  of  the  experimental 
apparatus  is  shown  in  Figure  3.1,  and  a  photograph  is  shown  in 
Figure  3.2.  The  apparatus  essentially  consists  of  8 
(components:  1)  a  Pyrex-glass  tee  for  the  pool  boiling  of  the 
R-114  liquid,  "the  evaporator";  (2)  a  Pyrex-glass  tee  for  the 
condensing  of  the  R-114  vapor,  "the  condenser";  (3)  a  R-114 


liquid  reservoir;  (4)  a  refrigerant /oil  subsystem  composed  of 
an  oil  reservoir  and  a  graduated  oil  cylinder;  (5)  a  cooling 
subsystem  composed  of  a  1/2  ton  R-502  refrigeration  plant,  2 
positive  displacement  pumps  and  a  30  gallon  water/ethylene 
glycol  sump;  (6)  a  vacuum  pump;  (7)  a  data  acquisition  and 
instrumentation  system;  (8)  an  aluminum/Plexi-glass  framework 
within  which  components  (1)  through  (4)  were  housed. 

The  apparatus  was  operated  as  follows:  liquid  R-114  was 
boiled  in  the  evaporator  (1)  ,  passed  as  vapor  through  the 
aluminum  L-shaped  connection  and  condensed  in  the  condenser 
(2) .  R-114  condensate  was  returned  by  gravity  via  the 
distribution  tube  to  the  condenser.  The  R-502  refrigeration 
plant  maintained  the  water/ethylene  glycol  mixture  in  the  sump 
between  -12  and  -18  °C.  Cooling  of  the  sump  was  accomplished 
by  a  countercurrent  flow  heat  exchanger  through  which  an  8-gpm 
turbine-type  pump  recirculated  the  water/ethylene  glycol 
mixture.  An  identical  positive  displacement  pump  circulated 
the  cooled  mixture  through  the  copper  condensing  coils  in  the 
condenser  via  the  control  valve  VC. 

Prior  to  dismantling  the  evaporator  in  order  to  change  out 
boiling  tubes,  the  liquid  R-114  was  boiled  out  of  the 
evaporator,  condensed  on  the  coils  of  the  condenser  then 
recovered  to  the  R-114  liquid  reservoir  rather  than  returned 
to  the  evaporator.  Upon  reassembly,  the  evaporator  was 
recharged  with  liquid  R-114  by  gravity  from  the  reservoir.  The 


level  in  the  reservoir  could  be  topped  off  via  the  condenser 
from  an  external  supply  cylinder  through  valve  V10. 

Exact  quantities  of  oil  could  be  added  to  the  liquid  R-114 
in  the  lower  glass  tee  from  the  graduated  cylinder  (5)  by 
gravity.  The  graduated  cylinder  was  filled  with  oil  by  gravity 
from  the  oil  reservoir  through  valves  V14  and  V2  .  Refilling 
the  oil  reservoir  was  done  through  valve  V3  located  at  the  top 
of  the  framework.  All  oil  left  over  from  the  recovery  of  the 
liquid  R-114  was  drained  out  of  the  lower  glass  tee  and 
disposed  of. 

The  boiling  tube  was  mounted  horizontally  in  the 
evaporator  and  held  in  place  by  Teflon  bushings  with  O-ring 
inserts.  The  whole  system  was  connected  to  the  vacuum  pump  (6) 
through  valve  V8  in  order  to  remove  noncondensible  gases.  A 
relief  valve  mounted  on  the  aluminum  L  connection  was  ser  zc 
a  gage  pressure  limit  of  20  psi,  which  is  50%  less  thai)  the 
manufacturer-recommended  working  pressure  limit  of  the  ?yrex 
glass  tees.  All  connections  were  assembled  with  Swagelok 
fittings  sealed  with  Teflon  ferrules  to  ensure  leak  tightness. 
Copper  tubing  of  3/8  inch  diameter  was  used  for  all  piping, 
except  in  the  oil  section,  where  1/4  inch  copper  tubing  was 
used. 
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B.   BOILING  TEST  SECTION 

1 .  Evaporator 

The  evaporator  is  a  T-shaped  container  made  of  Corning 
Pyrex  glass.  Figure  3.3  shows  a  schematic  of  the  evaporator 
with  dimensions.  Pyrex  glass  has  several  advantages:  it  is 
corrosion-resistent ,  transparent,  has  a  smooth  interior 
surface  (this  minimized  nucleate  boiling  at  the  inner  surface 
of  the  container)  and  is  stronger  with  temperature  and 
pressure  variations  compared  to  ordinary  glass.  The 
manufacturer  recommended  operating  pressure  limit  of  the  Pyrex 
glass  tee  is  30  psi  gage.  Each  end  of  the  evaporator  was 
fitted  with  a  cast-iron  flange  and  a  gasket.  A  detailed  sketch 
of  the  cast-iron  flange  is  shown  in  Figure  3.3  also.  To  reduce 
weight  upon  the  glass  tees,  two  aluminum  flanges,  210  mm  in 
diameter  and  12.7  mm  in  thickness,  were  bolted  to  the  cast- 
iron  flanges.  All  fittings  were  connected  through  the  aluminum 
flanges  to  the  Pyrex-glass  evaporator. 

2.  Test  Tube 

A  schematic  drawing  of  the  test  tube  is  shown  in 
Figure  3.4.  The  boiling  tube  was  held  in  place  by  two  Teflon 
bushings,  which  were  attached  to  the  aluminum  flanges  at  both 
ends  of  the  T-shaped  evaporator.  Four  studs  were  used  in  order 
to  attach  the  Teflon  bushing  on  the  aluminum  flange.  The 
Teflon  bushing  and  both  ends  of  the  test  tube  were  sealed  by 
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means  of  two  viton  O-rings  that  were  placed  between  the 
aluminum  flange  and  the  Teflon  plug  at  each  end. 
3 .  Boiling  Tubes 

A  smooth,  hard-copper  tube,  15.9  mm  (5/8  in.)  in  outer 
diameter,  12.7  mm  (1/2  in)  in  inside  diameter  and  431.3  mm  (17 
in.)  in  length,  was  used  to  provide  comparison  data  for  the 
nucleate  pool  boiling  heat-transfer  performance  of  the 
enhanced  tubes  tested.  All  of  the  tubes  were  heated  and 
instrumented  identically.  Dimensions  of  each  tube  are  provided 
in  Table  V  of  Appendix  E,  Sample  Calculations.  Wanniarachchi 
et  al.  [Refs.  8,  9]  and  Murphy  [Ref.  14]  provide  complete 
details  of  the  tube  surfaces. 

The  ten  tubes  tested  were: 

1 .  Smooth 

2.  GEWA-K  2  6  fpi  (Figure  3.5) 

3.  GEWA-K  40  fpi  (Figure  3.5) 

4.  GEWA-T  19  fpi  (Figure  3.5) 

5.  GEWA-T  2  6  fpi  (Figure  3.5) 

6.  GEWA-YX  26  fpi  (Figure  3.6) 

7.  High  Flux  (Figure  3.7) 

8.  Thermoexcel-HE  (Figure  3.8) 

9.  Thermoexcel-E  (Figure  3.9) 

10.  Turbo-B  (Figure  3.9) 
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The  heater  was  a  1000-Watt,  240-Volt  stainless-steel 
cartridge,  6.35  mm  (1/4  in.)  in  outer  diameter  and  203.2  mm  (8 
in)  in  length.  It  was  inserted  into  a  copper  sleeve,  which  was 
1/4  inch  in  inside  diameter,  1/2  inch  in  outer  diameter  and  8 
inches  in  length.  In  order  to  provide  a  uniform  heat  flux,  the 
cartridge  heater  and  the  copper  sleeve  were  soldered  together. 
To  minimize  the  thermal  contact  resistence,  the  copper  sleeve 
was  then  tinned  (the  cartridge  heater  was  used  to  maintain  the 
solder  in  molten  state)  and  inserted  as  a  unit  into  the  middle 
portion  of  the  test  tube.  The  active  boiling  length  of  the 
test  tube  was  therefore  eight  inches  in  the  middle  portion  of 
the  boiling  tube.  In  order  to  compute  the  actual  average  heat 
flux  in  the  heated  portion  of  the  tube,  a  suitable  correction 
was  applied  for  the  heat  lost  by  natural  convection  out  both 
ends.  See  Appendix  E,  Sample  Calculation  for  details. 

To  measure  the  wall  surface  temperature  of  the  boiling 
tube,  8  thermocouples  (Type-T  Teflon  coated  copper-constantan; 
0.01  inch  in  diameter)  were  inserted  into  eight  grooves 
machined  on  the  outside  of  the  copper  sleeve.  As  shown  in 
Figure  3.10,  these  thermocouples  were  located  at  different 
axial  and  circumferential  locations.  The  exact  locations  of 
the  thermocouples  and  dimensions  of  the  thermocouple  grooves 
are  given  in  Figure  3.10.  All  the  thermocouple  grooves  were 
axially  machined  from  the  location  of  the  thermocouple  hot 
junctions  to  the  nearest  end  of  the  copper  sleeve. 
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C.  CONDENSER  SECTION 

The  condenser  was  also  a  T-shaped  container  made  of 
Corning  Pyrex  glass.  It  was  identical  to  the  evaporator.  The 
position  of  the  condenser  can  be  seen  in  Figure  3.1.  R-114  was 
condensed  on  a  helical  copper  coil,  which  was  inserted  in  the 
Pyrex-glass  condenser.  Copper  tubing  (3/8  inch  outer  diameter) 
was  fabricated  into  a  three  inch  diameter  coil.  The  active 
condensation  length  was  estimated  to  be  15  ft. 

The  top  portion  of  the  condenser  was  connected  to  a 
portable,  mechanical  vacuum  pump  to  remove  noncondensible 
gases  from  the  apparatus.  The  bottom  of  the  condenser  was 
connected  to  the  evaporator  via  valve  V5  in  order  to  return 
the  condensed  R-114  liquid  to  either  the  evaporator  or  the 
reservoir.  The  coolant,  i.e.,  water/ethylene  glycol  mixture, 
entered  from  the  top  portion  of  the  condenser.  The  condenser 
was  placed  vertically  and  connected  to  the  vapor  outlet  of  the 
evaporator  using  a  fabricated,  L-shaped  aluminum  tube,  two 
inches  in  diameter.  A  Bourdon  gage  with  a  range  up  to  a  gage 
pressure  of  150  psi  and  a  relief  valve  which  was  set  to  20  psi 
gage  were  mounted  on  the  L-tube. 

D.  OIL  ADDING  SECTION 

To  study  the  boiling  performance  of  R-114/oil  mixtures,  a 
cylindrical  aluminum  oil  reservoir,  6  inches  in  diameter  and 
6  inches  in  height,  and  a  graduated  glass  oil  cylinder  were 
installed  above  the  evaporator  as  shown  in  Figure  3.1.  The 
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graduated  cylinder  was  355  mm  in  length  and  had  a  diameter  of 
25.4  mm.  A  resolution  of  0.5  ml  was  achieved  with  the 
graduated  cylinder.  The  oil  reservoir  was  connected  to  the  oil 
cylinder  through  valves  V14  and  V2 .  Refilling  the  oil 
reservoir  was  accomplished  through  valve  V3  located  on  top  of 
the  reservoir.  The  addition  of  oil  into  the  evaporator  was 
achieved  through  VI  by  gravity.  The  oil  cylinder  was  refilled 
by  gravity  from  the  oil  reservoir  located  above. 

E.   COOLING  SECTION 

1 .  Water-Ethylene  Glycol  Mixture  Tank 

In  order  to  store  the  water-ethylene  glycol  mixture, 
a  special  tank  was  manufactured.  The  total  volume  of  the  tank 
was  9243  in3  (18.84  in  x  18.84  in  x  26.04  in)  and  it  was  made 
of  0.5  inch  thick  Plexiglas  sheet.  All  sides  of  the  tank  were 
glued  together  with  methylene  chloride  solution.  The  joints 
were  held  together  with  small  screws  for  extra  strength.  The 
low  thermal  conductivity  of  Plexiglas  was  especially  suited  to 
minimize  heat  transfer  (from  room  to  water-ethylene  glycol 
mixture)  through  the  tank  walls.  The  tank  was  placed  on  top  of 
square  wood  platform  on  the  floor  and  all  sides  were  insulated 
with  7/8  in  thick  rubber  sheets.  The  cooling  mixture  contained 
13  gallons  of  ethylene-glycol  and  25  gallons  of  distilled 
water.  The  freezing  point  of  this  mixture  was  about  -25  deg  C. 
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2 .  R-502  Refrigeration  Plant 

A  1/2  ton  R-502  refrigeration  plant  was  installed  to 
cool  the  water/ethylene  glycol  mixture.  Figure  3.12  shows  a 
schematic  of  the  R-502  refrigeration  plant.  It  consists  of  a 
compact-type,  air-cooled  condenser,  a  compressor,  a  receiver, 
a  filter-drier  unit,  a  pressure  regulator,  a  temperature 
control  switch  and  a  thermostatic  expansion  valve.  The 
evaporator  of  the  R-502  refrigeration  plant  was  constructed 
using  a  3.8  inch  copper  tube  which  was  run  through  a 
countercurrent  heat  exchanger  opposite  the  water-ethylene 
glycol  mixture.  The  temperature  of  the  water/ethylene  glycol 
mixture  was  controlled  by  both  a  thermostatic  expansion  valve 
and  a  temperature  control  switch.  The  R-502  refrigeration 
plant  was  adjusted  to  keep  the  temperature  of  the  cooling 
liquid  at  about  -17  deg  C. 

3 .  Pump  and  Control  Valve 

An  8  gpm,  115  VAS  Burks  turbine-type,  positive- 
displacement  pump  was  installed  on  the  floor  and  the  1  inch 
diameter  suction  side  of  the  pump  was  directly  connected  to 
the  water/ethylene  glycol  tank.  Cooling  liquid  was  pumped  from 
the  tank  to  the  condenser  through  the  control  valve  VC .  Also, 
a  by-pass  valve  V9  was  placed  before  the  control  valve  VC  on 
the  discharge  line.  The  use  of  the  by-pass  line  served  two 
important  purposes:  (1)  it  avoided  overloading  the  pump  in  the 
event  valve  VC  is  completely  closed,  and  (2)  it  provided 
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proper  mixing  for  the   "warm"   stream  returning  from  the 

condenser.  The  by-pass  valve  required  adjustment  only  at  the 

highest  heat   fluxes  to   ^satisfy'   proper  cooling  of  the 
condenser . 

F.   R-114  RESERVOIR 

An  aluminum  cylindrical  reservoir,  9  inches  in  diameter 
and  10  inches  in  height,  was  placed  vertically  between  the 
evaporator  and  condenser  in  order  to  store  R-114  as  a  liquid. 
The  liquid  level  of  the  R-114  could  be  observed  by  means  of  a 
sight  glass  attached  on  the  reservoir.  The  R-114  reservoir  was 
connected  to  the  vapor  line  through  valve  V7  and  to  the  liquid 
line  through  valve  V6.  See  Figure  3.1  for  the  arrangement  of 
the  reservoir. 

6 .   CHAMBER 

An  aluminum  frame  (42.13  in  x  20.07  in  x  24.02  in)  was 
constructed  to  locate  all  the  parts  of  the  apparatus,  except 
the  cooling  section.  All  four  vertical  sides  of  the  frame  were 
covered  with  1/2  inch  thick  Plexiglas  sheets  and  both  left  and 
right  sides  were  provided  with  hinges  to  enable  easy  access  to 
the  components  of  the  apparatus.  Aluminum  and  plywood  plates 
were  used  to  cover  the  bottom  and  top  sides  of  the  frame, 
respectively.  The  valve  bodies  of  VI  through  V8  were  placed 
inside  of  the  front  Plexiglass  box.  The  whole  frame  was  placed 
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above  the  water-ethylene  glycol  tank  with  aluminum  supports  so 
that  the  system  was  very  compact . 

One  cf  the  main  advantages  of  this  chamber  was  that  the 
~er.pera~ure   surrounding  the  evaporator   could  re   reduced 

relative  to  the  ambient  temperature.  Also,  in  case  of 
emergency,  the  uhick  Plexiglas  chamber  provides  a  safety 
barrier  to  personnel  and  data  acquisition  equipment. 

H.   INSTRUMENTATION 

1 .  Power  Measurement 

A  240  volt  AC  source  was  used  as  the  power  supply,  and 
it  was  adjusted  by  a  variac  in  the  range  of  0-220  volts  and 
0-5  amperes  according  to  the  desired  heat  flux  at  the  surface 
of  the  boiling  tube.  Power  input  to  the  boiling  tube  was 
measured  with  an  AC  current  inductive  sensor  (output  in  volts) 
and  a  voltage  sensor  directly  coupled  to  the  variac  output. 
The  voltage  sensor  was  run  through  a  AC-DC  true  R.M.S. 
converter  which  output  a  proportional  signal  in  volts.  Figure 
3.12  shows  a  schematic  representation  of  the  power-measurement 
setup.  Both  the  AC  current  sensor  and  the  AC-DC  true  R.M.S 
converter  were  connected  to  the  data  acquisition/control  unit. 
2  .  Temperature  Measurement 

Various  temperatures  were  monitored  throughout  the 
system  to  include: 


1.   Boiling   tube   wall   (8   thermocouples   at   various 
longitudinal  and  circumferential  positions) 
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2.  Liquid  temperature  (one  thermocouple) 

3.  Vapor  temperature  (two  thermocouples),  and 

4.  Water/ethylene   glycol   mixture   temperature    (one 
thermocouple) 

The  locations  of  the  wall  thermocouples  in  the  sleeve  of  the 
boiling  tube  are  shown  in  Figure  3.5.  The  liquid  and  vapor 
thermocouples  were  inserted  into  the  three  specially 
manufactured  thermocouple  housings.  Figure  3.14  shows  a 
schematic  of  these  thermocouple  housings.  The  stainless-steel 
portion  minimizes  (owing  to  low  thermal  conductivity)  errors 
resulting  from  the  axial  conduction  of  heat  from  the 
surroundings.  The  copper  tip  of  the  housing  helps  to  minimize 
the  temperature  drop  from  the  area  being  measured  to  the 
thermocouple  location  (owing  to  the  high  thermal  conductivity 
of  copper) . 

All  the  temperature  measurements  were  accomplished  by 
30  gage  copper-constantan  thermocouples.  Each  thermocouple 
measurement  was  read  directly  by  a  Hewlett-Packard  34  97A  data 
acquisition  system,  which  was  controlled  by  a  Hewlett-Packard 
9826  computer.  Each  thermocouple  was  scanned  twenty  times  over 
5  seconds  and  the  readings  averaged  to  obtain  a  more  accurate 
measurement . 
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Figure  3.1   Schematic  of  Single  Tube  Apparatus 
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Figure  3.2   Photograph  of  Experimental  Apparatus 
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a)  Corning  Pyrex  Glass  Evaporator  (D  x  Dl  =  402x51  mm, 
L  =  178  mm,  Li=  127  mm) 
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b)  Cast  Iron  Flange  and  Gasket  (d^s  190mm,  d0 =  210mm, 
L  =  14mm,  A  =  21°  ) 


Figure  3.3   Schematic  of  the  Pyrex  Glass  Evaporat 
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GEWA-T  FIN 


(a)  Standard  fin  (GEWA-K)  and  GEWA-T  fin 
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(b)  GEWA-T  fin 
Figure  3.5   Schematic  and  Photograph  of  Fins 
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(a)    GEWA-YX   tube 


(b)  GEWA-YX  fin 
Figure  3.6   Photographs  of  GEWA-YX  Tube  and  Fin 
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VAPOR  BUBBLE   TRAPPED 
IN  ACTIVE  SITE 


Figure  3.7   Schematic  and  Photograph  of  High  Flux  Surface 
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Figure  3.8  Schematic  and  Photograph  of 
Thermoexcel-HE  Surface 


27 


(a)  Thermoexcel-E 


(b)  Turbo-B 


Figure  3.9   Photographs  of  Re-Entrant  Cavity  Surfaces 
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IV.  DATA  ACQUISITION  AND  REDUCTION 

A.  DATA  ACQUISITION  AND  STORAGE 

A  Hewlett-Packard  3852A  automatic  data  acquisition  unit 
was  used  to  read  temperatures  from  the  thermocouples  and  to 
read  current  and  voltage  values  of  the  cartridge  heater  from 
the  AC  current  sensor  and  the  AC-DC  true  RMS  converter.  A 
Hewlett-Packard  9000  Series  computer  was  used  to  control  the 
measurements  done  by  the  3852A,  analyze  the  collected  data  and 
store  it. 

The  iterative  data  collection/reduction  program  DRP8  was 
loaded  and  run.  Information  was  entered  through  the  keyboard 
to  prompt  the  data  acquisition  unit  when  to  take  data.  Channel 
assignments  are  listed  in  Table  I.  The  raw  data  was  processed 
and  transferred  to  a  file  on  the  computer's  hard  drive.  A 
printout  of  the  data  reduction  was  also  provided. 

B.  DATA  REDUCTION 

Following  data  acquisition  for  each  data  point,  results 
were  computed  according  to  the  stepwise  procedure  outlined  in 
the  next  section,  and  then  printed  out  using  a  Hewlett-Packard 
Inkjet  printer.  Graphs  were  plotted  using  commercial  software. 
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STEPWISE  DATA-COLLECTION  AND  SOLUTION  PROCEDURE 


1.  Input  name  of  user-specified  file  to  be  stored  on 
computer  hard  drive. 

2.  Select  tube  type  (all  dimensions  of  the  boiling  test 
tube  are  included) . 

3.  Input  desired  heat  flux  (W/m2)  and  saturation 
temperature  (°C)  of  the  boiling  liquid. 

4.  Set  desired  heat  flux  by  adjusting  VARIAC  rheostat. 

5.  Set  saturation  temperature  by  adjusting  flow  of  coolant 
through  condenser  coils  with  control  valve  VC. 

6.  Once  saturation  temperature  is  achieved,  wait  for 
steady-state  conditions  (at  least  5  minutes)  prior  to 
taking  data  readings. 

7.  Prompt  data  acquisition  unit  to  scan  all  channels 
listed  in  Table  I  (thermocouples,  AC  current  sensor  and 
AC-DC  true  RMS  voltage  converter) .  All  channel  readings 
are  made  in  volts  and  stored  in  user  specified  file. 

8.  Convert  these  voltage  readings  to  corresponding  units 

(temperature  in  °C,  current  in  amperes,  voltage  in 
volts) .  Karasabun  [Ref .  15]  describes  in  detail 
conversion  equation  for  temperature.  See  Appendix  B- 
Data  Acquisition  Apparatus  Calibration  for  details  of 
current  and  voltage  conversion. 

9.  Compute  the  heat-transfer  rate  from  the  cartridge 
heater.  See  Appendix  D-Sample  Calculations. 

10.  Compute  the  average  wall  temperature  of  the  boiling 
tube  and  calculate  the  wall  superheat  of  the  current 
data  set.  See  Appendix  D-Sample  Calculations. 

11.  Compute  the  physical  properties  of  R-114  using  given 
correlations  at  film  temperature.  See  Appendix  D-Sample 
Calculations . 

12.  Compute  the  natural-convection  heat-transfer 
coefficient  of  R-114  from  the  smooth  non-boiling  ends 
of  the  test  tube.  See  Appendix  D-Sample  Calculations. 

13.  Compute  heat  losses  from  the  non-boiling  ends  of  the 
boiling  tube.  See  Appendix  D-Sample  Calculations. 
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14.  Calculate  the  heat  flux  from  tube  to  the  refrigerant. 
See  Appendix  D-Sample  Calculations. 

15.  Calculate  the  heat  transfer  coefficient  of  the  R-114 
from  the  boiling  tube.  See  Appendix  D-Sample 
Calculations . 

16.  Store  the  heat  flux  versus  wall  superheat  values  for 
each  data  set  in  user-specified  file. 

17.  Use  commercial  software,  Harvard  Graphics,  to  plot 
data . 


TABLE  I.   HP  3497A  CHANNEL  ASSIGNMENTS 


Channel 

Assignment 

500-507 

Boiling  tube  wall  temperature 

508-509 

Boiling  liquid  temperature 

510 

Vapor  temperature 

511 

Sump  temperature 

512 

RMS  voltage  Tube  heaters 

513 

Tube  heater  current  sensor 

514 

Auxiliary  heater  current 
sensor 
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V.  EXPERIMENTAL  PROCEDURE 

A.   PREPARATION 

1 .  Vacuum  test  of  the  Apparatus 

Upon  changing  the  tube  and  reassembling  the  apparatus, 
the  system  was  evacuated  down  to  approximately  25  in.  Hg.  by 
a  portable  mechanical  vacuum  pump  through  valves  VI 1  and  V8 . 
The  apparatus  was  allowed  to  stand  for  30  minutes  under  vacuum 
to  observe  any  possible  leaks.  During  this  period, 
thermocouples  were  connected  to  the  data  acquisition  system, 
the  electrical  integrity  of  the  tube  heater  was  checked  and 
the  tube  heater  connected  to  the  variac  panel.  If  a  leak  was 
detected  (i.e.  an  observed  rise  in  pressure),  a  pressure  test 
was  conducted  to  isolate  the  source.  (See  V.A.2.) 

2 .  Pressure  Test  of  the  Apparatus 

In  order  to  detect  sources  of  leaks,  the  apparatus  was 
pressurized  with  air  up  to  15  psi  gage  and  a  soap/water  bubble 
test  was  carried  out  upon  suspect  joints.  All  detected  leaks 
were  successfully  and  systematically  fixed. 

3.  Charging  the  Apparatus  with  R-114 

Upon  completing  the  pressure  tests  (if  any) ,  the 
evaporator  of  the  apparatus  was  filled  with  R-114  liquid  up  to 
a  marked  position  (20  mm  above  the  top  of  the  boiling  tube) 
using  the  following  procedure. 
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1.  Valve  V7  was  opened  to  allow  the  pressures  within  the 
reservoir  and  the  evacuated  apparatus  to  equalize.  It 
was  noted  that  the  liquid  R-114  within  the  reservoir 
boiled  since  the  equilibrium  pressure  was  less  than  the 
R-114  vapor  pressure  at  room  temperature. 

2.  Valve  V6  was  opened  to  allow  the  liquid  R-114  to  drain 
into  the  evaporator  by  gravity. 

3.  At  the  designated  evaporator  level,  valves  V6  and  V7 
were  closed  to  isolate  the  reservoir. 


4.  Degassing  and  Data  Acquisition  Channel  Check 


1.  With  the  condenser  operating,  the  boiling  tube  was  run 
at  a  high  flux  of  about  100  kW/m2  for  ten  minutes  in 
order  to  de-gas  both  the  tube  surface  and  refrigerant. 
Care  was  taken  not  to  exceed  a  pressure  of  15  psi  gage 

(well  below  manufacturer  recommendation)  within  the 
apparatus . 

2.  Any  noncondensible  gases  present  collected  in  the  top 
of  the  condenser.  These  were  then  removed  using  the 
vacuum  pump  through  valves,  V8  and  VI 1. 

3.  The  program,  SETUP8,  was  run  to  check  output  on  all 
thermocouple  channels  as  well  as  power  output.  If  an 
erroneous  reading  was  detected,  then  appropriate  action 
was  taken.  In  the  case  of  a  faulty  boiling  tube  wall 
thermocouple,  the  faulty  thermocouple  was  tagged  for 
future  exclusion  in  data  acquisition  procedure.  (Each 
thermocouple  had  X  thermocouples  embedded  in  its  tube 
wall  and  the  loss  of  one  made  no  difference  to  the 
average  wall  temperature) . 

4.  The  apparatus  was  shut  down  and  the  boiling  tube 
allowed  to  soak  in  the  refrigerant  overnight.  This 
allowed  for  the  surface  to  become  wetted. 


B.   NORMAL  OPERATION 

The  following  procedure  was  followed  to  obtain  the  heat- 
transfer  coefficient  of  pure  R-114  and  R-114/oil  mixtures  from 
the  boiling  tube  being  tested: 
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1.  The  R-502  refrigeration  unit  was  operated  1-2  hours  in 
advance  in  order  to  reduce  the  temperature  of  the 
water-ethylene  glycol  sump  to  a  minimum  of  -8  °C. 

2.  The  data  acquisition/control  unit,  computer  and  variac 
panel  were  switched  on. 

3.  The  computer  program  Setup8  was  loaded  and  run: 

a.  All  the  data  acquisition  channels  were  rechecked. 

b.  A  power  output  of  approximately  5W  was  input  to 
the  boiling  tube.  This  has  been  measured  to 
produce  a  heat  flux  of  500  to  600  W/m2  in  all  tubes 
tested  which  facilitated  reaching  the  first  data 
point  on  the  increasing  flux  run. 

c.  The  temperature  of  the  refrigerant  was  slowly 
reduced  down  to  2.2  °C  by  circulating  a  small 
amount  of  coolant  through  the  condenser,  regulated 
by  control  valve,  VC.  This  reduction  was  monitored 
to  coincide  with  the  reduction  in  sump 
temperature . 

4.  The  iterative  computer  program  DRP8  was  loaded  and 
run . 

5.  Two  data  runs  were  made:  an  increasing  heat  flux  run 
and  a  decreasing  heat  flux  run  wherein  the  heat  flux 
was  incrementally  increased  or  decreased.  Ten  different 
heat  fluxes  (0.6,  1,  2,  3.5,  6,  10,  20,  35,  60,  and  100 
kW/m2)  were  selected  for  each  increasing  or  decreasing 
heat  flux  data  run. 

6.  The  desired  heat  flux  and  saturation  temperature  were 
input  to  the  program  as  reference  points  for  each  data 

point . 

7.  The  variac  was  adjusted  to  set  the  actual  heat  flux 
which  was  continuously  compared  to  the  desired  heat 
flux  by  the  program  until  they  agreed  within  2  percent. 

8.  The  control  valve  VC  was  adjusted  to  regulate  the  flow 
of  cooling  liquid  through  the  condenser  to  maintain 
nearly  constant  saturation  temperature  at  a  given  heat 
flux.  Desired  versus  actual  saturation  temperatures 
were  monitored  continuously  by  the  program  until  they 
agreed  ±  0.1  °C. 

9.  For  each  data  point,  conditions  in  the  evaporator  were 
allowed  to  stabilize  for  at  least  5  minutes  prior  to 
measuring  the  raw  data.  The  following  raw  data  were 
measured  and  stored  in  a  user  specified  file:  boiling 
tube   wall    thermocouple   readings,    liquid   bulk 
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thermocouple  readings,  vapor  thermocouple  readings, 
sump  thermocouple  readings,  current  sensor  readings, 
and  voltage  sensor  readings. 

10.  Two  data  points  were  taken  at  a  given  heat  flux  and 
saturation  temperature  to  comprise  a  data  set.  The 
following  processed  data  were  recorded  as  a  printout: 
wall  temperatures  of  the  boiling  section,  liquid  bulk 
temperatures,  vapor  temperature,  sump  temperature,  wall 
superheat,  heat  transfer  coefficient  of  the  R-114/R- 
114-oil  mixture  and  the  actual  heat  flux. 

11.  For  each  data  set,  the  above  procedure  beginning  with 
step  5  was  repeated. 

12.  Table  II  provides  a  listing  of  all  data  runs. 


TABLE  II.  LISTING  OF  DATA  RUNS 


Data 
Run  # 

Data  File 

Tube 

Heat 
Flux 

Oil 
(%) 

Purpose 

1 

DATA630I5 

HF 

I 

0 

CAL 

2 

DATA630D5 

HF 

D 

0 

CAL 

3 

DATA702I5 

HF 

I 

0 

CAL 

4 

DATA7  02D5 

HF 

D 

0 

CAL 

5 

DATA0705I5 

HF 

I 

0 

CALX DATA 

6 

DATA0705D5 

HF 

D 

0 

CALX DATA 

7 

DATA0709I4 

SM 

I 

0 

DATA 

8 

DATA07  0  9D4 

SM 

D 

0 

DATA 

9 

DAT0710I6 

GK-40 

I 

0 

DATA 

10 

DAT0710D6 

GK-40 

D 

0 

DATA 

11 

DAT0711I7 

GK-26 

I 

0 

DATA 

12 

DAT0711D7 

GK-26 

D 

0 

DATA 

13 

DAT0715I8 

GT-19 

I 

0 

DATA 

14 

DAT0715D8 

GT-19 

D 

0 

DATA 

15 

DAT0716I9 

GT-26 

I 

0 

DATA 

39 


Data 
Run  # 

Data  File 

Tube 

Heat 
Flux 

Oil 
(%) 

Purpose 

16 

DAT0716D9 

GT-26 

D 

0 

DATA 

17 

DAT0717I11 

T-HE 

I 

0 

DATA\DG 

18 

DAT0717D11 

T-HE 

D 

0 

DATA 

19 

DAT0722I12 

TB 

I 

0 

DATA 

20 

DAT0722D12 

TB 

D 

0 

DATA 

21 

DAT0723I09 

GY-26 

I 

0 

DATA 

22 

DAT0723D09 

GY-26 

D 

0 

DATA 

23 

DAT0724I11 

T-HE 

I 

0 

DATA\DG 

24 

DAT0724D11 

T-HE 

D 

0 

DATA\REP 

25 

DAT0803I4 

SM 

I 

0 

DATAXREP 

26 

DAT0803D4 

SM 

D 

0 

DATAXREP 

27 

DAT0804I4 

SM 

I 

3 

DATA 

28 

DAT0804D4 

SM 

D 

3 

DATA 

29 

DAT0805I4 

SM 

I 

10 

DATA 

30 

DAT0805D4 

SM 

D 

10 

DATA 

31 

DAT0806I7 

GK-2  6 

I 

0 

DATAXREP 

32 

DAT0806D7 

GK-26 

D 

0 

DATAXREP 

33 

DAT0807I7 

GK-26 

I 

3 

DATA 

34 

DAT0807D7 

GK-26 

D 

3 

DATA 

35 

DAT08  0  8I7 

GK-26 

I 

10 

DATA 

36 

DAT0808D7 

GK-26 

D 

10 

DATA 

37 

DAT0809I6 

GK-40 

I 

0 

DATAXREP 

38 

DAT0809D6 

GK-40 

D 

0 

DATAXREP 

39 

DAT0810I6 

GK-40 

I 

3 

DATA 

40 

DAT0810D6 

GK-4  0 

D 

3 

DATA 

41 

DAT0812I6 

GK-40 

I 

10 

DATA 

42 

DAT0812D6 

GK-4  0 

D 

10 

DATA 

43 

DAT0815I8 

GT-19 

I 

0 

DATAXREP 

40 


Data 
Run  # 

Data  File 

Tube 

Heat 
Flux 

Oil 
(%) 

Purpose 

44 

DAT0815D8 

GT-19 

D 

0 

DATA\REP 

45 

DAT0816I8 

GT-29 

I 

3 

DATA 

46 

DAT0  816D8 

GT-19 

D 

3 

DATA 

47 

DAT0820I8 

GT-19 

I 

10 

DATA 

48 

DAT0820D8 

GT-19 

D 

10 

DATA 

49 

DAT0821I9 

GT-26 

I 

3 

DATA 

50 

DAT0821D9 

GT-26 

D 

3 

DATA 

51 

DAT0822I9 

GT-26 

I 

10 

DATA 

52 

DAT0822D9 

GT-26 

D 

10 

DATA 

53 

DAT0823I10 

T-E 

I 

3 

DATA 

54 

DAT0823D10 

T-E 

D 

3 

DATA 

55 

DAT0824I10 

T-E 

I 

10 

DATA 

56 

DAT0824D10 

T-E 

D 

10 

DATA 

57 

DAT0825I9 

GY-26 

I 

3 

DATA 

58 

DAT0825D9 

GY-26 

D 

3 

DATA 

59 

DAT0826I9 

GY-26 

I 

10 

DATA 

60 

DAT0826D9 

GY-26 

D 

10 

DATA 

61 

DAT0827I5 

HF 

I 

3 

DATA 

62 

DAT0827D5 

HF 

D 

3 

DATA 

63 

DAT0828I5 

HF 

I 

10 

DATA 

64 

DAT0828D5 

HF 

D 

10 

DATA 

65 

DAT0829I12 

TB 

I 

3 

DATA 

66 

DAT0829D12 

TB 

D 

3 

DATA 

67 

DAT0830I12 

TB 

I 

10 

DATA 

68 

DAT0830D12 

TB 

D 

10 

DATA 

69 

DAT0831I11 

T-HE 

I 

3 

DATA 

70 

DAT0831D11 

T-HE 

D 

3 

DATA 

71 

DAT0901I11 

T-HE 

I 

10 

DATA 
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Data 
Run  # 

Data  File 

Tube 

Heat 
Flux 

Oil 
(%) 

Purpose 

72 

DAT0901D11 

T-HE 

D 

10 

DATA 

73 

DAT0906I10 

T-E 

I 

0 

DATA 

74 

DAT0906D10 

T-E 

D 

0 

DATA 

SM  =  Smooth   GK-2  6  =  GEWA-K  2  6  fpi  GK-40  =  GEWA-K  40  fpi 

HF  =  High  Flux      GT-19  =  GEWA-T  19  fpi  GT-2  6  =  GEWA-T  2  6  fpi 

TB  =  Turbo-B  GY-2  6  =  GEWA-YX  2  6  fpi     T-E  =  Thermoexcel-E 

T-HE  =  Thermoexcel-HE 

I  =  Increasing  Flux  D  =  Decreasing  Flux 


CAL  =  CALIBRATION  DATA 
REP  =  REPEATABILITY  DATA 


DATA  =  BASELINE  DATA 
DG  =  DEGASSING  DATA 
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VI.  RESULTS  AND  DISCUSSION 

A.  REPRODUCIBILITY 

To  test  the  reproducibility  of  the  experimental  data,  9 
data  runs  were  repeated  on  different  days.  All  data  runs  were 
conducted  at  a  saturation  temperature  of  2.2  °C.  Figure  6.1 
shows  a  comparison  of  two  pairs  of  runs,  13,  14,  and  43,  44 
for  increasing  and  decreasing  heat  flux.  It  can  be  seen  that 
there  is  excellent  agreement  between  both  the  increasing  and 
decreasing  heat  flux  runs. 

B.  COMPARISON  OF  DATA  WITH  PREDICTION 

To  check  the  general  validity  of  the  data  obtained,  four 
well-known  correlations,  two  for  natural  convection  and  two 
for  nucleate  pool  boiling,  were  plotted  and  compared  with  the 
smooth  tube  data  for  pure  R-114.  The  two  natural  convection 
correlations  used  were  Churchill  and  Chu  [Ref.  1]  and  a  more 
recent  update  Churchill  and  Usagi  [Ref.  2] .  The  two  nucleate 
pool  boiling  correlations  used  were:  Rohsenow  [Ref.  1]  and 
Stephan-Abdelsalam  [Ref.  3]  .  One  of  the  first  correlations 
developed,  the  Rohsenow  correlation  is  based  upon  the  physics 
of  pool  boiling  heat  transfer.  It  applies  only  for  clean 
surfaces  and  when  used  to  estimate  heat  flux,  errors  can 
amount  to  ±  100%.  The  Stephan-Abdelsalam  correlation  for 
refrigerants  was  developed  from  application  of  regression 
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analysis  to  a  large  amount  of  existing  data.  Application  to  a 
large  number  of  pure  substances  of  widely  differing  thermal 
properties  is  remarkably  simple.  Figure  6.2  shows  the  plotted 
data  and  correlations.  It  can  be  seen  from  the  plot  that  there 
is  excellent  agreement  for  natural  convection  and  not  bad 
agreement  in  the  nucleate  boiling  region.  The  slightly  worse 
agreement  for  nucleate  boiling  can  be  attributed  to  the  more 
complex  mechanisms  at  work  in  this  region. 

C.   EFFECTS  OF  DEGASSING 

Early  experiments  (Runs  1-6)  conducted  with  the  High  Flux 
tube  pointed  out  the  significant  effect  of  degassing  both  the 
tube  surface  and  refrigerant.  It  was  observed  that  tubes  that 
were  not  vigorously  boiled  after  placement  into  the  apparatus 
and  not  allowed  to  settle  overnight,  tended  to  nucleate  at 
very  low  heat  fluxes  the  following  day  when  conducting  an 
increasing  heat  flux  run.  Similar  effects  were  observed  if  the 
apparatus  was  allowed  to  stand  undisturbed  for  several  days. 
The  most  extreme  effects  were  seen  for  the  re-entrant  cavity 
tubes.  Figure  6.3  shows  the  change  in  behavior  (i.e.  for 
gassed  and  degassed)  for  a  Thermoexcel-E  tube  for  an 
increasing  heat  flux  run.  For  the  degassed  case,  an  extreme 
temperature  overshoot  was  observed,  which  indicates  a  large 
number  of  nucleation  sites  were  deactivated.  Once  all 
nucleation  sites  were  activated,  the  performance  of  both  cases 
was  identical.  It  can  be  surmised  that  the  entrapment  of  vapor 
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bubbles  in  surface  cavities  as  well  as  the  outgassing  from  the 

Teflon  O-rings  of  the  rig  over  time,  enhanced  the  nucleation 

on  the  tube  surfaces.  Vigorous  boiling  removed  those  gases 

from  both  surface  and  liquid. 

D.   BOILING  PERFORMANCE  OF  ENHANCED  TUBES  IN  PURE  R-114 
AND  R-114/OIL  MIXTURES 

The  results  of  the  74  data  runs  are  discussed  with 

references  to  the  figures   (Figures  6.1  to  6.80)  that  are 

provided  at  the  end  of  this  chapter.  Performance  of  each  tube 

will  be  presented  and  discussed  as  follows: 

1.  For  pure  R-114 

2.  For  R-114/oil  mixtures  (3%  and  10%  by  weight) 

3.  Comparison  of  pure  R-114  and  R-114/oil  mixtures  for 
increasing  heat  flux  runs. 

4.  Comparison  of  pure  R-114  and  R-114/oil  mixtures  for 
decreasing  heat  flux  runs. 

5.  Comparison  to  similar  tube  types  (i.e.  finned  -  as 
applicable  for  each  R-114/oil  mixture  for  both 
increasing  and  decreasing  heat  flux  runs) 

1 .  Boiling  Performance  of  Smooth  Tube 

Figures  6.4,  6.5  and  6.6  show  the  heat  transfer 
performance  of  pure  R-114,  R-114/3%  oil  mixture  and  R-114/10% 
oil  mixture  respectively,  boiling  from  a  smooth  copper  tube. 
Figures  6.7  and  6.8  show  the  comparative  performance  for  the 
three  refrigerant/oil  mixtures  on  increasing  and  decreasing 
heat  flux  respectively.  It  can  be  seen  that  in  Figure  6.7  in 
the  natural  convection  region  there  is  slight  degradation  with 
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increases  in  oil.  However,  in  the  nucleate  boiling  region 
there  is  significant  degradation  with  increasing  oil 
concentration.  Wanniarachi  et  al.  [Ref.  6]  reported  similar 
behavior  for  smooth  tubes.  This  behavior  was  attributed  to  two 
competing  effects:  (1)  the  ^choking'  of  the  tube  surface  by  an 
oil  rich  layer  which  develops  as  the  more  volatile  refrigerant 
is  boiled  off;  ()  the  diffusion  of  the  miscible  oil  rich  layer 
back  into  the  bulk  fluid.  In  the  natural  convection  region  it 
can  be  surmised  that  the  diffusion  effect  predominates  while 
in  the  nucleate  boiling  region  the  ^choking'  effect  eventually 
overwhelms  the  diffusion.  This  degradation  is  seen  again  in 
Figure  6.8.  It  can  be  seen  that  in  Figure  6.8  that  as 
conditions  approach  natural  convection  along  the  decreasing 
heat  flux  curve,  the  three  mixtures  perform  similarly  as 
nucleation  sites  die  out. 

2 .  Boiling  Performance  of  Finned/Modified  Finned  Tubes 
Figures  6.9,  6.10  and  6.11  show  the  heat  transfer 
performance  of  pure  R-114,  R-114/3%  oil  mixture  and  R-114/10% 
oil  mixture  respectively,  boiling  from  a  standard  26  fpi  tube, 
(GEWA-K)  .  At  a  heat  flux  of  35  kW/m2,  heat  transfer 
enhancement  compared  to  the  smooth  tube  for  pure  R-114, 
R-114/3%  oil  and  R-114/10%  oil  was  2.4,  2.9  and  3.2 
respectively.  Figures  6.12  and  6.13  show  the  comparative 
performance  of  the  three  ref rigerant /oil  mixtures  on 
increasing   heat   flux   and   decreasing   heat   flux   runs 
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respectively.  On  increasing  heat. flux  and  decreasing  heat  flux 
runs,  the  finned  tube  only  shows  a  small  degradation  in  both 
the  natural  convection  and  nucleate  boiling  regions.  The 
smooth  tube  only  shows  small  degradation  in  the  natural 
convection  region.  Also,  at  the  higher  heat  fluxes,  the 
R-114/3%  oil  mixture  shows  a  slightly  better  performance  than 
the  pure  R-114  and  R-114/10%  oil  mixture.  This  phenomenon  was 
observed  for  all  finned  tubes,  however,  the  reasons  for  such 
behavior  are  not  fully  understood;  it  may  be  due  to  bubble 
^scouring'  between  the  fins.  This  small  increase  in 
performance  at  low  oil  concentrations  agrees  with  results 
reported  by  Stephan  and  Mitrovic  [Ref .  15] . 

Figures  6.14,  6.15  and  6.16  show  the  heat  transfer 
performance  of  pure  R-114,  R-114/3%  oil  mixture  and  R-114/10% 
oil  mixture  respectively,  boiling  from  a  GEWA-K  40  fpi  tube. 
Figures  6.17  and  6.18  show  the  comparative  performance  of  the 
three  refrigerant/oil  mixtures  on  increasing  heat  flux  and 
decreasing  heat  flux  respectively.  Similar  behavior  to  the 
GEWA-K  2  6  fpi  tube  was  observed  including  the  improved 
performance  at  higher  fluxes  at  low  oil  concentrations  (3%) . 

Figures  6.19,  6.21,  and  6.23  show  the  comparative 
performance  of  the  three  R-114/oil  mixtures  boiling  from  GEWA- 
K  26  and  40  fpi  tubes  on  increasing  heat  flux.  For  all  three 
oil  concentrations,  the  GEWA-K  40  fpi  tube  outperforms  the 
GEWA-K  26  fpi  and  smooth  tubes  consistently  in  the  nucleate 
boiling  region.  At  a  heat  flux  of  35  kW/m2  the  GEWA-K  40  fpi 
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and  the  GEWA-K  26  fpi  had  heat  transfer  enhancements  of  4.1 
and  3.1  respectively  over  the  smooth  tube  for  pure  R-114. 
Similar  enhancements  were  seen  at  R-114/3%  oil  and  R-114/10% 
oil  mixtures.  See  Table  III  at  end  of  this  chapter.  In  the 
natural  convection  region,  the  results  were  mixed.  Neither 
tube  predominated  consistently  for  pure  R-114  and  R-114/oil 
mixtures.  Figures  6.20,  6.22  and  6.24  show  similar  comparisons 
on  decreasing  heat  flux.  It  can  be  seen  that  the  GEWA-K  4  0  fpi 
performs  best  in  all  three  cases  except  at  very  low  fluxes 
where  results  were  less  certain.  However,  the  calculation  for 
the  heat  flux  did  not  take  into  account  the  actual  increases 
in  surface  area  for  the  finned  tubes  over  the  smooth  tube.  A 
fin  root  diameter  was  used  to  calculate  surface  area.  The 
actual  area  increase  is  of  the  same  magnitude  as  the  heat 
transfer  enhancements  seen,  indicating  that  both  smooth  and 
finned  tubes  give  similar  heat  transfer  performance.  To 
demonstrate  this,  Figure  6.9,  the  GEWA-K  26  fpi  tube  in  pure 
R-114,  is  replotted  in  Figure  6.81  using  the  actual  fin 
surface  area.  The  ratio  of  actual  fin  surface  area  to  root 
diameter  area  was  a  factor  of  4.  Actual  surface  heat  flux  was 
calculated  by  dividing  the  root  diameter  surface  heat  flux  by 
four.  It  can  be  seen  in  Figure  6.81  that  in  the  natural 
convection  region  the  GEWA-K  2  6  fpi  tube  performs  comparably 
to  the  smooth  tube.  However,  in  the  mixed  boiling  and  nucleate 
boiling  regions  the  GEWA-K  2  6  fpi  tube  shows  an  enhanced 
performance  over  the  smooth  tube  up  to  three,  which  indicates 
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there  are  other  mechanisms  involved  besides  increased  surface 
area . 

Figures  6.25,  6.26  and  6.27  show  the  heat  transfer 
performance  of  pure  R-114,  R-114/3%  oil  mixture  and  R-114/10% 
oil  mixture  respectively,  boiling  from  a  modified  fin  tube, 
GEWA-T  19  fpi.  At  a  heat  flux  of  35  kW/m2,  the  heat  transfer 
enhancement  compared  to  the  smooth  tube  for  pure  R-114, 
R-114/3%  oil  mixture  and  R-114/10%  oil  mixture  was  3.1,  4.3 
and  5.2  respectively.  Figures  6.28  and  6.29  show  the 
comparative  performance  of  the  three  refrigerant/oil  mixtures 
on  increasing  and  decreasing  heat  fluxes  respectively.  Similar 
to  the  GEWA-K  type  tube,  only  a  slight  degradation  is  observed 
in  both  the  natural  convection  and  nucleate  boiling  regions. 
Again  at  higher  heat  fluxes,  the  R-114/3%  oil  mixture 
performed  best. 

Figures  6.30,  6.31  and  6.32  show  the  heat  transfer 
performance  of  pure  R-114,  R-114/3%  oil  mixture  and  R-114/10% 
oil  mixture  respectively,  boiling  from  a  GEWA-T  26  fpi  tube. 
At  a  heat  flux  of  35  kW/m2,  heat  transfer  enhancement  compared 
to  the  smooth  tube  for  pure  R-114,  R-114/3%  oil  mixture  and 
R-114/10%  oil  mixture  was  3.7,  4.6  and  4.8  respectively. 
Figures  6.33  and  6.34  show  the  comparative  performance  of  the 
three  refrigerant/oil  mixtures  on  increasing  and  decreasing 
oil  mixtures  respectively.  It  can  be  seen  that  at  higher  heat 
fluxes,  the  R-114/3%  oil  mixture  performed  best  and  there  was 
significant  degradation  in  performance  for  the  R-114/10%  oil 
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mixture  compared  to  the  smooth  tube.  At  lower  heat  fluxes,  the 
pure  R-114  outperformed  the  refrigerant/oil  mixtures.  It  can 
be  surmised  that  the  xchoking'  effect  predominates  for  the 
refrigerant/oil  mixtures. 

Figures  6.35,  6.37  and  6.39  show  the  comparative  heat 
transfer  performance  of  each  refrigerant/oil  mixture  boiling 
from  the  GEWA-T  19  and  26  fpi  tubes  for  increasing  heat 
fluxes.  Figures  6.36,  6.38  and  6.40  show  similar  comparisons 
for  decreasing  heat  fluxes.  It  can  be  seen  that  in  the 
nucleate  boiling  region  at  the  higher  heat  fluxes  that  the 
GEWA-T  26  fpi  tube  performed  best.  At  the  lower  heat  fluxes, 
in  the  natural  convection  region  for  an  increasing  heat  flux 
run  the  tubes  performed  equally.  For  a  decreasing  heat  flux 
run,  the  GEWA-T  19  fpi  performed  best.  Recalling  the  channel 
structure  of  the  GEWA-T  tube  suggests  more  complex  phenomena 
than  simple  increase  in  surface  area.  Indeed  Marto  et  al . 
[Ref.  16]  have  found  that  the  gap  width  with  GEWA-T  tubes  has 
a  significant  effect  on  boiling  performance.  Down  to  some 
limiting  value,  decreasing  the  gap  width  generally  increases 
heat  transfer  coefficient.  The  enhancement  can  be  attributed 
to  strong  interaction  between  heat  transfer  and  hydrodynamic 
effects  within  GEWA-T  channels. 

Figures  6.41,  6.42  and  6.43  show  the  heat  transfer 
performance  of  pure  R-114,  R-114/3%  oil  mixture  and  R-114/10% 
oil  mixture  respectively,  boiling  from  another  modified  fin 
tube,  GEWA-YX  26  fpi.  At  a  heat  flux  of  35  kW/m2,  heat 
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transfer  enhancement  compared  to  the  smooth  tube  for  pure 
R-114,  R-114/3%  oil  mixture  and  R-114/10%  oil  mixture  was  3.4, 
4.3  and  4.4  respectively.  Figures  6,44  and  6.45  show 
comparative  performance  for  the  three  refrigerant/oil  mixtures 
on  increasing  and  decreasing  heat  fluxes  respectively.  It  can 
be  seen  that  in  the  natural  convection  region  slight 
degradation  in  performance  occurs.  At  the  higher  fluxes  in 
the  nucleate  boiling  region,  the  R-114/3%  oil  mixture 
performed  best  as  seen  with  other  types  of  finned  tubes. 

Figures  6.46,  6.48  and  6.50  show  the  comparative 
performance  for  each  refrigerant/oil  mixture  boiling  from  the 
GEWA-K,  GEWA-T  and  GEWA-YX  2  6  fpi  tubes  for  increasing  heat 
fluxes.  Figures  6.47,  6.49  and  6.51  show  similar  comparisons 
for  decreasing  heat  fluxes.  Both  the  modified  finned  tubes, 
GEWA-T  and  GEWA-YX,  show  marked  increases  in  performance  over 
the  standard  finned  tube,  GEWA-K,  in  the  nucleate  boiling 
region.  Results  were  not  so  clear  however  in  the  natural 
convection  region,  possibly  due  to  variation  in  surface  area. 
3 .  Boiling  Performance  of  High  Flux  Tube 

Figures  6.52,  6.53  and  6.54  show  the  heat  transfer 
performance  of  pure  R-114,  R-114/3%  oil  mixture  and  R-114/10% 
oil  mixture  respectively,  boiling  from  the  High  Flux  tube.  At 
a  heat  flux  of  35  kW/m2,  heat  transfer  enhancement  compared  to 
the  smooth  tube  for  pure  R-114,  R-114/3%  oil  mixture  and 
R-114/10%  oil  mixture  was  6.4,  6.2  and  3.8.  The  High  Flux  tube 
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displays  a  much  larger  relative  temperature  overshoot  compared 
to  the  smooth  tube  and  its  point  of  incipience  occurs  at  a 
much  lower  heat  flux.  Accuracy  of  the  data  was  greatest  at 
high  heat  flux  and  large  superheat  values  (±  3.8%)  and  most 
uncertain  at  low  heat  flux  and  small  superheat  values 
(±  17.15%).  See  Uncertainty  Analysis-Appendix  E.  It  can  be 
seen  that  at  higher  heat  fluxes  that  the  High  Flux  tube  in  the 
R-114/10  oil  mixture  performed  almost  equally  to  the  smooth 
tube.  At  a  heat  flux  of  100  kW/m2  the  heat  transfer 
enhancement  compared  to  the  smooth  tube  was  1.1.  Wanniarachchi 
et  al.  [Ref.  7]  reported  similar  results.  The  marked  decrease 
in  boiling  performance  was  attributed  to  the  creation  of  an 
oil-rich  layer  within  the  boiling  pores.  The  superior  boiling 
performance  for  R-114/10%  oil  mixture  compared  to  the  smooth 
tube  is  completely  recovered  when  the  heat  flux  is  decreased. 
It  can  be  seen  that  in  the  natural  convection  region  the 
performance  of  the  High  Flux  tube  is  comparable  at  all  oil 
concentrations  to  the  smooth  tube  which  indicates  there  is  no 
effect  of  either  tube  type  or  oil  concentration  in  this 
region.  Figures  6.55  and  6.56  show  the  comparative  performance 
of  the  three  refrigerant/oil  mixtures  on  increasing  and 
decreasing  heat  fluxes  respectively.  The  extreme  rate  of 
degradation  of  performance  of  the  High  Flux  surface  at  10%  oil 
concentration  with  increasing  heat  flux  is  clearly  seen  as 
well  as  its  complete  recovery  as  the  heat  flux  is  decreased. 
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4 .  Boiling  Performance  of  the  Thermoexcel  Tubes 

Figures  6.57,  6.58  and  6.59  show  the  heat  transfer 
performance  of  pure  R-114,  R-114/3%  oil  mixture  and  R-114/10% 
oil  mixture  respectively,  boiling  from  the  Thermoexcel-E  Tube. 
At  a  heat  flux  of  35  kW/m2,  heat  transfer  enhancement  compared 
to  the  smooth  tube  for  pure  R-114,  R-114/3%  oil  mixture  and 
R-114/10%  oil  mixture  was  5.7,  5.3  and  5.2  respectively. 
Figures  6.60  and  6.61  show  the  comparative  performance  of  the 
three  refrigerant /oil  mixtures  on  increasing  and  decreasing 
heat  fluxes  respectively.  As  seen  with  the  High  Flux  tube,  the 
Thermoexcel-E  tube  performed  equally  to  the  smooth  tube  in  the 
natural  convection  region  for  all  refrigerant/oil  mixtures. 

Figures  6.62,  6.63  and  6.64  show  the  heat  transfer 
performance  of  pure  R-114,  R-114/3%  oil  mixture  and  R-114/10% 
oil  mixture  respectively,  boiling  from  the  Thermoexcel-HE 
tube.  At  a  heat  flux  of  35  kW/m2,  heat  transfer  enhancement 
compared  to  the  smooth  tube  for  pure  R-114,  R-114/3%  oil 
mixture  and  R-114/10%  oil  mixture  was  4.8,  5.1  and  4.8. 
Similar  performance  to  the  Thermoexcel-E  tube  was  observed. 
Figures  6.65  and  6.66  show  the  comparative  performance  of  the 
three  refrigerant/oil  mixtures  on  increasing  and  decreasing 
heat  flux  respectively.  At  a  higher  heat  flux  of  100  kW/m2  the 
heat  transfer  enhancements  compared  to  the  smooth  tube  for  the 
R-114/3%  oil  mixture  and  R-114/10%  oil  mixture  was  reduced  to 
2.2  and  1.9. 
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Figures  6.61,  6.69  and  6.71  show  the  comparative 
performance  of  each  refrigerant/oil  mixture  boiling  from  both 
the  Thermoexcel-E  and  HE  tubes  for  increasing  heat  fluxes. 
Figures  6.68,  6.70  and  6.72  show  similar  comparisons  on 
decreasing  heat  fluxes.  As  discussed  earlier  both  tubes  showed 
improved  performance  over  the  smooth  tube  in  the  nucleate 
boiling  region.  It  can  be  seen  that  the  Thermoexcel-HE  tube 
performed  marginally  better  than  the  Thermoexcel-E  at  high 
heat  fluxes.  In  the  natural-convection  region  both  tubes 
performed  comparably  to  the  smooth  tube.  Wanniarachchi  et  al 
[Ref.  8]  reported  similar  behavior  for  Thermoexcel-E  and 
Thermoexcel-HE  tubes.  It  can  be  surmised  that  in  the  natural 
convection  region  the  small  re-entrant  cavities  of  both  tubes 
do  not  increase  the  total  boiling  surface  area  like  fins.  The 
Thermoexcel-E  and  HE  tubes  ^appear'  smooth  and  would  be 
expected  to  perform  comparably  to  the  smooth  tube  which  was 
indeed  seen. 

5 .  Boiling  Performance  of  the  Turbo-B  Tube 

Figures  6.73,  6.74  and  6.75  show  heat  transfer 
performance  of  pure  R-114,  R-114/3%  oil  mixture  and  R-114/10% 
oil  mixture  respectively,  boiling  from  the  Turbo-B  tube.  At  a 
heat  flux  of  35  kW/m2,  heat  transfer  enhancement  compared  to 
the  smooth  tube  for  pure  R-114,  R-114/oil  mixture  and 
R-114/10%  oil  mixture  was  6.2,  6.7  and  5.6  respectively.  As 
seen  with  all  previous  re-entrant  tubes  tested,  comparable 
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performance  with  the  smooth  tube  is  observed  in  the  natural 
convection  region  and  marked  improved  performance  over  the 
smooth  tube  is  observed  in  the  nucleate  boiling  region. 
Figures  6.76  and  6.77  show  the  comparative  performance  of  each 
refrigerant /oil  mixtures  boiling  from  the  Turbo-B  tube  on 
increasing  and  decreasing  heat  fluxes  respectively.  Similar  to 
the  other  re-entrant  cavity  surfaces  tested,  the  R-114/10%  oil 
mixture  had  the  greatest  rate  of  degradation  for  increasing 
heat  flux  in  the  nucleate  boiling  region. 

6.  Overall  Performance  of  Smooth  and  Enhanced  Surfaces 
Figures  6.78,  6.79  and  6.80  show  the  comparative 
performance  of  8  of  the  10  tube  surfaces  tested.  The  other 
two,  Thermoexcel-E  and  GEWA-T  19  fpi,  performed  similarly  to 
the  Thermoexcel-HE  and  GEWA-T  26  fpi  respectively,  and  so  are 
not  included.  Three  distinct  groups  of  tubes  can  be  identified 
based  upon  boiling  heat  transfer  performance  at  all  oil 
concentrations.  These  are:  (1)  smooth,  (2)  finned/modified 
finned  and  (3)  re-entrant  cavity.  Within  the  finned/modified 
finned  group,  it  can  be  seen  that  the  GEWA-T  is  the  best 
performer.  Its  enhanced  performance  was  attributed  to  its 
channel  geometry.  The  High  Flux  tube  in  R-114/10%  oil  mixture 
performed  comparably  to  the  smooth  tube  (heat  transfer 
enhancement  was  1.1)  at  a  heat  flux  of  100  kW/m2.  The  severe 
degradation  in  performance  for  the  High  Flux  tube  can  be 
attributed  to  the  xchoking'  of  the  porous  surface.  Table  III 
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shows  numerically  the  enhancement  in  performance  over  the 
smooth  tube  for  all  three  refrigerant/oil  mixtures  boiling 
from  each  of  the  tubes  for  three  heat  fluxes  of  10,  35  and  100 
kW/m?.  It  can  be  seen  at  the  moderate  heat  flux  of  35  kW/m2, 
the  finned/modified  finned  tubes  show  enhanced  performances  of 
2.4  to  5.2  over  the  smooth  tube  while  the  re-entrant  cavity 
tubes  show  enhancements  of  3.8  to  6.7  over  the  smooth  tube. 
The  effect  of  oil  concentration  is  significant  at  high  heat 
fluxes.  At  a  heat  flux  of  100  kW/m?  for  pure  R-114,  the  finned 
and  re-entrant  cavity  tubes  show  comparable  enhancement 
compared  to  the  smooth  tube.  For  R-114/10%  oil  mixtures  at  a 
heat  flux  of  100  kW/m?,  the  finned  tubes  performed  best. 
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VII.  CONCLUSIONS 


1.  Accurate  and  repeatable  boiling  heat-transfer  data 
using  R-114  and  R-114/oil  mixtures  was  obtained  on  ten 
enhanced  tube  surfaces. 

2.  Three  distinct  groups  of  boiling  tubes  can  be 
identified  by  structure  and  boiling  performance: 
smooth,  finned  and  reentrant  cavity. 

3.  Significant  increase  in  performance  was  observed  for 
reentrant  cavity  surfaces  versus  either  the  finned  or 
smooth  surfaces  within  the  nucleating  regime  at  all 
heat  fluxes  tested. 

4.  Degradation  in  boiling  performance  accompanied 
increases  in  oil  concentration  within  the  nucleating 
regime.  Most  notably,  the  high  flux  surface  performed 
comparably  to  a  smooth  surface  at  the  highest  flux 
tested. 

5.  Within  the  convection  regime  all  tube  surfaces  tested 
performed  equally  well.  No  significant  effect  of  oil 
concentration  can  be  observed  up  to  incipient 
nucleation. 

6.  Distinct  temperature  overshoot  within  the  convection 
regime  occurs  repeatably  in  all  tube  surfaces  tested. 
The  reentrant  cavity  surfaces  exhibit  the  most  extreme 
overshoot  versus  either  the  finned  or  smooth  surfaces. 
Oil  concentration  has  no  significant  effect. 
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VIII.  RECOMMENDATIONS 


1.  Due  to  the  way  a  real  flooded  evaporator  operates,  it 
would  be  preferable  to  heat  the  tube  within  the  single 
tube  apparatus  with  hot  water. 

2.  As  a  complement  to  present  single  tube  and  bundle  tube 
studies,  the  present  single  tube  apparatus  should  be 
modified  to  accommodate  two  vertically  aligned  boiling 
tubes  in  order  to  study  the  effects  of  multiple  tubes 
and  variation  of  tube  pitch  to  diameter  ratios. 

3.  The  physical  properties  of  the  R-114/oil  mixtures 
tested  should  be  measured  for  comparison  to  the 
predictions  used  in  the  data  reduction  program  and  for 
better  understanding  of  the  mechanisms  driving  the 
degradation  in  heat  transfer  performance. 

4.  Neutrally  buoyant  particles  should  be  placed  in  the 
pool  to  facilitate  study  of  the  circulation  patterns 
within  the  flooded  evaporator  in  more  detail. 

5.  A  high  speed  camera  should  be  used  to  study  the 
nucleation  process  and  circulation  patterns  in  more 
detail . 
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APPENDIX  A:  THERMOPHYSICAL  PROPERTIES  OF  R-114 

The  following  thermophysical  properties  of  saturated  R-114 
are  plotted  versus  temperature  °C: 


1. 

Liquid 

Density  (kg/m3) 

2. 

Liquid 

Dynamic  Viscosity  (Pa*s) 

3. 

Liquid 

Thermal  Conductivity  (W/m*K) 

4. 

Liquid 

Specific  Heat  (J/kg*K) 

5. 

Latent 

Heat  (J/kg) 
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APPENDIX  B:  DATA  ACQUISITION  APPARATUS  CALIBRATION 

A.   BACKGROUND 

Early  data  runs  revealed  a  disparity  in  calculated  values 
between  previously  reported  theses  and  those  presently 
obtained  for  heat  flux  and  heat-transfer  coefficients.  The 
thermocouple  channels  were  checked  for  validity  of  temperature 
readings  using  other  programs  and  data  acquisition  equipment. 
Their  readings  were  found  to  be  as  accurate  as  possible 
(±  3%) .  Voltage  readings  originating  from  the  RMS  converter 
were  also  checked  out  as  accurate.  The  RMS  converter  converts 
true  voltage  measured  into  a  proportional  voltage  output  (0-10 
volts)  to  the  data  acquisition  unit.  The  proportionality 
constant  to  convert  RMS  voltage  readings,  Vs,  to  real  tube 
heater  voltage,  V,  was  equal  to  25.  This  was  verified  using  a 
standard  voltmeter  across  the  power  leads  of  the  cartridge 
heater  measured  over  several  power  settings.  The  discrepancy 
was  identified  to  be  the  measurement  of  the  current  supplied 
to  the  tube  cartridge  heater  by  the  VARIAC  unit  as  measured  by 
the  inductive  pickup  (Figure  3.13)  .  The  inductive  measurement 
of  current  (in  proportional  volts)  is  no  longer  run  through 
the  RMS  converter.  Instead,  it  is  run  directly  to  the  Data 
Acquisition  Unit.  A  procedure  to  measure  the  true  current  and 
calculate  the  proportionality  constant  for  converting  the 
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inductive  current,  Is,  to  the  true  current,  I,  was  established 
as  follows: 

B.  EQUIPMENT  USED 

1 .  Voltmeter 

2.  2  ohm  resistor 

3.  HP-3852A  Data  Acquisition  Unit 

4.  Computer  -  HP-9300  Series 

5.  Current  Sensor  (Inductive  Pickup) 

6.  Voltage  Sensor /RMS  Converter 

7.  Cartridge  Heater  -  lkW,  220V  rating 

8.  VARIAC  -  power  supply  220V  AC 

9.  SETUP8  program 

C.  CALIBRATION  PROCEDURE 

1.  Ensure  all  power  is  turned  off  to  the  VARIAC. 

2.  Insert  the  standard  2  ohm  resistor  in  series  with  the 
cartridge  heater  as  shown  in  Figure  C.l. 

3.  Connect  leads  of  voltmeter  across  the  2  ohm  resistor  to 
measure  the  drop  in  voltage.  At  any  given  power 
setting,  the  current  through  the  resistor  will  equal 
the  current  through  the  cartridge  heater.  This  assumes 
the  resistances  of  both  heater  and  resistor  are 
constant.  This  was  verified  as  shown  in  Figure  C.2.  The 
total  resistance  of  the  heater/resistor  rapidly  and 
asymptotically  approached  a  constant  value. 

4.  Connect  power  to  the  VARIAC. 

5.  Load  and  run  SETUP8  program  to  provide  the  voltage  and 
current  readings  made  by  the  data  acquisition  unit. 
Note  the  SETUP8  program  used  in  calibration  differed 
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from  that  listed  in  Appendix  G.  A  proportionality 
constant  of  2.0  was  assumed  for  calibration  then 
modified  once  the  true  value  was  calculated. 

6.  Starting  with  the  VARIAC  set  at  zero  volts, 
incrementally  increase  the  voltage  up  to  the  maximum  of 
220V. 

7.  At  each  voltage,  measure  the  drop  in  voltage  across  the 
resistor.  The  current  (amps)  is  calculated: 

I  =  Drop  in  voltage  (volts)  /  2  ohms 

8.  Note  and  compare  the  current  provided  by  the  program. 

9.  At  the  maximum  220V,  incrementally  decrease  voltage 
down  to  0,  taking  comparative  current  readings  as 
before,  in  order  to  detect  any  hysteresis  effect.  None 
was  observed. 

10.  Plot  current  (amps)  vs.  voltage  (volts)  .  See  Figure 
C.3.  Note  the  linearity  of  both  data  sets.  The  true 
meter  reading  is  some  factor  smaller  than  the  data 
acquisition  readings. 

11.  The  calibration  constant, C,  is  calculated  by: 
C  =  (I  Meter  current) *  2 

(E  Data  acquisition  current) 
C  =  1.9182 


D.   DISCUSSION 

A  disparity  between  the  previous  reported  theses  values 
and  those  presently  obtained  still  exists.  Calculated  values 
reported  in  previous  theses  for  the  heat  flux  and  the  heat 
transfer  coefficient  are  consistently  higher.  No  previous 
thesis  reports  such  a  calibration  as  carried  out  above  and  it 
is  thought  that  this  is  the  cause  of  the  discrepancy. 
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APPENDIX  C:  AN  EXAMPLE  OF  REPRESENTATIVE  DATA  RUN 


Month,  date  and  tine  :  5  Jul  1991        17:21:51 

NOTE:  Prooran    name    :  DRP8 
Da  si  nunber   •  01 
New  file  name:    DATA070SI5 
TC  is  defective  at  location   2 
Tube  Nunber:    5 

Data  Set  Nunber  -    I   Bulk  Oil  X   •   0.0  2.11545538805 

TC  No:     I  2      3  4  5      £      7       8 

Tenp  :  5.11  0.00   5.15  5.01    5.12   5.08   5.09   5.11 

Twa  TjjQd  Tl.tqd2   Tvapr  Psat    Taunp 

5.08  2. 36  2. 26    3.91  -1.65   -13.8 

Thetab  Htube      Odp 

2.770  1. 590E+02   4.406E+02 

Data  Set  Nunber  -   2  Bulk  Oil  X  -  0.0      2.11545539033 
TC  No:     1      2      3      4      5      £      7      8 

Tenp  :  4.94    0.00   5.04    4.87  4.99    4.95    4.96    4.98 

Twa  Tllqd   Tliqd2   Tvapr  Peat    Ttunp 

4.95  2.23    2.18    3.29  -1.76  -13.7 

Thetab  Htube      Odp 

2.740  I.603E+02  4.391E+02 

Data  Set  Nunber  -   3  Bulk  Oil  X  -  0.0      2. It 545539756 

TC  No:     1      2      3      4  5      6      7      8 

Tenp  :    6.08   0.00   6.19   5.96  6.14   6.09   6.11    6.15 

Twa     Tliqd    Tliqd2   Tvapr   Paet  Tsunp 

6.07   2.38    2.16    3.27  -1.69  -13.6 
Thetab   Htube      Odp 
3.806   1.E94E+02  6.446E+02 

Dete  Set  Nunber  -    4   Bulk  Oil  X  -   0.0      2.11545539780 
TC  No:     I       2      3      4       5      6      7      8 
Tenp  :    6.12   0.00   6.23   6.02   6.18   6.14   6.15   6.18 
Twe     Tllqd    Tllqd2   Tvepr   Paet    Tsunp 

6.12  2.39    2.19    3.29  -1.67   -13.6 
Thetab   Htube      Odp 

3.828   1.682E+02   6.439E+02 

Data  Sat  Nunber  •   5  Bulk  Oil  X  -  0.0      2.11545540628 
TC  No:     1      2      3      4      5      6      7      8 
Tenp  :    7.13   0.00   7.27   6.97   7.24   7.15   7.17   7.25 
Twe    TJiqd   Tllqd2  Tvapr  Psat    Tsunp 

7.13  2.24    2.16    2.92  -1.77  -13.6 
Thetab  Htube      Odp 

4.934   1.791E+02   8.834E+82 

Dete  Set  Nunber  •   6  Bulk  Oil  X  -  0.0      2.11545540655 
TC  No:     1      2      3      4      S      6      7      8 

Tenp  :    7.09   0.00   7.25   6.95  7.21    7.13   7.15   7.20 
Twe     Tllqd    Tliqd2   Tvepr   Pset    Tsunp 
7.11   2.20    2.12    2.91   -I. 81   -13.5 
Theteb  Htube      Odp 
4.945   1.786E+02   8.834E+02 

Oete  Set  Nunber  -   7  Bulk  Oil  X  •  0.0  2.11545541023 

TC  No:     1       2       3  4       5       6       7       8 

Tenp  :    9.56   0.00   9.74  9.36   9.71    9.60   9.62   9.75 

Twe    Tliqd   Tliqd2  Tvepr  Peet    Tsunp 

9.57   2.47    2.22    2.94  -1.61   -13.5 
Thetab  Htube     Odp 
7.223   I.863E+02   1.346E+03 
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Data  Set  Number  •    8   Bull  Oil  I  -   0.0  2. I  154554 II GS 

TC  No:     I  2  3  4      5      E      7      8 

Tenp  :  9.33  0.00    9.52  9.11    9.49    9.37    9.39    9.50 

Twa  Tliqd  Tliod2   Tvapr  Peat    Taunp 

9.33  2.2)  2.13    2.82  -1.80  -13.5 

Thetab  Htube  Qdp 

7.IG5  1.874E+02   1.343E+03 

Data  Set  Nunber  -    9   Bulk  Oil  X  -  0.0  2 .  1  154S54154S 

TC  No:     I       2      3      4  5  E      7      8 

Tenp  :  10.99   0.00   11.53   11.09  11.50   11.38   11.40   11.40 

Twe  Tliod    TliodZ   Tvapr   Peat  Taunp 

11. 2E  2.30    2. IB    2.GB   -1.73  -13.5 

Thetab  Htube      Odp 
9.02E   1.963E+02   1.772E+03 

Data  Set  Nunber  -   10  Bulk  Oil  X  -  0.0      2.11545541589 
TC  No:     I       2      3      4      5      E      7      8 

Tenp  :   11.00   0.00   11.52   11.09  11.50   11.37   11.40   11.41 

Twa     Tliqd    Tliod2   Tvapr   Paat    Taunp 

11. 2E   2.32    2.22    2.E7  -I.B9  -13.5 

Thetab  Htube      Odp 
8.988   1.977E+02   1.777E+03 

Data  Set  Nunber  •   11   Bulk  Oil  X  •  0.0      2.11545541889 

TC  No:     1      2      3      4  5      E      7      8 

Tenp  :   11.45    0.00   12.32   11.93  12.30   12.18   12.21   12.10 

Twa     TllQd    TliooT   Tvapr   Paat  Taunp 

11.99   2.24    2.17    2.B7  -1.76   -13.5 

Thetab   Htube      Odp 
9.783  2.074E+02  2.029E+03 

Data  Set  Nunber  •   12   Bulk  Oil  X  -   0.0      2. 1 1S4S541 913 
TC  No:     I      2      3      4      5      E      7      8 
Tanp  :   11.50   0.00   12.38   11.97   12.35   12.24   12.25   12.14 

Twa     TJiod    Tliod2   Tvapr   Paat    Taunp 

12.04   2.2B    2.19    2.B9  -1.74   -13.5 

Thetab  Htube      Odp 
9.814  2.0SSE+02  2.027E+03 

Data  Sat  Nunber  •   13  Bulk  Oil  X  -   0.0      2.11545542290 
TC  No:     1      2      3      4      5      E      7      8 
Tenp  :   12.08    0.00   13.1)   12.71   13.09   12.98   12.99   12.82 

Twa     Tligd    Tlio,d2   Tvapr   Paat    Taunp 

12.74   2.22    2.20    2.96   -1.75  -13.5 

Thetab  Htube      Odp 

10.527  2.I12E+02  2.223E+03 

Data  Set  Nunber  •   14   Bulk  Oil  X  -   0.0      2.115455423)5 
TC  No:     I       2      3      4      5      6      7      8 
Tanp  :   12.02   0.00   13.05   12.67   13.03   12.92   12.94   12.75 

Twe     Tliod    TllQd2   Tvepr   Paat    Taunp 

12.68   2.20    2.17    2.96   -1.78   -13.5 

Thetab  Htube      Odp 

10.494   2.110E+02   2.214E+03 

Data  Set  Nunber  •   15  Bulk  Oil  X  -   0.0      2.11545542511 
TC  No:     1      2      3      4      5      B      7      8 
Tenp  :   12.55   0.00   13.63   13.53   13.82   13.72   13.73   13.43 

Twe     Tliod    Tliod2   Tvapr   Paat    Taunp 

13.42   2.19    2. IE    2.92   -1.79  -13. S 

Theteb   Htube      Odp 

11.239  2.I82E+02  2.452E+03 
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Data  Set  Number  -   IS   Bull'  Oil  X  -  0.0      2.11545542534 
TC  No:     I       2      3      4       5      S      7       8 

Teno  :   12.S2   0.00   13.85   13.61  13.83   13.75   13.75   13.47 

Twe     Tliqd    Tiiod?   Tveor   Psat    Tauno 

13. 4S    2.26     2.22     2.91   -1.72   -13.5 

The  tat.   Htube       Odo 

11.217   2.183E+02   2.449E+03 

Data  Set  Nunber  -   17   Bulk  Oil  Z  -  0.0      2.11545543206 

TC  No;     I       2      3      4  5      6      7       8 

Temp  :    2.89   0.00   2.74   2.69  2.69   2.59   2.57   2.68 
Twe     Tliqd    TJiod2   Tvepr   Psat    Taunp 

2.58   2.44    2.19    2.01   -1.64  -13.5 
The  tab   Htube      Odo 

.270   1.126E+04   3.036E+03 

Data  Set  Nunber  -   18   Bulk  Oil  X  -  0.0      2. 1 154S543229 

TC  No:     I       2      3      4  5      C      7      8 

Teno  :    2.86   0.00   2.71    2.66  2.67   2.56   2.55   2.65 

Twa     Tiiod    Tllqd2   Tvapr   Paat  Taunp 

2.56   2.39    2.15    1.98  -1.69  -13.5 
Thetab   Htube      Odp 

.288   I.0S2E+04   3.029E+03 

Data  Sat  Nunber  -   19  Bulk  Oil  S  -  0.0      2.11545543435 

TC  No:     I      2      3      4  5      6      7      8 

Tenp  :    3.15   0.00   2.97   2.92  2.92   2.79   2.75   2.91 

Twa    Tiiod   Tliqd2   Tvapr  Piat  Taunp 

2.76   2.41     2.17    2.09  -1.B7  -13.5 
Thetab  Htube      Odp 

.472   9.S76E+03  4.S23E+03 

Data  Set  Nunber  •   20  Bulk  Oil  X  -  0.0      2.11545543508 

TC  No:     1       2      3      4  S      6      7       8 

Tenp  :    3.12   0.00   2.93   2.89  2.89   2.75   2.72   2.88 

Twe     Tliqd    Tliod2   Tvapr   Paat  Taunp 

2.73   2.47    2.19    2.08  -1.62  -13.5 
Thetab  Htube      Odp 

.394   I.I44E+04   4.S06E+03 

Date  Set  Nunber  -   21   Bulk  Oil  X  -  0.0      2.11545543579 

TC  No:     1      2      3      4  5      E      7      8 

Tenp  :    3.02   0,00   2.83   2.79  2.79   2.B4   2.62   2.77 

Twa     Tliqd    Tiiod?   Tvapr   Paat  Taunp 

2.62   2.37    2.09    1.98   -1.74  -13.5 
Thetab  Htube      Odp 

.394   1.I49E+04   4.S2SE+03 

Data  Set  Nunber  -  22  Bulk  Oi 1  X  -  0.0      2.11545543829 

TC  Mo:     1       2      3      4  5      E      7      8 

Tenp  :   3.45   0.00   3.20   3. IE  3.15   2.9B   2.91   3.14 

Twa    Tliqd   Tliqd?   Tvapr  Paat  Taunp 

2.89   2.41    2. IB    2.7E  -1.E7  -13.5 
Theteb   Htube      Odp 

.E03   I.178E+04  7.099E+03 

Data  Sat  Nunber  -   23   Bulk  Oil  X  -  0.0      2.11545543855 

TC  No:     I      2      3      4  5      E      7      8 

Tenp  :    3.44   0.00   3.19   3.15  3.14   2.9E   2.90   3.13 

Twa     Tiiod    Tliqd2   Tvapr   Paat  Taunp 

2.88   2.39    2.15    2.87   -1.E9  -13.5 
Theteb   Htube      Odp 

.6  1?   I.15BE+04   7.092E+03 
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Data  Set  Number  •   24   Bull  Oil  X  -  0.0      2.11545944709 
TC  No:     1      2      3      4      5      6      7      8 

Tenp  :    3.95   0.00   3.S6   3.55  3.50   3.29   3.20   3.52 
Twe     TJiod    TJiod2   Tveor   Pest    Taunp 
3.12    2.36     2.17    2.71   -1.69   -13.6 
Theteb   Htube      Qdo 

.857   I.279E+04   I.09EE+04 

Date  Set  Nunber  -   25   Bull  Oil  Z  -  0.0      2.11545544732 
TC  No:     I       2      3      4      S      6      7      8 

Tenp  :    3.95    0.00    3.56    3.55  3.49    3.28    3.19   3.52 
Twe     Tljod    TliodT   Tvapr   Paat    Taunp 
3.12   2.36    2.16    2.75   -1.70   -13. E 
Thetab   Htube      Qdo 

.859   1.279E+04   1.09BE+04 

Data  Set  Nunber  •   26   Bulk  Oil  Z  •  0.0      2.11545545094 
TC  No:     1       2      3      4       5      6      7      8 

Tenp  :    4.48    0.00   3.92    3.95  3.85    3. 62    3.49    3. 92 
Twa     TJiod    Tliqd2   Tvapr   Paat    Taunp 
3.36   2.38    2.19    2.87  -1.67   -13.6 
Thetab   Htube      Qdo 
1.073   1.4116*04   1.S14E+04 

Data  Set  Nunber  -   27  Bulk  Oil  Z  •  0.0      2.11545545119 
TC  No:     I       2      3      4       5      6      7      8 

Tenp  :    4.45   0.00   3.89   3.93  3.83    3.60    3.46    3.89 
Twa     Tliad    Tliad2   I.ior   Paat    Taunp 
3.33   2.33     2.15    2.04   -1.72   -13.6 
Thetab   Htube      Odo 
1.092   1.38SE+04   1.S12E+04 

Data  Set  Nunber  -   28   Bulk  Oil  Z  -  0.0      2. 1 I5455454E0 

TC  No:     1       2      3      4  5      6      7      8 

Tenp  :    5.12    0.00   4.37   4.43  4.27    4.04    3.86    4.40 
Twa     Tiiod    TliQd2   Tvapr   Paat    Taunp 

3.65    2.35     2.17    2.08   -1.70  -13.6 
Thetab   Htube      Odo 
1.390   1.4376*04   1.996E+04 

Deta  Set  Nunber  -   29  Bulk  Oil  Z  -  0.0      2.11545545484 

TC  No:     I       2      3      4  5      6      7      8 

Tenp  :    5.12    0.00   4.3E   4.42  4.26    4.03    3.84    4.39 

Twa     Tliod    TliQd2   Tvapr   Paat  Taunp 

3.64    2.34     2.16     2.07   -1.71  -13.6 
Thetab   Htube      Odp 
1.393   I.433E+04   1.99SE+04 

Data  Set  Nunber  -   30  Bulk  Oil  Z  -  0.0      2.11545546220 

TC  No:     I      2      3      4  5      6      7      8 

Tanp  :    6.47   0.00   5.25    5.39  5.11    4.90   4.58    5.33 

Twa     Tliad    Tliod2   Tvapr   Paat  Taunp 

4.23   2.31    2.18    2.10  -1.72  -13.6 
Thetab   Htube      Odp 
1.989   1.505E+04   2.993E+04 

Data  Set  Nunber  -   31   Bulk  Oil  Z  •  0.0      2.11545546248 

TCNo:l       2      3      4  5      6      7      8 

Tano  :    6.45   0.00   5.23    5.37  5.09   4.88   4.57    5.31 

Twe     Iliad    Tliad2   Tvapr   Paat  Taunp 

4.22    2.27     2.15     2.08   -1.75  -13.6 
Thetab   Htube       Odp 
2.003   1.497E+04   2.998E+04 
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Data  Set  Number  -   32   Bull  Oil  X  -  0.0      2.  I  IS4SS4EE59 
TC  No:     I       2      3      4       5      G      7      8 

Temp  :    8.99   £.00    G.6E    £.92  C.39   G.17    5.57    G.G3 

Twe    Tliqd   Tliod2   Tvapr   Paet    Tsunp 
5.16   2.33    2.18    2.13   -1.71   -13.4 

Thetab  Htube      Odp 

2.910  I.SS9E+04      4.S3SE+04 

Data  Set  Number  -   33   Bulk  Oil  I    -  0.0      2. 1 1545546682 

TC  No:     1       2      3      4  5      6      7      8 

Tenp  :  8.97   0.00   G.G4   6.90  5.37   6.14   5. 56   6.62 

Tua  Tli ad    TlioaT   Tvapr   Paat  Tiunp 

5.15  2.30    2.17    2.10   -1.73  -13.3 

Thetab  Htube      Qdp 

2.911  1.SS6E+04  4.S31E+04 

Data  Set  Number  -   34   Bulk  Oil  2  -   0.0  2.11545546951 

TC  No:     I      2      3  4  5  6      7      8 

Tenp  :  12.55  0.00   8.98  9.39   8.54   8.06   6.96   8.56 

Tua  Tliod  Tliodl'   Tvapr  Paat    Taump 

6.53   2.29    2.16    2.15  -1.74   -13.2 
Thetab  Htube      Odp 
4.305   1.E32E+04   7.028E+04 

Data  Sat  Number  ■   35   Bulk  Oil  I  -   0.0  2.11545546986 

TC  No:     1      2      3  4  5      6      7      8 

Tamp  :  12.59   0.00   9.01  9.42   8.57   8.09   6.99   8.60 

Twe  Tliod    Tiled;   Tvapr  Paat    Taump 

6.S6  2.31    2.19    2.18  -1.71   -13.2 

Thatab  Htube      Odp 

4.. 3)5  1.627E+04   7.0196+04 

Data  Set  Number  -  36  Bulk  Oil  J  ■  0.0      2.11545547257 

TC  No:     I      2      3      4  5      6      7      8 

Temp  :   16.06   0.00   11.35   11.79  10.57   10.04   8.26   10.43 

Twe     Tliqd    Tliqd2   Tvepr   Paet  Taump 
7.84   2.23    2.15    2.14   -1.77   -12.9 
Thetab  Htube      Odp 
5.E49   1.694E+04   9.S69E+04 

Oete  Set  Number  •  37  Bulk  Oil  Z  -  0.0      2.11545547281 
TC  No:     I      2      3      4      5      6      7      8 
Tamp  :   16.11    0.00   11.39   11.81   10.59   10.07   8.28   10.46 
Twe     Tllqd    Tllqd2   Tvepr   Paet    Taump 
7.87   2.26    2.18    2.17  -1.75   -12.8 
Thetab  Htube      Qdo 
5.656   1.G92E+04   9.572E+04 
NOTE:  37  dete  rum  were  atorad  in  file  OATA0705IS 

NOTE:  37  X-Y  palra  were  atorad  in  plot  data  file  PLOT0705I5 
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APPENDIX  D:  SAMPLE  CALCULATION 

Data  run  number  9  (saturation  temperature  was  2.22  °C  and 
heat  flux  was  20.14  kW/m2)  was  chosen  for  the  sample 
calculation . 

A.   TEST-SECTION  DIMENSIONS 


Do  = 

0 

.01588 

m 

Di  = 

0 

.01270 

m 

Dl  = 

0, 

.01143 

m 

D2  = 

0, 

.01588 

m 

L 

0, 

.20320 

m 

Lu  = 

0. 

.07620 

m 

B.   MEASURED  PARAMETERS 

V  =  112.08  volts 

I  =  1.94  volts 

Tl  =  16.51  °C 

T2  =  15.80  °C 

T3  =  15.88  °C 

T4  =  15.93  °C 

T5  not  read  due  to  defective  thermocouple 

T6  =  15.94  °C 

T7  =  15.63  "C 

T8  =  16.16  °C 

Tsat  =  2.22  °C 

kc  =  344.00  (W/m-K) 


C.   OUTER  WALL  TEMPERATURE  OF  THE  BOILING  TUBE 

p   =  ij-Do  =  3.1416  ■  0.01588  =  0.0499  (m) 

Ac  =  7i  (Do2-Di2) /4  =  3.1416  ■   [  (0  .  01588) 2- (0  .  0127)  2] /4 

Ac  =  7.14E-5  (m2) 

Qh  =  VI  =  112.08  •  1.94  =  217.44  (W) 

Tavg  =   £  Tn  /n  =  15.97  (°C) 

Two  =   Tavg  -  Qh-  [  ln(D2/Dl)  /  (2-n-L-kc)  ] 

Two  =   15.97  -  217.44- 

[In  (15.  88/11. 43) / (2-  3.1416- 0.203-344) ] 
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Two 

e 


15.80  CC) 

Two  -  Tsat  =  15.80  -  2.22  =  13.58  °C 


PROPERTIES  OF  R-114  AT  FILM  TEMPERATURE 


Tf 
H 

Tc 
Tf 

J 
P 

P 
v 

k 
k 

Tf 
Cp 

Cp 
a 

P 

Pr 


(Two+Tsat)/2  =  (15.80+2.22) /2  =  9.01 
exp[-4.4636  +  (1011 . 47/Tf ) ]  •  1E-3 
exp[-4. 4636+ (1011.47/9.01)  ]   •   1E-3 
(N-  s/m2) 

Critical  Temperature  (R)  =  753.95  (R) 
9.01  CC)  =  507.89  (R) 
1  -  Tf(R)/Tc(R)  =  1  -  507.89/753.95 


4.16E-4 


581.77  +  984.15-  j(1/3)  +  263.02-j  + 
(kg/m3) 


=  0.33 
279.99' 


1/2 


17.94- j 

1504.22  (kg/m3) 

|i/p  =  4.16E-4/1504 

7.10E-2  -  (2.61E-4 

7.10E-2  -  (2.61E-4 

9.01  CC)  =  282.16 

400  +  (1.65  •  Tf) 

Tf3) 

970.43  (J/kg-K) 

k/p-Cp  =  6. 86E-2/ (1504.22- 970 

-(Ap/AT)/p  =  1.92E-3  (1/K) 

v/o  =  5.88 


22  =  2.77E-7 

•  TfCC)) 

•  9.01)  =  6.86E-2 
(K) 

+  (1.51E-3  •  Tf2) 


(W/m-  K) 
-  (6.68E-7 


43)  =  4.7E-8  (m2/s) 


E.   HEAT-FLUX  CALCULATION 

Average  natural-convection  heat  transfer  coefficient  at 
non-boiling  ends: 


h  =  k_ 

Do 


0.60+0.387 


g*  ft-  (Do3)  ;  8-  Tanh  (m-  Lu) 
v  a-  Lu*  m 

[l+(0.559/Pr)9/16]8/27 
m   =    [h-p/  (kc-Ac)  ]1/2 


1/6 


h  =  235.4  (W/m2-K) 


The  value  for  h  was  computed  iteratively.  A  value  for  h  of  190 
(W/m2-K)  was  assumed  at  start  of  iteration. 


Heat-transfer  rate  through  non-boiling  ends: 


m 
Qf 


21.87  (1/m) 

(h-p- kc-Ac)  1/2- 8- Tanh  (m- Lu)  =  6.79  (W) 
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Heat  flux  through  active  boiling  surface 

Q  =  Qh  -  2-Qf  =  217.44  -  2  ■  6.79  =  203.85  (W) 

Ab  =  it-Do-L  -  3.1416  •  0.015875  •  0.2032  =  1.01E-2  (m2) 

q"  =  Q/Ab  =  203. 85/1. 02E-2  =  20,108  (W/m2) 

h  =  q/0  =  20,108/13.59  =  1480  (W/m2-K) 

The  following  are  the  results  obtained  from  the  computer 
by  running  the  data  reduction  program  DRP8  (See  Appendix  H) . 

q"  =  20,  140  (W/m2) 
0   =  13.64  CO 
h   =  1477  (W/m2-K) 
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TABLE  IV.  DIMENSIONS  OF  BOILING  TUBES 


Dl 
(mm) 

D2 
(mm) 

Di 
(mm) 

Do 
(mm) 

L 
(mm) 

Lu 
(mm) 

Kcu 
(W) 

m2-  K 

Smooth 

12.44 

15.88 

12.70 

15.88 

203.2 

76.20 

344 

GEWA-K 
26  fpi 

10.1 

12.7 

11.11 

12.7 

203.2 

76.2 

344 

GEWA-K 
40  fpi 

10.1 

12.7 

11.11 

12.7 

203.2 

76.2 

344 

GEWA-T 
19  fpi 

11.57 

13.8 

11.8 

13.31 

203.2 

76.2 

398 

GEWA-T 
26  fpi 

11.57 

13.8 

11.8 

13.31 

203.2 

76.2 

398 

GEWA-YX 
2  6  fpi 

11.57 

13.8 

11.8 

13.31 

203.2 

76.2 

398 

Thermoexcel- 
E 

11.57 

13.8 

11.8 

13.31 

203.2 

76.2 

398 

Thermoexcel- 
HE 

11.57 

13.8 

11.8 

13.31 

203.2 

76.2 

398 

Turbo-B 

11.57 

13.8 

11.8 

15.8 

203.2 

76.2 

398 

High  Flux 
(Copper/ 
Nickel 
95/5) 

12.95 

15.82 

13.2 

15.82 

203.2 

76.2 

45 
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APPENDIX  E:  UNCERTAINTY  ANALYSES 

Four  data  points  from  two  data  runs  were  chosen  for  the 
uncertainty  analysis:  data  run  number  8  (smooth  tube,  0%  oil, 
100,000  &  3500  W/m2)  and  data  run  number  74  (Thermoexcel-E,  0% 
oil,  6000  &  100,000  W/m2.  The  dimensions  of  the  boiling  tubes 
can  be  found  in  Appendix  D,  Table  IV.  The  measured  and 
calculated  parameters  used  in  this  analysis  were  obtained  as 
shown  in  Appendix  D  and  are  listed  in  Table  V  of  this 
appendix.  All  uncertainties  are  presented  as  a  percentage  of 
the  calculated  parameter.  What  follows  is  a  sample  calculation 
for  data  run  number  8:  smooth  tube,  0%  oil  at  a  heat  flux  of 
93,130  W/m2  using  the  methods  outlined  by  Kline  and  McClintock 
[Ref .  17] .  All  other  data  point  uncertainties  were  calculated 
similarly,  the  results  of  which  are  listed  in  Table  V  of  this 
appendix . 

A.   UNCERTAINTY  IN  SOURCE  HEAT-TRANSFER  RATE 

Qh  =  VI  W 

Is  =  2.122  volts       61   =  ±  0.025  amp 

Vs  =  9.45  volts        6V   =  ±  0.05  volts 

I  =  1.9182- Is  =  4.07  amp 

V  =  25-Vs  =  236.26  volts 

where : 

6   =   uncertainty  in  measurement  and  calculation 

6Qh/Qh  =   ((6V/Vs)2  +  (6l/Is)2)1/2 

6Qh/Qh  =    ( (0.05/9. 45)2  +  (0 . 025/2  .  122) 2)  K2 

6Qh/Qh  =   1.29  percent 
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B.   UNCERTAINTY  IN  SURFACE  AREA 


Ab  =  n-Do-L 

Do  =  15.88  (mm)        6Do  =    0.1  (mm) 

L   =  203.20  (mm)       6L   =    0.1  (mm) 

6Ab/Ab  =   ((5Do/Do)2  +  (6L/L)2)1/2 

6Ab/Ab  =   (  (0.1/15. 88)2  +  (0  . 1/203  .  2) 2)  1/2 

6Ab/Ab  =   0.63  percent 


C.  UNCERTAINTY  IN  WALL  SUPERHEAT 

AT  =   Two  -  Tsat 

Tsat    =   2.22  °C       6T   =    0.01  "C 

Two  =   Tavg  -  Qh [  (In (D2/D1) ) / (2-  ir  L- kc) ] 

Tavg  =   (ETn/8)  where  n  =  1  to  8 

Tn  =  thermocouple  readings 

Tl  =  23.86  °C         T2   =  19.62  °C     T3  =  20.38  T 

T4  =  19.87  °C         T5   =  Defective  thermocouple 

T6  =  20.18  °C          T7   =  20.16  °C      T8  =  20.20  °C 

Tavg  =   2  0.18  °C 

S.D.    =   (  (I  (Tn  -  Tavg)2) /7)1/2  =  0.19  °C 
where  S.D.  =  standard  deviation 

The  logarithmic  term  in  equation  of  Two  is  very  small 
compared  to  the  standard  deviation,  this  term  can  be  neglected 
for  the  uncertainty  analysis. 

Two     =   Tavg  =  18.72  °C       6Two  =  S.D.  =  0.19  °C 

AT      =   17.96  °C 

6AT/AT  =   ((6Two/AT)2  +  (-6Tsat/AT)2) 1/2 
6AT/AT  =   (  (0.93/16. 5)2  +  (0  .  01/16  .  5) 2)  1/2 
6AT/AT  =   1.06  percent 

D.  UNCERTAINTY  IN  HEAT  FLUX 

q   =    (Qh-2Qf)/Ab 

Qh  =   961.58  W  6Qh  =    12.42  W 
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Assuming  the  same  proportion  in  the  uncertainty  for  Qf : 

Qf  =   8.99  W  5Qf  =    0.12  W 

Qh  -  2Qf   =   943.6  W 

6q/q    =    [  (5Qh/ (Qh-2Qf) )2  +     (2dQf / (Qh-2Qf ) )  2)     + 

(SAb/Ab)2)  ]1/2 

6q/q    =    [  (12. 42/943. 6)2  +  (0  .  24/943  .  6)  2  +  (0  .  0063)  2]  1/2 
6q/q    =   1.46  percent 


UNCERTAINTY  IN  BOILING  HEAT-TRANSFER  COEFFICIENT 

h   =   q/AT 

6h/h        =        [(&q/q)2    +    (6AT/AT) 2] 1/2 
6h/h         =        [(0.0146)2    +     (0.01.06)2]1/2 
6h/h        =       1.8   percent 
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TABLE  V 

.  UNCERTAINTY 

ANALYSIS  OF  FOUR  DATA  POINTS 

PARAMETERS 

SMOOTH 
100  kW/m2 

SMOOTH 
3500  W/m2 

THERMO-E 
6000  W/m2 

THERMO-E 
100  kW/m? 

Is  (amp) 

2.122 

0.454 

0.494 

2.06 

Vs  (volts) 

9.45 

2.09 

2.20 

8.88 

I  (amp) 

4.07 

0.90 

0.98 

3.95 

6 1  (amp) 

0.025 

0.025 

0.025 

0.025 

V  (volts) 

236.26 

52.29 

55.28 

222.05 

SV  (volts) 

0.05 

0.05 

0.05 

0.05 

Qh  (W) 

961.58 

47.06 

54.17 

877.1 

6Qh  (W) 

12.42 

2.83 

3.00 

7.43 

6Qh/Qh  (%) 

1.29 

6.00 

5.55 

0.85 

Do  (mm) 

15.88 

15.88 

13.31 

13.31 

6Do  (mm) 

0.1 

0.1 

0.1 

0.1 

L  (mm) 

203.2 

203.2 

203.2 

203.2 

6  L  ( mm ) 

0.1 

0.1 

0.1 

0.1 

Ab  (m2) 

1  .014E-2 

1 .014E-2 

8.50E-3 

8.50E-3 

6Ab/Ab  (%) 

0.63 

0.63 

0.75 

0.75 

STsat  CC) 

0.01 

0.01 

0.01 

0.01 

Two  CO 

20.18 

12.92 

2.77 

11.13 

6Two  CO 

0.19 

0.28 

8.53E-2 

0.36 

AT  CC) 

17.96 

10.70 

0.53 

9.80 

6AT/AT  (%) 

1.06 

2.62 

16.20 

3.67 

Qf  (W) 

8.99 

5.02 

5.69 

2.91 

5Qf 

0.12 

0.03 

0.15 

2.48E-2 

Qh-2Qf  (W) 

943.60 

37.02 

54.06 

871.3 

q  (W/m2) 

93,081 

3664 

6109 

102,542 

6q/q  (%) 

1.46 

7.67 

5.63 

1.14 

h  (W/m2»K) 

5182 

342.4 

11,200 

10,463 

6h/h  (%) 

1.80 

8.10 

17.15 

3.84 
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APPENDIX  F:  SETUP  PROGRAM 

The  following  is  a  listing  of  the  complete  setup  program 
used  in  the  preparation  of  the  apparatus  for  testing.  It 
consists  of  five  sections: 

1 .  Monitor  the  Sump  temperature 

2.  Monitor  the  evaporator  liquid  temperature 

3.  Measurement  and  readout  of  all  thermocouple  channels. 

4.  Measurement  and  readout  of  the  power  supplied  to  the 
tube  cartridge  heater. 

5.  Measurement  and  readout  of  the  power  supplied  to  the 
auxilliary  heaters  if  used. 

The  program  was  written  in  Hewlett-Packard  Basic  5.0  for  both 
the  Hewlett-Packard  9300  series  computer  and  9852A  series  data 
acquisition/control  unit. 
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1  I  PROGRAM:  SETUP  PR06RAM  FOR  HP  385ZA 

Z  I  DATE:   AUGUST  13.199) 

3   I  PROGRAMMER:    LT  DEAN  SU6IYAMA 

13  REAL    T<0: 1Z).Su«<0: lZ).Te*p<0: 1Z>  ,E<  1:3) 

23  ON  KEY  0.15  GOTO  27 

25  PRINTER   IS    1 

27  PRINT 

28  PRINT 

30  PRINT  USING   "4X. ""SELECT  OPTION 

31  PRINT  USING   "6X.""0-HONITOR  SUMP 

32  PRINT  USING  "6X,""1 -MONITOR  LIQUID 

33  PRINT  USING   "6X.""Z-CHECK   THERMOCOUPLES 

34  PRINT   USIN6   "6X.""3-CHECK   MAIN  HEATER 

35  PRINT  USING   "6X,""4-CHECK   AUX  HEATERS 

36  PRINT  USIN6   nBX.""5-EXIT  PR06RAH""" 

37  PRINT  USING   "4X,""N0TE:    KEY  0  -  ESCAPE 

36  BEEP 

40  INPUT   Ido 

41  IF   Ido>5  THEN  Ido-S 

42  IF   Ido-0  THEN  50 

43  IF    I  do- I    THEN   155 

44  IF  Ido-Z  THEN  165 

45  IF  Ido-3  THEN  209 

46  IF  I  do- 4  THEN  209 

47  IF   Ido-S  THEN  252 

48  PRINT 
491 

50  PRINT 

51  PRINT   "SUMP   TEMPERATURE   0E6  C   " 

53  PRINT 

54  OUTPUT  709i"RST" 

55  Tot-0 

56  FOR  J-1    TO  5 

60  OUTPUT  709i"CONFMEAS   TEMPT. 511 .USE   600" 

90  ENTER  709 i A 

100  Tot-Tot+A 

110  NEXT  J 

120  Tavg-Tot/5 

140  PRINT  USING  "4X,MD0.D0"|Tbv0 

141  BEEP 

142  PRINT 

150  WAIT  60 

151  60T0  50 
1521 

155  PRINT 

156  PRINT   "LIQUID  TEMPERATURE   0E6  C" 

158  PRINT 

159  OUTPUT  709i"RST" 

160  FOR  1-1    TO  2 

162  Sun<I>-0 

163  T(I>-0 
165  NEXT   I 

1B6  FOR  J-1    TO  5 

167  OUTPUT  709«"CONFMEAS  TEMPT, 508-509 .USE   600" 

168  FOR  1-1    TO  2 

169  ENTER  709iT(I) 

170  Sun(I)-Su«<I)+T(I) 

171  NEXT   I 

172  NEXT  J 

173  FOR  1-1   TO  2 
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174  Te«p<I>-Suw<I>/5 

175  NEXT   I 

176  Tavg-( To«p<  1 )+Te«p(2) )/Z 

188  PRINT   USING   "4X.HD0.D0" jTbvq 
181  BEEP 

1BZ  WAIT  &t 

183  GOTO    155 
1841 

185  PRINT 

186  PRINT   "CHANNEL  TEMPERATURE   DEG  C" 

187  OUTPUT  709t"RST" 
168  FOR   1-1    TO    12 

189  Suw<I>-0 

190  T(I)-0 
19Z  NEXT  I 

193  FOR  J-1    TO  5 

194  OUTPUT  709i"CONFHEAS  TEMPT, 500-51 1 .USE   600" 

195  FOR  1-1    TO    1Z 

196  ENTER  709|T<I> 

197  SunU  >-Su«(I>  +  T<  I  ) 

198  NEXT  I 
ZM  NEXT  J 

201  FOR  1-1  TO  1Z 

202  Tenp<I)-Sun<I)/S 

Z63  PRINT  TAB<3)|IiTAB(15)|Tmp(I> 

Z04  NEXT   I 

205  BEEP 

zee  writ  5 

Ze7  GOTO    185 

Z08I 

Z09  PRINT 

Zie  OUTPUT  7e9i"RST" 

211  FOR  1-1    TO   3 

Z1Z  Su«<I>-0 

Z15  NEXT   I 

Z16  FOR  J-1    TO  5 

Z18  OUTPUT  709."CONFMEAS   DCV.51Z-514.USE   600" 

Z19  FOR  1-1    TO   3 

ZZ1  ENTER  709|E(I) 

ZZZ  Sun(I)-Sun(I>+E<I> 

ZZ3  NEXT  I 

ZZ5  NEXT  J 

ZZ6  Anp-0 

ZZ8  FOR  1-1    TO   3 

ZZ9  IF   1-1    THEN  Volt-Sun< I >/5 

Z30  IF  I-Z   AND  Ido-3  THEN 

Z31  PRINT   "HAKE   SURE   VOLTAGE   BOX   IS  SET  TO  MAIN  HEATERS' 

Z3Z  A*p-Su«U>/5 

233  END  IF 

Z34  IF   1-3  AND  Ido-4  THEN 

Z35  PRINT  "MAKE  SURE  VOLTAGE  BOX  IS  SET  TO  AUX  HEATERS" 

Z36  Anp-Sun<I>/5 

237  END  IF 

Z38  NEXT  I 

Z  39  Anp-ABS<  Anp* 1 . 9 1 82 ) 

240  Volt-ABS(VoI-t«25) 

241  Power-VoIt*Anp 

242  Realstence-Volt/Anp 

243  IF  Anp-0  THEN  Reeietence-0 
245  PRINT 
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246 

BEEP 

Z47 

PRINT   "           VOLTAGE(V)        CURRENT(fl) 

248 

PRINT 

249 

PRINT   USING   '•1X.5<MDOOOOO.DOO,3X>" 

250 

UAIT  5 

251 

GOTO  209 

252 

BEEP 

253 

PRINT 

254 

PRINT 

255 

END 

R£SISTENCE(oh«e)        POUER(U)1 
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APPENDIX  G:  DATA  REDUCTION  PROGRAM 

The  following  is  a  listing  of  the  complete  data-reduction 
program,  DRP8,  written  in  Hewlett-Packard  Basic  5.0  for  both 
the  Hewlett-Packard  9300  series  computer  and  9852A  series  data 
acquisition/control  unit. 
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1000!  FILE  NAME:  DRP8 

1004!  DATE:       October  19,  1984 

1006!  REVISED:    August  14.1991  by  LT  Dean  Sugiyana 

1012! 

1016  COM  /I dp/  Idp 

102®  PRINTER  IS  1 

1024  CALL  Select 

1028  INPUT  "UANT  TO  SELECT  ANOTHER  OPTION  < I -Y ,0-N>?" .Iael 

1032  IF  I*el-1  THEN  GOTO  1024 

1036  BEEP 

1040  BEEP 

1044  PRINTER  IS  1 

1048  PRINT  "DATA  COLLECTION/REPROCESSING  COMPLETED" 

1052  END 

1056  SUB  Main 

106®  COM  /Idp/  Idp 

1064  COM  /Cc/  C(7),Ical 

1068  COM  /Uil/   D2,Di,Do.L.Lu,Kcu 

1072  DIM   Enf  (  12) ,T(  12) ,D1a(  13) ,D2a(  13) .Olo< 13).Doa(  1 3) ,La< 1 3) ,Lua< 13),Kcua< 13). 

EtK20>,EtZ<Z0),Tn«<4>[  151  ,Sun(Z0).Tt<  12  )  ,Ta<  12  )  .Vt<  2  >  ,Va(  2  >  ,Eliq<  1Z).E<  12) 

1076  DATA  0.1 008609 1.25727. 94369, -7B7345. 8295.78025595.81 

1080  DATA   -9247486589, 6. 97668E+11 ,-Z . 66192E-H  3 , 3.9407BE+ 1 4 

1084  READ  C<«) 

10881  DATA  "Smooth", "High  Flux" ,"Ther*oexel-E" ."Ther«oexel-HE" 

1092  DATA  Smooth, High  Flux ,Turbo-B, High  Flux  Mod.Turbo-B  Mod 

1096  READ  Tn«(*> 

1100  PRINTER  IS  701 

1 104  BEEP 

1108  IF    Idp-4   THEN  PRINTER   IS    1 

1112  IF    Idp-4  THEN  60T0  2032 

1116  OUTPUT  709«"RST" 

1120  OUTPUT  709«"SET   TIMEDATE"|TIMEDATE 

1124  OUTPUT  709 i "TIHEDATE" 

1128  ENTER  709 ;Dt 

1132  PRINT 

1136  PRINT  "  Month,  date  and  tine  I" |DATE«( Dt ) ,TIME«(Dt ) 

1140  PRINT 

1144  PRINT  USING  " 10X .""NOTE:  Progran  name    I  DRP8 

1148  BEEP 

1152  INPUT  "ENTER  DISK  NUMBER". On 

1156  PRINT  USING  " 1BX,""Di»k  nunber   -  "".ZZ"iDn 

1160  BEEP 

1164  INPUT   "ENTER  INPUT  MODE   (0-3852A. 1-FILE )" .In 

1168  BEEP 

1172  INPUT   "SELECT  HEATING  MODE   <0-ELECi    1-UATER)" ,Ih* 

1176  BEEP 

1180  INPUT  "ENTER  THERMOCOUPLE  TYPE  (0-NEU. 1 -OLD)" .leal 

11 84  IF  I«-«  THEN 

1188  BEEP 

1 192  INPUT  "GIVE  A  NAME  FOR  THE  RAU  DATA  FILE" ,DZ_f lie* 

1196  PRINT  USING  "16X.""New  file  nane 14A" iD2_f ile« 

1Z00  Size1-Z0 

1Z04  CREATE   BDAT  DZ_f ilet.Sizel 

1Z08  ASSI6N  CFlleZ  TO  DZ.flle* 

1Z1ZI 

1216!  DUMMY  FILE  UNTIL  Nrun  KNOUN 

1ZZ0  D1_file»-"OUMMY" 

1Z24  CREATE  BDAT  D1_f ilet.Sizel 

1ZZ8  ASSI6N  CFilel  TO  OI_f lie* 
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1232  OUTPUT  «File1iDt 

123S  IF  Ih«-0  THEN 

1240  BEEP 

1244  INPUT  "ENTER  NUMBER  OF  DEFECTIVE  TCS  (©-DEFAULT )" , Idtc 

1248  IF  Idtc-0  THEN 

1252  Ldtcl-0 

125B  Ldtc2-0 

1260  PRINT  USING  "16X,""No  defective  TCo  exist 

12B4  END  IF 

1268  IF  Idtc-1  THEN 

1272  BEEP 

1276  INPUT  "ENTER  DEFECTIVE  TC  LOCATION" .Ldtcl 

1280  PRINT  USING  "16X.""TC  Is  defective  at  location  "".DD'sLdtcl 

1284  Ldtc2-0 

1288  END  IF 

1292  IF  Idtc-2  THEN 

1296  BEEP 

1300  INPUT  "ENTER  DEFECTIVE  TC  LOCATIONS" .Ldtcl .Ldtc2 

1304  PRINT  USING  •,16X.""TC  are  defective  at  location*  "" .DD.4X.DD" iLdtcl ,Ldtc2 

1308  END  IF 

1312  IF    Idtc>2    THEN 

1316  BEEP 

1320  PRINTER   IS    I 

1324  BEEP 

1328  PRINT   "INVALID  ENTRY" 

1332  PRINTER   IS  701 

1336  60T0    1240 

1340  END  IF 

1344  ENO  IF 

1348  OUTPUT  •File1»Ldtc1.Ldtc2 

13521  In-1    option 

1 356  ELSE 

1360  BEEP 

1364  INPUT   "6IVE   THE   NAME   OF   THE   EXISTING  DATA  FILE" .D2_f ileS 

1368  PRINT  USING   " 16X,M"01d  file   nane:    "" . UA" iD2_f lie* 

1372  ASSIGN  «Fll«2   TO   D2_file» 

1376  ENTER  •File2iNrun 

1380  ENTER  tFile2iDt 

1384  PRINT  USING  "16X.""Thie  data  eet  taken  on  :  "" . 11A" iDATEtC Dt ) 

1388  ENTER  0File2iLdtc1 ,Ldtc2 

1392  IF  Ldtcl >0  OR  Ldtc2>0  THEN 

1396  PRINT  USING  " 16X,"" Thermocouples  uere  defective  at  location*: "" ,2< 3D. 4X>" j 

Ldtcl ,Ldtc2 

1400  END  IF 

1404  ENTER  «File2iItt 

1408  ENO  IF 

1409  Idtc-0 

1411  IF  Ldtc1>0  THEN  Idtc-Idtc+1 

1412  IF  Ldtc2>0  THEN  Idtc-Idtc+1 
14131  IF  I«-0  AND  Ih«-1  THEN  1595 
1416  BEEP 

1420  INPUT  "UANT  TO  CREATE  A  PLOT  FILE7  (0-N. 1-Y>" .Iplot 

1424  IF  Iplot-1  THEN 

1428  BEEP 

1432  INPUT   "GIVE   NAME   FOR  PLOT   FILE".P   file* 

1436  CREATE   BOAT  P_file*,4 

1440  ASSI6N  BPlot    TO  P_flle« 

1444  END  IF 

1448  IF    Ih«-1    THEN 

1452  BEEP 
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USB 

146© 

1464 

1468 

1472 

1476 

1480 

1484 

1488 

I4S2 

1496 

1500 

1504 

1508 

1  51  Z 

1516 

1520 

1524 

1528 

1532 

1536 

1540 

1544 

1548 

1552 

1556 

1560 

1564 

1568 

1572 

1576 

1580 

1584 

1585 

1586 

1588 

1592 

1596 

1600 

1604 

1608 

1612 

1616 

1620 

1624 

1628 

1632 

1636 

1640 

1644 

1648 

1652 

1656 

16601 

16641 

1668! 

16721 

1676 

1680 

1684 


INPUT  "UANT  TO  CREATE  Uo  FILE?  (0-N, 1-Y>" . Iuf 

IF  Iuf-t  THEN 

BEEP 

INPUT  "ENTER  Uo  FILE  NAME".UfileS 

CREATE  BOAT  Ufile$.4 

ASSI6N  BUf lie  TO  UfileB 

END  IF 

BEEP 

INPUT  "UANT  TO  CREATE  Re  FILE?  (0-N, 1-Y>" .Ire 

IF  Ire-1  THEN 

BEEP 

INPUT  "ENTER  Re  FILE  NAME" ,Ref list 

CREATE  BOAT  Refilet.10 

ASSI6N  CRefile  TO  RefileB 

END  IF 

END  IF 

PRINTER  IS  1 

IF  Im-0  THEN 

BEEP 

PRINT  USIN6  "'4X, ""Select  tube  number 

IF  Ihm-0  THEN 

Smooth  4  inch  Ref 

Smooth  4  inch  Cu  <Prees/Sl ido )' 

Soft  Solder  4  inch  Cu""" 

Soft  Solder  4  Inch  HI6H  FLUX 

Uleland  Hard  8  inch 

HI6H  FLUX  8  inch 

6EUA-K  40  Fin*/ in 

6EUA-K  26  Fine/ in 

6EUA-T  19  Fine/in""" 

6EUA-T  OR  6EUA-TY  26  Fine/In 

10  THERHOEXCEL-E 

11  THERHOEXCEL-HE 


EX. 

••-« 

ex. 

..., 

EX. 

""2 

EX. 

••"3 

EX. 

""4 

EX. 

••-5 

EX. 

M-g 

EX. 

""7 

EX. 

--8 

6X. 

•"•9 

6X. 

""\( 

EX. 

""1 

•6X.""12   TURBO-B 

'6X.""13  6EUA-K    19  Fine/in" 


PRINT   USIN6 

PRINT   USING 

PRINT   USIN6 

PRINT   USING 

PRINT   USING 

PRINT  USING 

PRINT   USING 

PRINT  USING 

PRINT  USING 

PRINT  USIN6 

PRINT   USING 

PRINT  USING 

PRIN*  USING 

PRINT  USING 

ELSE 

PRINT   USING 

PRINT  USING 

PRINT  USING 

PRINT   USING 

PRINT   USING 

END  IF 

INPUT   Itt 

OUTPUT  0File1iItt 

END  IF 

PRINTER  IS  701 

IF    Itt<10  THEN  PRINT   USING 

IF   Itt>9   THEN  PRINT  USING 

IF    Ihm-1    THEN  PRINT   USING 

BEEP 

INPUT  "ENTER  OUTPUT  VERSION  <0-LON6. 1 -SHORT. 2-NONE)" .Iov 

BEEP 

INPUT  "SELECT  (0-LIQ. 1-VAP.2-(LIQ*VftP)/2 )" .Ilqv 

DIMENSIONS  FO  TESTED  TUBES 

ELECTRIC  HEATED  HOOE 

D1 -Diameter  at  thermocouple  poeltlone 

DATA  .01 11 125.. 01 11 125.. 01 11 125,. 0129540,. 012446.. 0129540.. 0100965 

DATA  .0100965.. 01 157.. 01 157.. 01 157.. 01 157.. 01 157.. 0100965 

READ  OlotO 


EX. 

-»© 

Smooth  tube"" 

EX, 

...., 

High  Flux 

EX. 

""2 

Turbo-B""" 

EX. 

■•••3 

High  Flux  Hod' 

EX, 

""4 

Turbo-B  Hod""' 

"1  EX. "•'Tube  Number: 
1  EX.  ""'Tube   Number: 
16X.""Tube  Type: 


"".D"iltt 
". DO" i Itt 
",1SA"iTn*<Itt> 
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1668  Dl-Dlo(Itt) 

I6SZI 

16961  Di-Diane-tor  of  'test  section  to  the  base  of  fins 

170©  DATA  . 015875.. ©15875,. ©15875, .©15824 . .©15875. . ©15824, .©1270 

17©4  DATA  .©127, .©138. .©138, .©138, .©138, .©138, .©127 

17©8  READ  D2o(0 

17121 

1716!  Dl-Inslde  diameter  of  unenhanced  ends 

172©  DATA  .©127, .©127. .©127, .©132. .©127. .©132,. ©111125. .©111125 

1721  DATA  .©118..©118..«118,.©118,.©1I8..©1111Z5 

1724  READ  Dia(«> 

17281 

17321  Do-Out  side  diameter  of  unenhanced  ends 

1736  DATA  .©15875, .©15875, .©15875. .©15824. .©15875. .©15824, .©127©, .©127© 

1737  DATA  .©1 331 . .©1 331 . .©1 331 . .©1 331 . .©158. .©127 
174©  READ  Doa<*> 

17441 

17481  L -Length  of  enhanced  surface 

1 75Z  DATA  .1©16..1©16..1©16..1©16.. 2©32 . . Z©32 . . 2©32 , .2©32 . . 2©32 . . 2©32 . . Z©32 . . 2© 

32,.Z©32..2©32 

1756  READ  La(  • ) 

17601 

17641  Lu'Length  of  unenhanced  surface  at  the  ends 

1768  DATA  .©Z54. .©Z54. .©254 , .©254. .©762 . .©762 . .©76Z . .©76Z . .©76Z . .0762 .,0762  .  .©7 

62,. 0762.. ©762 

177Z  READ  Lua<»> 

1776! 

178© I  Kcu- Thermal  Conductivity  of  tube 

1784  DATA  398.344.344.45,344.45.344.344.398.398.398.398.398,344 

1788  READ  Kcua(«) 

179Z  IF  Ibn-1  THEN 

1796! 

16©©!  Data  statements   for   water   heating  node 

18041 

1808  DATA  ©.©15875. ©.©15875, ©.©169. ©.©138. ©.©169. ©.©.©.©.© 

1812  READ  DZafO 

1816  DATA  ©. ©127, ©. ©1Z7, ©. ©145, ©.©1Z7,©. ©145, ©,©.©,©.© 

182©  READ  Dia<»> 

1824  DATA  ©.©15875. ©.©15875. ©.©169 ,©.©15875, ©.©169 ,©,©,©.©.© 

18Z8  READ  DoaO) 

183Z  DATA  0.3048.0. 3048. 0.3048. 0.3048. 0.3048. ©.©.©.©.© 

1836  READ  La(  • > 

1840  DATA  0.0Z54.©.©Z54.©.0Z54.0.©Z54,©.©Z54, ©,0,0.©.© 

1844  READ  Lua< • ) 

1848  DATA  398,45,398.45.398.©.©.©.©.© 

185Z  READ  Kcua(t) 

1856  END  IF 

186©  DZ-OZa(Itt) 

1864  Di-Dia(Itt) 

1868  Do-Ooa(Itt) 

187Z  L-La(Itt) 

1B76  Lu-Lua(Itt) 

188©  Kcu-Kcua<Itt> 

1684  Kn-.B 

1888  Fr-.3 

189Z  IF  Itt-©  THEN  Cf-1.7©E+9 

1896  IF  Itt>©  THEN  Cf-3.7©37E+1© 

19©©  A-PI»(Do*Z-Di*Z)/4 

1904  P-PI»Do 

1908  IF  Ihm-1  THEN 
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1312  BEEP 

1916  INPUT  "TUBE  INITIATION  MODE.  ( 1-HOT  WATER. Z-STEAH.3-C0LD  WATER)", Itin 

19Z0  IF  Itin-1  THEN  PRINT  USING  " 1BX,""Tube  Initiate:  Hot  Water 

1924  IF    Itiw-2    THEN  PRINT   USIN6   " 1BX.""Tube   Initiate:    Steam 

1928  IF  Itin-3  THEN  PRINT  USIN6  "1BX.""Tube  Initiate:  Cold  Water 

193Z  INPUT  "TEHP/VEL  MODE:  <0-T-CONST,V-DECi 1-T-DEC .V-CONSTj  Z-T-INC.V-CONST)1 
Itv 

193B  IF  Itv-0  THEN  PRINT  USING  "16X,""Temp/Vel  Mode:  Constant/Decreasing 

1940  IF  Itv-1  THEN  PRINT  USING  " 1BX . ""Tenp/Vel  Mode:  Decreasing/Constant 

1944  IF  Itv-2  THEN  PRINT  USING  " 1BX  .""Tenp/Vel  Mode:  Increasing/Constant 

1948  INPUT  "WANT  TO  RUN  WILSON  PLOT?  < 1-Y ,0-N)" , Iuil 

195Z  IF  Ih«-1  AND  Iuil-0  THEN 

1956  IF  Itt-0  THEN  Ci-.03Z 

196©  IF  Itt-1  OR  Itt-3  THEN  C1-.059 

1964  IF  Itt-Z  OR  Itt-4  THEN  Ci-.06Z 

1968  BEEP 

197Z  INPUT  "ENTER  CI  (DEF:  WH-.032  .HF-.059  .TB-.062  )"  ,Ci 

1976  PRINT  USING  " 1BX,""Sieder-Tate  """ 

1980  PRINT  USIN6  "I6X,""   Constant   -  "" .Z.4D" iCi 

1984  END  IF 

1988  END  IF 

199Z  IF  Ih*-1  AND  Im-1  AND  Iuil-1  THEN 

1996  IF  Itt-0  THEN  Ci-.032 

Z000  IF  Itt-1  OR  Itt-3  THEN  Ci-.059 

Z004  IF  Itt-Z  OR  Itt-4  THEN  Ci-.06Z 

2W8  ASSI6N  CFileZ  TO  • 

201 2  CALL  Wllson(Cf .Ci) 

2016  ASSI6N  SFileZ   TO  OZ.filev 

Z0Z0  ENTER  CFileZtNrun.Doldt.Ldtcl.LdtcZ.Itt 

2024  END  IF 

20ZB  Neub-0 

Z03Z  IF    Idp-4   THEN   Ihw-1 

Z036  IF    Ih««1    THEN  Nsub-8 

Z040  Ntc-6 

Z041  IF  Ih«-0  THEN  Ntc-IZ 

2044  J-1 

2046  sx-e 

205Z  Sy-0 

Z056  Sxe-0 

2060  S*y-0 

2064  Repeat:  I 

2068  IF  In-0  THEN 

2072  Dtld-Z.ZZ 

2076  Ido-Z 

Z080  ON  KEY  0,15  RECOVER  Z064 

Z084  PRINTER  IS  1 

Z0B8  PRINT  USING  "4X. ""SELECT  OPTION 

Z09Z  PRINT  USIN6   "BX.""0-TAKE   DATA""" 

Z096  IF  Ih*-0  THEN  PRINT  USING  "BX.""1-SET  HEAT  FLUX""" 

Z100  IF  Ihn-1  THEN  PRINT  USING  "BX.""1-SET  WATER  FLOW  RATE 

Z104  PRINT  USING  "BX.""Z-SET  Teat""" 

Z108  PRINT  USING  "4X.""N0TE:  KEY  0  -  ESCAPE 

2112  BEEP 

2116  INPUT  Ido 

2120  IF  Ido>2  THEN  Ido-Z 

Z1Z4  IF  Ido-0  THEN  Z724 

Z1Z8I 

Z13Z!  LOOP  TO  SET  HEAT  FLUX  OR  FLOWMETER  SETTING 

Z136  IF  Ido-1  THEN 

Z140  IF  Ihn-0  THEN 
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21 48  BEEP 

Z15Z  INPUT  "ENTER  DESIRED  Odp",Dqdp 

2156  PRINT  USING  "4X. ""DESIRED  Qdp  ACTUAL  Odp 

216©  Err-1000 

2164  FOR  1-1  TO  2 

2165  Sun(I)-0 

zi66  vtm-e 

2166  NEXT  I 

2169  OUTPUT  709;"RST" 

Z176  FOR  Ji-1  TO  5 

2180  OUTPUT  709j"CONFMEAS  DCV, 512-51 3, USE  600" 

2181  FOR  1-1  TO  2 

21 82  ENTER  709>Vt<I) 

2184  Su«<I)-Su«<I>+Vt(I) 

2185  NEXT  I 

2188  NEXT  Ji 

2189  FOR  1-1  TO  2 

2192  IF  1-1  THEN  Volt-Su«< I >/5 

2196  IF  1-2  THEN  A«p-Su«< I  )/5 

2200  NEXT  I 

Z201  A«p«ABS< A«p« 1 . 9 1 82 ) 

2Z0Z  Vol-t-ABS(VoH»Z5) 

ZZ04  Aqdp-Volt»A«p/<PI»DZ«L) 

ZZ08  IF  ABS<Aqdp-Dqdp)>Err  THEN 

ZZ1Z  IF    Aqdp>Dqdp   THEN 

ZZ16  BEEP   4000,. Z 

ZZZ0  BEEP  4000.. Z  -"* 

ZZZ4  BEEP  4000 ..Z 

ZZZ8  ELSE 

ZZ3Z  BEEP  Z50..Z 

ZZ36  BEEP   Z50..Z 

ZZ40  BEEP   Z50..Z 

ZZ44  END  IF 

ZZ48  PRINT   USING   "4X.HZ.3DE.ZX.nZ. 3DE" lOqdp.Aqdp 

ZZ5Z  WAIT  Z 

ZZ56  60TO  Z164 

ZZ60  ELSE 

ZZ64  BEEP 

ZZ68  PRINT  USING   "4X.HZ.3DE.ZX.HZ. 3DE" lOqdp.Aqdp 

ZZ7Z  Err-500 

ZZ76  UAIT  Z 

ZZ80  60T0  Z164 

ZZ84  END  IF 

ZZ88  ELSE 

ZZ9Z  BEEP 

ZZ96  INPUT  "ENTER  FLOWMETER  SETTING". F«s 

Z300  60T0  Z088 

Z304  END  IF 

Z308  END  IF 

2312! 

Z316I  LOOP   TO  SET   Teat 

Z3Z0  IF    Ido-Z   THEN 

2324  IF   Ikdt-1    THEN  2344 

Z3Z8  BEEP 

Z33Z  INPUT   "ENTER  DESIRED  Teaf.Dtld 

Z336!  PRINT   USING   "4X.""    OTsat      AT»at        Rate  Tv  Rate"' 

Z340  Ikdt-1 

Z344  Oldl-0 

Z348  OldZ-0 

Z35Z  Nn-l 
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2356  Nrs-Nn  HOD  15 

2360  Nn-Nn+1 

2364  IF  Nra-1  THEN 

2365  IF    Ih«-0  THEN  PRINT  USING  "4X,""        Tsat  Tldl  Tld2  Tv  Taun 
P 

2368  IF    IHn-1    THEN  PRINT   USIN6   "4X.""      Teat        Tldl  Tld2  Tv        Taunp      Tinle 

t         Tpile      Tout 

2372  END   IF 

2373  FOR   1-1    TO   6 

2374  Su«<I)-0 

2375  NEXT   I 

2376  OUTPUT  709rRST" 

2377  FOR  Jl-1    TO  20 

2  379  IF    Ihn-0  THEN  OUTPUT  709rCONFHEAS   TEMPT ,508-51 1 .USE   600" 

238®  IF    Ihn-1    THEN  OUTPUT  709t "CONFMEAS   TEMPT, 300-305, USE   600" 

2384  FOR  1-1    TO  6 

2388  IF   Ih*-0  AND   I>4  THEN   2401 

2392  ENTER  709iEllq(I> 

239G  Sund  >-Sim<  I  >+Eliq<  I  ) 

2400  NEXT   I 

2401  NEXT   Jl 

2402  FOR   1-1    TO  6 

2403  IF    Ihn-0  AND  I>4   THEN   2496 
2416  Eliq(I)-Sun(I)/20 

2424  IF    1-1    THEN  Tld1-Eliq<I> 

2428  IF    1-2    THEN   TldZ-Eliq(I) 

2432  IF   1-3  THEN  Tv-Eliq<I) 

Z436  IF  1-4  THEN  Teunp-El lq( I ) 

2440  IF    1-5   THEN  Tinlet-Eliq< I > 

2444  IF    1-6  THEN   Tout-Eliq< I ) 

2448  NEXT   I 

2452  IF    Ih*-1    THEN 

2456  Sun 1-0 

2460  OUTPUT  709j"RST" 

2468  FOR  Kk-1    TO  20 

2469  OUTPUT  709 i "CONFMEAS  TEMPT. 320. USE   600" 
2472  ENTER  709|E1 

2476  Sunl-Sunl+EI 

2480  NEXT  Kk 

2484  Enf(7)-ABS(Sun1/20) 

2488  Tpile-Epif(7)/3.96E-4 

2492  END   IF 

2496  Atld-(Tldl+Tld2>«.5 

Z500  IF   ABS(Atld-Dtld>>.2   THEN 

2504  IF   Atld>Dtld  THEN 

2508  BEEP  4000.. 2 

2512  BEEP  4000.. 2 

2516  BEEP  4000,. 2 

2520  ELSE 

2524  BEEP   250.. 2 

2528  BEEP   250.. 2 

2532  BEEP   250.. 2 

2536  END  IF 

2540  Errl-Atld-Oldl 

2544  Oldl-A-tld 

2548  Err2-Tv-01d2 

2552  01d2-Tv 

2553  IF   Tldl > 100.    THEN  2572 

2556  IF    Ihw-0  THEN  PRINT  USING  "4X.5(HDD0.DD.2X)"»Dtld.Tld1 .Tld2 .Tv.Tbupip 

2560  IF   Ihn-1    AND  Idp-0  THEN  PRINT  USING   "4X,7<  HD0.DD.2X),,iDild.Tld1  ,Tld2  .Tv.Tb 
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unp.T inlet , Tpile 

ZS64  IF    In*-)    AND   Idp-4   THEN  PRINT   USING    "4X ,5( MDO.  DO.ZX ) , 3( M3D. DD.2X >" tDtld. Tl 

dl  ,TldZ  ,Tv,T«upip,T  inlet  .Tpile, Tout 

Z568  WAIT    Z 

Z572  60T0  Z35B 

Z57B  ELSE 

Z5B0  IF   ABS< Atld-Dtld>>.1    THEN 

Z584  IF   Atld>Dtld   THEN 

Z5B8  BEEP   3000.. Z 

Z59Z  BEEP    300®.. Z 

Z596  ELSE 

Z600  BEEP   800,. 2 

ZB04  BEEP   800.. Z 

ZE08  END  IF 

Z61Z  Errl-Atld-Oldl 

ZE1E  Oldl-Atld 

ZBZ0  ErrZ-Tv-OldZ 

ZBZ4  OldZ-Tv 

Z6Z8  IF    Ih»-0  THEN  PRINT  USIN6  "4X,5(  HDDO. DO.ZX)"  iDtld.Tldl  .TldZ  ,Tv  .Teunp 

ZB3Z  IF    Ih«-1    THEN  PRINT   USING   "4X,5< MOO. DO.ZX) ,3<M30. DO. 1X>" jDtld.Tldl .TldZ .Tv 

, Tsunp,T inlet .Tpile, Tout 

ZB36  UflIT   Z 

ZB40  GOTO  Z356 

Z644  ELSE 

Z648  BEEP 

ZB5Z  Errl-Atld-Oldl 

2656  Old  I -At  Id 

ZBB0  ErrZ-Tv-OldZ 

Z664  OldZ-Tv 

ZBB8  IF   Ihn-0  THEN  PRINT   USING   "4X,5<  MDOO. DO.ZX)" iDtld.Tldl , TldZ ,Tv .Tsuwp 

ZB7Z  IF    Ih«-1    THEN  PRINT   USING   "4X.8< MDO. DO.ZX)" iDtld.Tldl .TldZ .Tv.Taunp.Tinlet 

.Tpile .Tout 

ZB76  UAIT  2 

ZE80  60T0  2356 

Z684  END  IF 

Z688  END  IF 

ZBSZ  END  IF 

2696!  ERROR  TRAP  FOR  Ido  OUT  OF  BOUNDS 

2700  IF  Ido>2  THEN 

2704  BEEP 

2708  60T0  2088 

2712  END  IF 

27161 

2720!  TAKE  DATA  IF  I«-0  LOOP 

27Z4  IF  Ikol-1  THEN  Z737 

2728  BEEP 

2732  INPUT  "ENTER  BULK  OIL  X",Bop 

2736  Ikol-1 

Z737  FOR  1-1  TO  Ntc 

Z738  E<I)-0 

2739  Sun< I )-0 

2740  NEXT  I 

2741  OUTPUT  709t"RST" 

2742  FOR  Ji-1  TO  20 

2744  IF  Ihw-0  THEN  OUTPUT  709 t"CONFHEAS  TEMPT .500-511 .USE  600" 

Z745  IF  Ih«-1  THEN  OUTPUT  709i "CONFHEAS  TEMPT . 300- 305. USE  600" 

2748  IF  Ih«-0  THEN  Ntc-12 

2752  FOR  1-1  TO  Ntc 

2756  ENTER  709«E(I> 

2760  Su«(I)-Sun(I)4E<I) 
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2776  IF    I-<9-Nsub>    THEN  Et 1 < Ji- 1  )-E( I > 

2777  IF    I-<10-Nsub)    THEN  EtZ<  Ji-  1  )-E(  I ) 

2778  NEXT    I 
Z780  NEXT  Ji 

2784  Kdl-0 

2785  FOR  1-1    TO  Ntc 

2788  IF    I-O-Naub)    OR   I-<  10-Nsub>    THEN 

2732  Eave-Sum(I>/Z0 

2797  Suw(I>-0 

2798  END  IF 

2800  IF    I-O-Nsub)    THEN 

2801  FOR  Jk-0  TO    19 

2804  IF   ABS<EtKJk)-EaveK.1    THEN 

2808  Su«< I >-Su«< I )+Eil< Jk ) 

2812  ELSE 

2816  Kdl-Kdl+1 

2817  END  IF 

2819  NEXT  Jk 

2820  END  IF 

2822  IF    I-(  10-N»ub>    THEN 

2823  FOR  Jk-e  TO    19 

2825  IF   ABS<Et2<  Jk>-EaveK.I    THEN 

282E  Sun< I)-Su«< I >+EtZ< Jk) 

2828  ELSE 

2829  Kdl-Kdl-H 

2830  END  IF 

2832  NEXT  Jk 

2833  END  IF 

2836  IF    I-<9-Neub>   OR   I-( 10-Nsub)   THEN  PRINT  USIN6   "4X.""Kdl    -   "" .DO" iKdl 

2837  IF   Kdl>10  THEN 

2838  BEEP 
2840  BEEP 

2844  PRINT   USING   "4K.""Too  mien  scattering   In  data  -    repeat   data  set 

2848  60T0  2084 

2852  END  IF 

2860  Enf(I)-Sun< I )/(20-Kdl  ) 

2864  NEXT   I 

2868  IF   Ihn-1    THEN 

2872  Su«( 1 )-0 

2876  OUTPUT  709i"CONFHEftS   TEff T. 320 .USE   600" 

2884  FOR  Kk-1    TO   20 

2888  ENTER  709iE(Kk) 

2892  Sun< 1 )-Sun(  1  >+E(Kk  ) 

2896  NEXT  Kk 

2900  E«f<7)-ftBS<Su«<1 ))/20 

2904  END  IF 

2908  IF  Ihn-0  THEN 

2909  FOR  1-1  TO  2 

2910  Su«(I)-0 

2911  NEXT  I 

29 1Z  OUTPUT  709i"RST" 

2913  FOR  Ji-1  TO  5 

2915  OUTPUT  709t"CONFMEn*S  DCV, 512-513. USE  600" 

2916  FOR  1-1  TO  2 
2920  ENTER  709»E(I) 

2936  Sun(I)-Sun(I)4E(I) 

2937  NEXT  I 

2940  NEXT  Ji 

2941  FOR  1-1  TO  2 

2944  IF  1-1  THEN  Vr-Su«<I>/5 
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2948  IF  1-2  THEN  Ir-Su«( I )• 1 .9182/5 

2952  NEXT  I 

2956  END  IF 

2960  ELSE 

2964  IF  Ihw-0  THEN  ENTER  «File2 iBop.Told.Enf < • ) .Vr , Ir 

2968  IF  Ihn-1  THEN  ENTER  8File2iBop.ToId.Enf  (  • ) .Fns 

2972  END  IF 

29761 

29801  CONVERT  e«f' S  TO  TEMP .VOLT. CURRENT 

2984  Tuo-0 

2988  FOR  1-1  TO  Ntc 

2992  IF  ldtc>0  THEN 

2996  IF  I-Ldtcl  OR  I-Ld"tc2  THEN 

3000  T(I)-0 

3004  60T0  3044 

3008  END  IF 

3012  END  IF 

3016  IF  Itt<4  AND  Ih«-0  THEN 

3020  IF  I>4  AND  I<9  THEN 

3024  T(I>-0 

3028  60T0  3044 

3032  END  IF 

3036  END  IF 

3040  T( I >-Enf i I > 

3044  NEXT  I 

3048  IF  Itt<4  THEN 

3052  FOR  1-1  TO  4 

3056  IF   I-Ldtcl    OR   I-Ldtc2    THEN 

3060  Tuo-Tuo 

3064  ELSE 

3068  Tua-Tua+T< I ) 

3072  END  IF 

3076  NEXT  I 

3080  Tua-Tya/(4-Idtc> 

3084  ELSE 

3088  IF   IhM-1    THEN   3128 

3092  FOR   1-1    TO  8 

3096  IF   I-Ldicl    OR  I-Ldic2   THEN 

3100  Tua-Tuo 

3104  ELSE 

3108  Tua-Tua+TU)  • 

3112  END  IF 

3116  NEXT  I 

3120  Tu-Twa/O-Idtc) 

3124  END  IF 

3128  Tld-T<9-Nsub> 

3132  Tld2-T(10-Naub) 

3136  Tlda-(Tld+TldZ)« .5 

3140  Tv-T( 11-Nsub) 

3158  Teunp-T<12-Nsub> 

3160  IF  Ih«-0  THEN  3176 

3168  Tinlei-T( 13-Nsub) 

3172  Tout-T( 14-Neub) 

3176  IF  Ih«-0  THEN 

3180  Anp-ABS(Ir) 

3184  VoH-flBS(Vr)«25 

3188  Q-Volt«A«p 

3192  END  IF 

3196  IF  Itt-0  AND  Ihn-0  THEN 

3200  Kcu-FNKcu(Tu) 


181 


3204  ELSE 

3208  Kcu-KcuaUtt ) 

3212  END  IF 

3216! 

3220!  FOURIER  CONDUCTION  EQUATION  WITH  CONTACT  RESISTANCE  NEGLECTED 

3224  IF  Ih*-0  THEN  Ty-Ty-Q»LOG< D2/D1 )/< 2«PI»Kcu»L > 

3228  IF  Ilqv-0  THEN  Taat-Tlda 

3232  IF   Ilqv-1    THEN  Tsai-( Tlda+Tv )• . 5 

3236  IF  Ilqv-2  THEN  Teat-Tv 

3244  IF  Ih«-1  THEN 

3248  Tavg-Tinlet 

3252  6rad-37. 9853+. 104388« Tavg 

3256  Tdrop-ABS(Ewf(7)>»1.E*6/< 10»Grad> 

326®  Tavgc-Tinle+.-Tdrop«.5 

3264  IF  ABS<Tavg-TavQc)>.01  THEN 

3268  Tavo-< Tavg+Tavgc )• .5 

3272  60T0  3252 

3276  END  IF 

32801 

32841  COMPUTE  UATER  PROPERTIES 

3288  IF  Ihn-1  THEN 

3292  Kw-FNKy(Tavg) 

3296  Huuo-FNMuu(  Tavg) 

3300  Cpy-FNCpy(Tavg) 

3304  Pr y-FNPry<  Tavg ) 

3308  Rhow-FNRhou(Tavg) 

3312  Tyi-Tavg 

3316! 

3320!  Compute    MOOT 

3324  Mdo-t-3.9657E-3+FM»<3.6t9S5E-3-Fn6»<B.820O6E-6-F«s»<  1 .23688E-7-F«s»4. 31897 

E-10))) 

33281  Ndo+.-Hdot«<  1 .0365-Tinle+.M  1 .96644E-3-Tinle+,«5.252E-6)  >/1  .0037 

3332  Kdt-0 

3336  Q-Mdot«Cpu»Tdrop 

3340  L«td-Tdrop/L06< ( Tinlet-Teat )/( Tinlet-Tdrop-Tsat. ) ) 

3344  Uo-Q/(PI»0o«L«L«td> 

3348  Ru-Do«LOG(Do/Oi )/(Z.«Kcu) 

3352  Ty-T»at*Fr«Lw*d 

3356  Vy-Ndot/<Rhoy«PI»Di*Z/4> 

3360  Reu-Rhou*Vu»Dl/Muuo 

3364  Hi-Cl«Ku/Di»Reu*.8»Pru-(  1/3.  )•( Huua/FNMuy( Tyl )  ) " .  14 

3368  Twic-TBvg-Q/(PI«Do«L«Hi) 

3372  IF  ABS<Tui-Tyic)>.01  THEN 

3376  Tui-<Twi+Tuic)«.5 

3380  60T0   3364 

3384  END  IF 

3388  Tyl-(Tui+Tuic)«.5 

3392  Ho-1/<1/Uo-Do/<0i«Hi)-Ry> 

3396  END  IF 

3397  END  IF 

3399  IF    Ih«-1    THEN 

3400  Thetab-Q7<Ho»PI»Do«L> 

3404  Ty-T»at+Thatab 

3405  ELSE 

3408  Thetab-Ty-Tsat 

3411  END  IF 

3413  IF  Thetab<0  THEN 

3416  BEEP 

3420  INPUT  "TUALL<TSAT  (0-CONTINUE.  1-END)",Iav 

3424  IF  lev-®  THEN  60T0  2068 
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34Z8  IF    Iev-1    THEN   383Z 

3432  END  IF 

34401 

3444!  COMPUTE   VARIOUS   PROPERTIES 

345Z  Tf il*-(Tu+Teat)«.5 

3456  Rho-FNRho(Tfiln) 

3460  Hu-FNMu<Tfilm) 

3464  K-FNK(Tfiln) 

3468  Cp-FNCp(Tfilm) 

347Z  Beta-FNBeta<Tf iln) 

3476  Hfg-FNHfg(Tsat) 

3480  Ni-Mu/Rho 

3484  Alpha-K/(Rho»Cp) 

3488  Pr-Ni/Alpha 

3492  Pea-t-FNPBai(  Tsat  > 

3496! 

35001  COMPUTE  NATURAL-CONVECTIVE  HEAT-TRANSFER  COEFFICIENT 

3504!  FOR  UNENHANCED  ENO(S) 

3508  Hbar-190 

351Z  Fe-(Hbar«P/(Kcu«A))*.5«Lu 

3516  Tanh-FNTanh(Fe) 

35Z0  Theta-The"tab«Tanh/Fe 

3SZ4  Xn-< 9. 81«Bata»Thetab»Do*3«Tanh/<Fa»Ni«Alpha>>*. 166667 

35Z8  Yy-<1+< .5S9/Pr >*< 9/16) >* < 8/Z7) 

353Z  Hbarc-K/Do«(  .6-».387»Xx/Yy)*2 

3536  IF  ABS<<Hbar-Hbarc>/Hbarc)>.001  THEN 

3540  Hbar-(Hbar+Hbarc>*.5 

3544  GOTO  3512 

3S48  END  IF 

355ZI 

3556!  COMPUTE  HEAT  LOSS  RATE  THR0U6H  UNENHANCED  END(S) 

3560  01 -( Hbar«P«Kcu« A ) * . 5*  Thetab«  Tanh 

3564  Qc-Q-Z»Q1 

3568  Ae-PI*DZ«L 

357Z! 

3576!  COMPUTE  ACTUAL  HEAT  FLUX  AND  BOILING  COEFFICIENT 

3580  Qdp-Qc/Ae 

3584  Mube-Qdp/Thatab 

3588  Cef-<Cp»Thetab/Hfg)/<Qdp/<Hu»Hfg>»<  .014/< 9.81»Rho>* .5)* <  1/3. >»Pr*1 .7) 

35SZI 

35961  RECORD  TIME  OF  DATA  TAKING 

3600  IF  In-0  THEN 

3604  OUTPUT  709 i" TIME DATE" 

3608  ENTER  709 » Told 

3612  END  IF 

36161 

36Z0I  OUTPUT  DATA  TO  PRINTER 

36Z4  PRINTER   IS  701 

36Z8  IF    Iov-0  THEN 

363Z  PRINT 

3636  PRINT   USING  " 10X.""Data  Sat   Nunber   -   "" .DOO.ZX.-'Bulk   Oil    X  -   "" .00.0.5X. 1 

4A"iJ.Bop.TIHE»(Told) 

3640  IF   Ihn-0  THEN 

3644  PRINT   USING   ,'10X.""TC   No:  1  Z  3  4  5  6  7 

8""" 

3648  PRINT  USING   " 10X,""Te«p   :"" ,8< 1X.HD0.D0>" |T< 1 > ,T< Z > ,T( 3) ,T<4> ,T< 5> ,T< 6) .T< 
7).T<8) 

365Z  PRINT   USIN6   "10X.""   Tub  Tliqd        TliqdZ      Tvapr     Peai        T»u«p 

3656  PRINT   USING   " 10X,Z< HDO.DD. IX) , IX, MOO. DO. IX. Z< 1X.MD0.D0) ,ZX,HDD.D"iTu,Tld.T 
ldZ .Tv.Psat .Tsunp 
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3660  PRINT  USING  "10X,""  Thetab   Htube      Odp 

3GG4  PRINT  USING  " 10X .HDD. 3D, 1 X.HZ. 3DE . IX ,MZ. 3DE" j Thetab .Htube .Odp 

3GG8  ELSE 

3G7Z  PRINT  USING  "lex,""  Fns    Vu    Tsat   Tinl  Tdrop  Thetab     q         Uo 

Ho 

3G7G  PRINT  USING  " 10X ,4< ZD.OO, 1 X) ,Z . 3D. 1 X ,00. DO, I  X , 3( MZ. 3DE , 1 X >" jFns , Vu .Tsat ,Ti 

nlet , Tdrop. Thetab. Odp, Uo. Ho 

3680  ENO  IF 

3GB4  END  IF 

3688  IF  Iov-1  THEN 

369Z  IF  J-1  THEN 

3696  PRINT 

370®  IF  Ih«-0  THEN 

3704  PRINT  USING  "10X.""  RUN  No  OilX  Tsat    Htube      Qdp       Thetab 

3708  ELSE 

37IZ  PRINT   USING   "10X.""      FMS  OILX     TSAT  HTUBE  QDP  THETAB 

3716  END  IF 

37Z0  END  IF 

37Z4  IF    Ihn-0   THEN 

37Z8  PRINT   USING   " 1ZX.3D.4X.DO,ZX,HOO.D0.3<  1X.HZ.3DE >"i J. Bop. Tsat .Htube. Qdp. The 
tab 

373Z  ELSE 

3736  PRINT  USING  " 1ZX  .3D.4X,D0.ZX,HD0.D0.3< 1X.HZ.3DE>" iF«e. Bop, Tsat .Htube. Qdp, T 
hetab 

3740  ENO  IF 

3744  END  IF 

3748  IF  I«-0  THEN 

375Z  BEEP 

3756  INPUT  "OK  TO  STORE  THIS  DATA  SET  ( 1-Y.O-N)?" ,0k 

37B0  END  IF 

3764  IF  Ok-1  OR  I«-1  THEN  J-J+1 

3768  IF  Ok-1  AND  In-0  THEN 

377Z  IF  Ih*-0  THEN  OUTPUT  •File! iBop. Told. E«f< •) ,Vr ,Ir 

3776  IF  Ihn-1  THEN  OUTPUT  «File1 »Bop.Told.E«f  (  • > ,Fns 

3780  END  IF 

3784  IF  Iuf-1  THEN  OUTPUT  eUfileiVu.Uo 

3788  IF  Ire-1  THEN  OUTPUT  «Ref ileiFns.Reu 

379Z  IF  <I*-1  OR  Ok-1)  AND  Iplot-1  THEN  OUTPUT  ePlot i Qdp. Thetab 

3796  IF  I«-0  THEN 

3800  BEEP 

3804  INPUT  "UILL  THERE  BE  ANOTHER  RUN  ( 1-Y.O-N)?" .Go_on 

3808  Nrun-J 

381 Z  IF  Go_on-0  THEN  383Z 

3816  IF  6o_on<>0  THEN  Repeat 

38ZO  ELSE 

38Z4  IF  J<Nrun+l  THEN  Repeat 

38Z8  END  IF 

3B3Z  IF  I«-0  THEN 

3836  BEEP 

3840  PRINT  USING  M10X.M"NOTE:  "".ZZ.""  data  rune  were  stored  in  file  "" . 10A" iJ- 
I.DZ.filef 

3844  ASSIGN  fFilel  TO  • 

3848  OUTPUT  eFileZiNrun- 1 

385Z  ASSIGN  •Fllel  TO  D1_file» 

3856  ENTER  "Filel iDt .Ldtcl .LdtcZ .Itt 

3860  OUTPUT  BFileZiOt .Ldtcl .LdtcZ ,Itt 

3864  FOR  1-1  TO  Nrun- 1 

3868  IF  Ihw-0  THEN 

387Z  ENTER  CFilel iBop.Told.Enf ( • ) ,Vr ,Ir 

3876  OUTPfT  «FileZiBop.Told,Enf < • > .Vr .Ir 
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3880 
3884 
3888 
3892 
3896 
3900 
3904 
3908 
3912 
39  IB 
3920 
plot 
3924 
3928 
3932 
39  36 
3940 
3944 
3948 
3952 
3956 
3960! 
39641 
3968 
3972  I 
3976 
3980 
3984 
3988 
3992 
39961 
4000 
4004 
4008 
4012 
4016 
40201 
4024 
4028 
4032 
4036 
4040 
4044 
4048 
4052 
4056 
4060 
4064 
4068 
4072 
4076 
4080 
4084 
4088 
4092! 
4096 
4100 
4104 
4108 
4112 


ELSE 

ENTER  «File1iBop.Told.E*f<O.F«s 

OUTPUT   •File2iBop.Told.Eiif  (•  )  .Fms 

END  IF 

NEXT   I 

ASSIGN  «File1  TO  • 

PUR6E  "DUMMY" 

END  IF 

BEEP 

PRINT 

IF  Iplot-1  THEN  PRINT  USING  "  10X.""NOTE ,2Z.""  X-Y  pairs  uere  stored  in 

data  file  "" . 10A"i J- 1 ,P_f ile$ 

ASSIGN  BFileZ  TO  • 

ASSI6N  Wlot  TO  • 

IF  Iuf-1  THEN  ASSIGN  0Ufile  TO  • 

IF  Ire-1  THEN  ASSI6N  0Refile  TO  • 

CALL  Stats 

BEEP 

INPUT  "LIKE  TO  PLOT  DATA  < 1-Y.0-N)?" .Ok 

IF  Ok-1  THEN  CALL  Plot 

SUBEND 

CURVE  FITS  OF  PROPERTY  FUNCTIONS 

DEF  FNKcu<T) 

OFHC   COPPER  250  TO  300  K 

Tk-T+273.1S  IC   TO  K 

K-434-.112«Tk 

RETURN  K 

FNEND 

DEF  FNMu(T) 

170  TO   360  K   CURVE   FIT  OF   VISCOUSITY 

Tk-T+273.15  IC  TO  K 

Mu-EXP( -4.4636-K 101 1 .47/Tk ) )• 1 .0E-3 

RETURN  Mu 

FNEND 

DEF   FNCp(T) 

180  TO  400  K   CURVE   FIT   OF   Cp 

Tk-T+273.15  IC   TO  K 

Cp-.40188+1.65007E-3»Tk+1.S1494E-6«Tk*2-6.67853E-10«TkA3 

Cp-Cp«1000 

RETURN  Cp 

FNEND 

DEF  FNRho<T) 

Tk-T+273.15  IC   TO   K 

X-1-<1.8»Tk/753.9S)      IK    TOR 

Ro-36. 32*61. U6414«X*< 1/3)4 16. 418015* X417.476838«X* 

Ro-Ro/. 062428 

RETURN  Ro 

FNEND 

DEF   FNPr<T> 

Pr-FNCp<T)»FNMu(T)/FNK(T) 

RETURN  Pr 

FNEND 

DEF  FNK(T) 

T<360  K   WITH  T   IN  C 

K-. 071 -. 000261 «T 

RETURN  K 

FNEND 

DEF  FNTanh<X> 

P-EXP(X) 


.5+1.119828«X*2 
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411 B  0-1 /P 

412®  Tanh-(P-0)/<P+0) 

41 24  RETURN   Tanh 

4128  FNENO 

4138! 

4184  DEF    FNBeto(T) 

4188  Rop-FNRho<T+. 1 > 

4192  Ron-FNRho< T- . 1 > 

4136  Beta--2/<Rop+Row>»<Rop-Ro«>/.2 

42®©  RETURN  Beta 

4204  FNEND 

4208  D€F   FNHfg<T) 

4212  Hf g-1. 3741 344E+5-T»< 3. 3094361 E+2+T* 1.2 1651 43) 

4216  RETURN   Hfg 

422®  FNENO 

4224  DEF   FNPaat(Tc) 

42281  ®  TO   B«  dag  F   CURVE   FIT   OF   Paat 

4232  Tf-1.8«Tc+32 

4236  Pa-5.945525+Tf»< . 153S2082+Tf •< 1 .4840963E-3+Tf •9.615«671E-6) ) 

424®  Pq-Pb-14.7 

4244  IF  Pg>0  THEN       I  ♦-PSI6.— in  Hg 

4248  Peat-Pg 

4252  ELSE 

4256  PBat-Pg»29.92/14.7 

426®  END  IF 

4264  RETURN  P»at 

4268  FNEND 

4272  DEF   FNHa«ooth< X ,8op . Iaat ) 

4276  DIM  fl(5) ,B(5).C(5),D(5) 

428®  DATA  . 2®52B .. 25322 .. 31 9®48 , .55322 . .799®9 . 1 .®0258 

4284  DATA  .74515. .72992 , .73189 , .71 225 . .68472 . .64197 

4288  DATA  .41®92 ,. 17726, . 25142 , .548®6, .819 16, 1 .0845 

4292  DATA  .714®3, .729 1 3, .72565, .69669 1 , .665867 ,. 61889 

4296  READ  A( • ) ,B( • ) ,C< • ) ,D< • ) 

43®®  IF  Bop<6  THEN  I-Bop 

43®4  IF  Bop-6  THEN  1-4 

43®8  IF    Bop-1®  THEN   1-5 

4312  IF  Ieat-1  THEN 

4316  H»-EXP<A(I)+B<I >»L0G<X> > 

432®  ELSE 

4324  Hs-EXP<C<I)*D<I)»L06<X>) 

4328  END  IF 

4332  RETURN   Ha 

4336  FNEND 

434®  DEF    FNPoly(X) 

4344  COM  /Cply/   A< 10, 10) ,C< 10) ,B< 5) .Nop .Iprnt ,0po , I  log, Ifn, I Join.Nj oin 

4348  Xl-X 

4352  Poly-B<®> 

4356  FOR   1-1    TO  Nop 

436®  IF    I  log- I    THEN   Xl-LOG(X) 

4364  Poly-Poly+B< I)«X1*I 

4368  NEXT   I 

4372  IF    Ilog-1    THEN  Poly-EXP<Poly ) 

4376  RETURN  Poly 

438®  FNEND 

4384  SUB  Poly 

4388  DIM   R( 10). S( 1®).Sy( 12).Sx( 12).Xx<  100), Yy< 1®®) 

4392  COM   /Cply/    A< 10, 10) ,C< 10) ,B< 5) ,N. Iprnt .Opo.Ilog. If n, IJoln.Njoin 

4396  COM   /Xxyy/    Xp( 25) . Yp( 25) 

44®®  FOR   I-®   TO  4 
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4404  BU1-0 

4408  NEXT  I 

4412  BEEP 

4416  INPUT  "SELECT  (  ©-FILE . 1-KEYB0ARD.2-PROGRAH)" .1* 

44Z0  In-IMI 

44Z4  BEEP 

4428  INPUT  "ENTER  NUMBER  OF  X-Y  PAIRS". Np 

4432  IF  I«-1  THEN 

4436  BEEP 

4440  INPUT  "ENTER  DATA  FILE  NAME" .D_fi let 

4444  BEEP 

4448  INPUT  "LIKE  TO  EXCLUDE  DATA  PAIRS  < l-Y.0-N>7" . lad 

4452  IF  led- I  THEN 

445B  BEEP 

44B0  INPUT  "ENTER  NUMBER  OF  PAIRS  TO  BE  EXCLUDED" .Ipex 

4464  END  IF 

4468  ASSI6N  fiFlle  TO  D_file« 

4472  ELSE 

4476  BEEP 

4480  INPUT  "UANT  TO  CREATE  A  DATA  FILE  (  1-Y.0-N)?" .Yee 

4484  IF  Yee-1  THEN 

4488  BEEP 

4492  INPUT  "6IVE  A  NAME  FOR  DATA  FILE" ,D_f ile« 

4496  CREATE  BOAT  D_file».5 

4500  ASSIGN  BFlle  TO  D_file» 

4504  END  IF 

4508  END  IF 

4512  BEEP 

451B  INPUT  "ENTER  THE  ORDER  OF  POLYNOMIAL" ,N 

4520  FOR  1-0  TO  N«2 

4524  Sy< I >-0 

4528  Sx(I>-0 

4532  NEXT  I 

4536  IF  Ied-l  AND  In-1  THEN 

4540  FOR  1-1  TO  Ipex 

4544  ENTER  «File»X,Y 

4548  NEXT  I 

4552  END  IF 

4556  FOR  1-1  TO  Np 

4560  IF  Ir-1  THEN 

45B4  IF  Opo-2  THEN  ENTER  BFileiX.Y 

4568  IF  Opo<2  THEN  ENTER  CFileiY.X 

4572  IF  Opo-1  THEN  Y-Y/X 

4576  IF  Ilog-1  THEN 

4580  IF  Opo-2  THEN  Xt-X/Y 

4584  X-L06<X) 

4588  IF  Opo-Z  THEN  Y-L06< Xt ) 

4592  IF  Opo<2  THEN  Y-L06(Y) 

4596  END  IF 

4600  END  IF 

4604  IF  In-2  THEN 

4608  BEEP 

4612  INPUT  "ENTER  NEXT  X-Y  PAIR".X.Y 

4616  IF  Ye«-1  THEN  OUTPUT  «FileiX.Y 

4620  END  IF 

4624  IF  I«<3  THEN 

4628  Xx<I>-X 

4632  Yy<I)-Y 

4636  ELSE 

4640  X-Xp(I-l) 
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4644  Y-Yp(I-1 )  ' 

4648  END  IF 

4652  R<0)-Y 

4656  Sy(0)-Sy<0)+Y 

4660  S( 1  >«X 

4664  Sx< 1 >-Sx( 1>  +  X 

4668  FOR  J- I  TO  N 

467Z  R( J)-R< J-1 >«X 

4676  Sy( J)-Sy< J)+R( J) 

4680  NEXT  J 

4684  FOR  J-2  TO  N»2 

4688  S< J)-S< J-1)«X 

4692  Sx< J)-Sx( J)+S< J) 

4696  NEXT  J 

4700  NEXT  I 

4704  IF  Yee-1  AND  I«-2  THEN 

4708  BEEP 

4712  PRINT  USING  "12X.DO,""  X-Y  pairs  were  stored  in  file  "" . 10A";Np.D_f ilet 

4716  END  IF 

4720  Sx<0)-Np 

4724  FOR  1-0  TO  N 

4728  C(I)-Syd) 

4732  FOR  J-0  TO  N 

4736  Ad.J)-Sxd+J> 

4740  NEXT  J 

4744  NEXT  I 

4748  FOR  1-0  TO  N- 1 

4752  CALL  Divided) 

4756  CALL  Subtract (I +1 > 

4760  NEXT  I 

4764  B(N)-C(N)/A(N.N) 

4768  FOR  1-0  TO  N- 1 

4772  B(N-1-I)-C(N-1-I) 

4776  FOR  J-0  TO  I 

4780  B<  N-  1  - 1  >-B<  N- 1  - 1  >  -  A<  N-  1  - 1 .  N-  J  >  »B<  N- J  > 

4784  NEXT  J 

4788  B(N-1-I)-B(N-1-I)/A(N-1-I.N-1-I) 

4792  NEXT  I 

4796  I  PRINTER  IS  701 

4800  'PRINT  B<*> 

4804  I  PRINTER  IS  705 

4808  IF  Iprnt-0  THEN 

4812  PRINT  USING  " 1 2 X." "EXPONENT    COEFFICIENT""" 

4816  FOR  1-0  TO  N 

4820  PRINT  USING  " 15X.D0.5X.M0.7DE" i I ,B(I ) 

4824  NEXT  I 

4828  PRINT  "  " 

4832  PRINT  USING  "12X.""DATA  POINT      X  Y      Y< CALCULATED)  DISCR 

EPANCY 

4836  FOR  1-1  TO  Np 

4840  Yc-B(0) 

4844  FOR  J-1  TO  N 

4848  Yc-Yc+B<  J)»Xxd)*J 

4852  NEXT  J 

4856  D-Yyd)-Yc 

4860  PRINT  USING  "15X. 3D.4X,4(MD.5DE , IX)" i I ,Xx( I ) ,Yyd ) , Yc.D 

4864  NEXT  I 

4868  END  IF 

4872  ASSIGN  CFile  TO  • 

4876  SUBEND 
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4880  SUB  Divide* M) 

4884  COM  /Cply/  A< 10. 10). C< 10>,B<5> .N.Iprnt  ,Opo,  Ilog.Hn.IJoin.Njoiri 

4888  FOR  I-M  TO  N 

4892  Ao-A(I.M) 

4896  FOR  J-M  TO  N 

4900  A(I,J>-A(I,J)/Ao 

4904  NEXT  J 

4908  C(I>-C(I >/Ao 

49)2  NEXT  I 

49  IE  SUBEND 

4920  SUB  Subtract ( K) 

4924  COM  /Cply/  A< 10. 10) ,C< 10) ,B< 5) .N.Iprnt ,Opo, Ilog.If n.IJoin.Njoin 

4928  FOR  I-K  TO  N 

4932  FOR  J-K-1  TO  N 

493B  A(I.J)-A<K-1 .J)-A(I.J) 

4940  NEXT  J 

4944  C(I)-C<K-1)-C<I  ) 

4948  NEXT  I 

4952  SUBEND 

495B  SUB  PI  in 

4960  COM  /Cply/  A<  10, 10) .C<  10) ,B< 5) .N.Iprnt .Opo.Ilog, Hn, Ijoin.Njoin 

49B4  COM  /Xxyy/  Xx< 25) . Yy< Z5> 

4968  PRINTER  IS  705 

4972  BEEP 

4976  INPUT  "UANT  TO  PLOT  Uo  vs  Vy7  < I-Y.0-N)" , Iuo 

4980  IF  Iuo-0  THEN 

4984  BEEP 

4988  INPUT   "SELECT   <0-h/h0X   sane   tube.  1-h(HF  >/h<  siO"  ,Irt 

4992  BEEP 

4996  INPUT  "SELECT  h/h  RATIO  < 1-FILE ,0-COMPUTED)" .Ihrat 

5000  IF  Ihrat-0  THEN 

5004  BEEP 

5008  INPUT  "UHICH  Teat  < 1-6.7,0— 2 .2 )" .Ieat 

5012  END  IF 

5016  Xnln-0 

5020  Xnax-10 

5024  Xetep-2 

5028  IF  Irt-0  THEN 

5032  Ynin-0 

5036  Ynax-1.4 

5040  Yetep-.Z 

5044  ELSE 

5048  Ynin-0 

5052  Y*ax-15 

5056  Yetep-5 

5060  END  IF 

5064  ELSE 

5068  Opo-2 

5072  Y«in-0 

5076  Ynax-12 

5080  Ystep-3 

5084  Xmin-0 

5088  Xnax-4 

5092  Xetep-I 

5096  END  IF 

5100  IF  Ihrat-1  THEN 

5104  Ynin-0 

5108  Ynax-18 

5112  Ystep-3 

5116  Xmn-0 
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5120  Xnax-9 

5124  Xstep-Z 

51Z8  ENO  IF 

5132  BEEP 

5136  PRINT  "IN;SP1iIP  2 300, 2200. 8300, 68O0i " 

5140  PRINT  "SC  0, 100.0, 100iTL  2,0)" 

5144  Sfx-100/(X«ax-X*ln> 

5148  Sfy-IOO/(Y*ax-Y*in) 

5152  PRINT  "PU  0.0  PD" 

5156  FOR  Xa-X«in  TO  Xnai  STEP  Xatep 

5160  X-<Xa-X*in)»Sfx 

5164  PRINT  "PA"iX.",Ot  XT;" 

5168  NEXT  Xa 

5172  PRINT  "PA  IOO.OiPUi" 

5176  PRINT  "PU  PA  0,0  PD" 

5180  FOR  Ya-Y«in  TO  Ynax  STEP  Yetep 

5184  Y-(Ya-Y«in)»Sfy 

5188  PRINT  "PA  0."|Y,"YT" 

5192  NEXT  Ya 

5196  PRINT  "PA  0.100  TL  0  2" 

5200  FOR  Xa-Xnln  TO  Xnax  STEP  Xetep 

5204  X-<Xa-X«in>«Sfx 

5208  PRINT  "PA"iX.",100i  XT" 

5212  NEXT  Xa 

5216  PRINT  "PA  100,100  PU  PA  100.0  PD" 

5220  FOR  Ya-Ymn  TO  Ynax  STEP  Yetep 

5224  Y-< Ya-Y«in)»Sfy 

5228  PRINT  "PD  PA  100," ,Y."YT" 

5232  NEXT  Ya 

5236  PRINT  "PA  100.100  PU" 

5240  PRINT  "PA  0.-2  SR  1.5,2" 

5244  FOR  Xa-XMin  TO  Xnax  STEP  Xstep 

5248  X-<Xa-X«in)«Sfx 

5252  PRINT  "PA"iX,".0i" 

5256  IF  Iuo-0  THEN  PRINT  "CP  -2 .-1 iLB" «Xai"" 

5260  IF  Iuo-1  THEN  PRINT  "CP  -  1 .5.-1 iLB" iXai "" 

5264  NEXT  Xa 

5268  PRINT  "PU  PA  0,0" 

5272  FOR  Ya-Ywin  TO  Ywax  STEP  Yetep 

5276  IF   ABS<YaX1.E-5  THEN   Ya-0 

5280  Y-(Ya-Y«in>»Sfy 

5284  PRINT  "PA  0,"jY,"" 

5288  IF  Iuo-0  THEN  PRINT  "CP  -4,-.25lLB"iY«l"M 

5292  IF  Iuo-1  THEN  PRINT  "CP  -3.- .25iLB" i Yai "" 

5296  NEXT  Ya 

5300  Xlabel»-"Oil    Percent" 

5304  IF    Iuo-0   THEN 

5308  IF    Irt-0  THEN 

5312  Ylabel»-"h/hOX" 

5316  ELSE 

5320  Ylabel»-"h/he«ooth" 

5324  END  IF 

5328  PRINT   "SR    1.5.2tPU  PA  50. -10  CP"i -LEN( Xlabelt)/2i "OiLB" sXlabelSi"" 

5332  PRINT   "PA  -11,50  CP  0." i -LEN(  Ylabeli )/2»5/6» "DI   0. 1 »LB" i YlabelS; "" 

5336  PRINT   "CP  0,0" 

5340  ELSE 

5344  PRINT    "SP0«SP2" 

5348  PRINT   "SR    1.2.2.4iPU  PA  -8.35|DI   O.liLBUiPR   1.0.5»LBoiPR  -1,0.5iLB   (kU/n 

5352  PRINT   "PR  -1.0.5iSR    1.1.5jLB2iSR    1.5,2iPR    .5,.5»LB.»PR    .S.OiLBK)" 
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5356  PRINT  "PA4Z.-10;DI  1,0iLBV;PR  .4,-1tLBuiPR  1  ,  .  Si  LB<  w/e  >' 

5360  PRINT  "SP0;SP1" 

53G4  END  IF 

5368  Ipn-0 

537Z  BEEP 

5376  INPUT  "UANT  TO  PLOT  DATA  FROM  A  FILE  < 1-Y ,0-N>?" .Okp 

538©  Icn-0 

5384  IF  Okp-1  THEN 

5388  BEEP 

539Z  INPUT  "ENTER  THE  NAME  OF  THE  DATA  FILE" ,D_f lleS 

53361  IF  Iuo-0  THEN 

5400  BEEP 

5404  INPUT  "SELECT  (0-LINEAR.  1-L06( X . Y )" , Ilog 

5408!  END  IF 

541Z  ASSIGN  «File  TO  D_file« 

5416  BEEP 

54Z0  INPUT  "ENTER  THE  BE6INNIN6  RUN  NUMBER", Md 

54Z4  BEEP 

54Z8  INPUT   "ENTER   THE   NUMBER  OF   X-Y  PAIRS   STORED" .Npaira 

543Z  IF   Iuo-0  AND   Ihrat-0  THEN 

5436  BEEP 

5440  INPUT   "ENTER  DESIRED  HEAT   FLUX",Q 

5444  END  IF 

5448  BEEP 

545Z  PRINTER   IS    1 

5456  PRINT  USING  "4X. ""Select   a  aynbol : 

5460  PRINT  USING   "4X.""1      Star        2      Plue   aign 

5464  PRINT   USING   "4X,""3     Circle  4      Square 

5468  PRINT   USIN6   "4X,""5     Ronbua 

5472  PRINT   USING   "4X.""6     Right-eide-up  triangle 

5476  PRINT  USING   "4X,""7     Up-eide-doun  triangle 

5480  INPUT   Syn 

5484  PRINTER   IS  705 

5488  PRINT   "PU  DI" 

5492  IF  Syn-1  THEN  PRINT  "SH«" 

5496  IF  Sym-Z  THEN  PRINT  "SH+" 

5500  IF  Syn-3  THEN  PRINT  "SHo" 

5504  Nn-4 

5508  IF  Ilog-1  THEN  Nn- I 

5512  IF  Hd>1  THEN 

5516  FOR  1-1  TO  (  Md- 1 ) 

5520  ENTER  CFileiXa.Ya 

5524  NEXT  I 

5528  END  IF 

5532  IF  Ihrat-0  THEN 

5536  Ql-Q 

5540  IF  Ilog-1  THEN  O-LOG(Q) 

5544  END  IF 

5548  FOR  1-1  TO  Npaira 

5552  IF  Iuo-0  AND  Ihrat-0  THEN 

5556  ENTER  «File«Xa.B< • ) 

5560  Ya-B(0) 

5564  FOR  K-1  TO  Nn 

5568  Ya-Ya+B(K)»Q*K 

5572  NEXT  K 

5576  END  IF 

5580  IF  Iuo-1  OR  Ihrat-1  THEN 

5584  ENTER  CFilesXa.Ya 

5588  IF  Iuo-1  THEN  Ya-Ya/1000 

5592  END  IF 
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5S96  IF  Iuo-0  AND  Ihrat-0  THEN 

5600  IF  Ilog-1  THEN  Ya-EXP(Yo) 

5604  IF  Ilog-0  THEN  Ya-Ql/Ya 

5608  IF  Irt-0  THEN 

5612  IF  Xa-0  THEN 

5616  Yo-Ya 

5GZ0  Ya-1 

5624  ELSE 

5628  Ya-Ya/Yo 

5632  END  IF 

5636  ELSE 

5640  Han-FNHsnooth(Q.Xa.lBat) 

5644  Ya-Ya/Hen 

5648  END  IF 

5652  END  IF 

5656  Xx<I-1)-Xa 

5660  Yy(I-1)-Ya 

5664  X-<Xa-X«ln>»Sfx 

5668  Y-(Ya-Y«in)«Sfy 

5672  IF  Sy«>3  THEN  PRINT  "SH" 

5676  IF  Sy«<4  THEN  PRINT  "SR  1.4.2.4" 

5680  PRINT  "PA",X.Y."" 

5684  IF  Syn>3  THEN  PRINT  "SR  1.2,1.6" 

5688  IF  Sy«-4  THEN  PRINT  "002.4.99,0.-8,-4.0,0.8.4.0.1" 

5692  IF  Sy«-5  THEN  PRINT  "UC3. 0.99 .-3.-6.-3.6. 3.6.3. -6i" 

5696  IF  Syn-6  THEN  PRINT  "UC0.5. 3,99 . 3 .-B.-6.0. 3,8|" 

5700  IF  Sy*-7  THEN  PRINT  "UC0.-5. 3.99 .-3.8.6,0.-3,-8»" 

5704  NEXT  I 

5708  BEEP 

5712  ASSI6N  •File  TO  • 

5716  END  IF 

5720  PRINT  "PU  SM" 

5724  BEEP 

5728  INPUT  "UANT  TO  PLOT  A  POLYNOMIAL  < 1-Y .0-N>7" ,0kp 

5732  IF  Okp-1  THEN 

5736  BEEP 

5740  PRINTER  IS  1 

5744  PRINT  USING  "4X. ""Select  line  type 

5748  PRINT  USIN6  "6X,""0     Solid  line 

5752  PRINT  USIN6  "6X.""1     Dashed 

5756  PRINT  USING  "6X,""2...S  Longer  line  -  daah 

5760  INPUT  Ipn 

5764  PRINTER  IS  70S 

5768  BEEP 

5772  INPUT  "SELECT  (0-LINEAR.  1-L06(  X.Y) )"  .Hog 

5776  Iprnt-1 

5780  CALL  Poly 

5784  IF  Iuo-1  THEN 

5788  BEEP 

5792  INPUT  "DESIRE  TO  SET  X  Lower  and  Upper  Liwit  (  1-Y.0-N^7" .Ixlin 

5796  IF  Ixli«-0  THEN 

5800  XI 1-0 

5804  Xul-7 

5808  END  IF 

5812  IF  Ixlin-1  THEN 

5816  BEEP 

5820  INPUT   "ENTER  X   Louer    Linlf.Xll 

5824  BEEP 

5828  INPUT  "ENTER  X  Upper  Li«it".Xul 

5832  END  IF 
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5836  END  IF 

5840  FOR  Xa-Xll  TO  Xul  STEP  X«tep/ZS 

5844  Icn-Icn+1 

5848  Ye-FNPoly<Xa) 

5852  IF  Iuo-1  THEN  Ya-Ya/1000 

585B  Y«< Ya-Y«in>»Sfy 

5860  X-<Xa-Xwin)«Sfx 

5864  IF  Y<0  THEN  Y-0 

5868  IF  Y>100  THEN  60T0  5908 

5872  Pu-e 

5876  IF  Ipn-1  THEN  Idf-Icn  NOD  Z 

5880  IF  Ipn-2  THEN  Idf-Icn  MOD  4 

5884  IF  Ipn-3  THEN  Idf-Icn  HOD  8 

5888  IF  Ipn-4  THEN  Idf-Icn  MOD  16 

5892  IF  Ipn-5  THEN  Idf-Icn  HOD  32 

5896  IF  Idf-1  THEN  Pu-1 

5900  IF  Pu-0  THEN  PRINT  "PA" .X.Y,"PD" 

5904  IF  Pu-1  THEN  PRINT  "Pfl" .X.Y."PU" 

5908  NEXT  Xa 

5912  PRINT  "PU" 

5916  60T0  5372 

5920  END  IF 

S924  BEEP 

5928  INPUT  "UflNT  TO  QUIT  ( 1-Y.0-N)?" .Iqult 

5932  IF  Iqult- 1  THEN  5940 

5936  60T0  5372 

5940  PRINT  "PU  SP0" 

5944  SUBEND 

5948  SUB  Stats 

5952  PRINTER  IS  701 

5956  J-0 

5960  K-0 

5964  BEEP 

5968  IF  Iplot-1  THEN  ASSI6N  (File  TO  P.fllet 

5972  BEEP 

5976  INPUT  "LAST  RUN  No7(0-QUIT>" ,Nn 

5980  IF  Nn-0  THEN  6124 

5984  Nn-Nn-J 

5988  Sx-0 

5992  Sy-0 

5996  Sz-0 

6000  Sxs-0 

6004  Sys-0 

6008  Sze-0 

6012  FOR  1-1  TO  Nn 

6016  J-J+1 

6020  ENTER  CFileiQ.T 

6024  H-Q/T 

6028  Sx-Sx+Q 

6032  Sxe-Sxs+CTZ 

6036  Sy-Sy+T 

6040  Sya-Sy8+T*Z 

6044  Sz-Sz+H 

6048  Sza-5z8+H-2 

6052  NEXT  I 

6056  Oave-Sx/Nn 

6060  Tave-Sy/Nn 

6064  Have-Sz/Nn 

6068  Sdevq-SQR(  flBS<  <  Nn»SxB-Sx*Z  >/<  Nn»<  Nn- 1 ) ) ) ) 

607Z  Sdevt-SOR<  flBS(  ( Nn»Sya-Sy*Z >/<  Nn»<  Nn- 1 ) ) ) ) 
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6076  Sdevh-SQR<ABS<(Nn»Szs-Sz*2>/<Nn«(Nn-1 )>>) 

6®60  Sh-100«Sdevh/Have 

6084  Sq-100»Sdevq/Qave 

6088  St-100»Sdevt/Tave 

6092  IF  K-1  THEN  El  16 

609B  PRINT 

6100  PRINT  USING  " II X , ""DATA  FILE: "" . 14A" i Filet 

6104  PRINT 

6108  PRINT  USING  " 1 IX, ""RUN  Htube     SdevH   Qdp        SdevQ   Thetab  SdevT"" 

6112  K-1 

6116  PRINT  USING  " 1 1 X .DD,2< 2X.D. 3DE . IX, 3D. 2D> ,2X ,00. 3D, 1 X.3D. 2D" \ J .Have.Sh.Qave 

.Sq.Tave.St 

6120  60T0  B972 

6124  ASSIGN  BFilel  TO  • 

6128  PRINTER  IS  1 

6132  SUBEND 

6136  SUB  Coef 

6140  COM  /Cply/  A( 10.10), C( 10). B(5>.N.Iprnt,0po,Ilog,Ifn.IJoin.Njoin 

6144  BEEP 

6148  INPUT  "6IVE  A  NAME  FOR  CROSS-PLOT  FILE",Cpf« 

6162  CREATE  BOAT  Cpf»,2 

6156  ASSI6N  (File  TO  Cpf« 

6160  BEEP 

6164  INPUT   "SELECT   (0-LINEAR, 1-L06(  X, Y) )" ,Ilog 

6168  BEEP 

6172  INPUT   "ENTER  OIL  PERCENT   < - 1-STOP )" .Bop 

6176  IF  Bop<0  THEN  6192 

6180  CALL  Poly 

6184  OUTPUT   CFileiBop.BC  > 

6188  GOTO  B1B8 

6192  ASSI6N  "File  TO  • 

6196  SUBEND 

6200  SUB  UileonCCf ,Ci ) 

6204  COM  /Uil/  02.01. Do. L.Lu.Kcu 

6208  Din  Enf<12) 

62 12 1  ULISON  PLOT  SUBROUTINE  DETERMINE  CF  AND  CI 

6216  BEEP 

6220  INPUT  "ENTER  DATA  FILE  NAME". File* 

6224  BEEP 

6228  PRINTER  IS    1 

6232  PRINT  USING   "4X. ""Select   option: 

6236  PRINT  USING   "4X.""      9  Vary  Cf   and     Cl""" 

6240  PRINT   USING   "4X.""       1    Fix      Cf    Vary  Ci 

6244  PRINT   USING   "4X,""      2    Vary   Cf    Fix      Ci 

6248  INPUT   "ENTER  OPTION" .Icfix 

6252  PRINTER  IS  701 

625B  IF   Icfix-0  THEN  6272 

6260  IF   lcfix>0  THEN  BEEP 

6264  IF   Icfix-1    THEN   INPUT   "ENTER  Cf.Cef 

6268  IF   Icfix-2   THEN   INPUT   "ENTER  CI".Ci 

6272  PRINTER  IS    1 

6276  INPUT   "Want    To  Vary  Coef f ?< 1-Y.0-N)" .Iccoef 

6280  IF    Iccoef-1    THEN   INPUT   "ENTER  COEFF'.R 

6284  PRINTER  IS  701 

6288  IF   Icfix-0  OR   Icfix-2    THEN  Cfa-.004 

6292  IF  Icflx-1  THEN  Cfa-Cef 

629B  Cla-Ci 

6300  Xn-.8 

B304  Fp-.3 
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6308  Jj-» 

6312  Rr-3. 

B31B  IF    Iccoef-I    THEN   Rr-R 

6320  PRINTER   IS    1 

6324  PRINT   Do.Di.Kcu 

6328  ASSIGN   CFile   TO   Filet 

6332  ENTER  ftTilejNrun.Dt  .Ldtcl .Ldtc2 .Itt 

6336  Ru-Do*L06<Do/Di>/(2»Kcu) 

6340  S*-0 

6344  Sy-0 

6348  Sxy-0 

6352  SxZ-0 

6356  Sy2-0 

6360  FOR   1-1    TO  Nrun 

6364  ENTER   «F ilei Bop . TIHEtC Told ) ,Enf  <•> ,F«t 

63681  CONVERT  EHF'S  TO   TEMPERATURE 

6372  FOR  J- I    TO  5 

6376  T<J>-FNTvsv<E*f<J>> 

6380  NEXT  J 

6384  Tent- (T(  1  )  +  T(2))t.5 

6388  Tavg-T(5) 

6392  6r ad- 37. 9853+. 1 04388* Tavg 

6396  Tdrop-Enf <7>«1 .E+6/( 10.»6rad) 

6400  Tavgc-T<S>-Tdrop«.5 

6404  IF   ABS<Tavg-Tavgc>>.01    THEN 

6408  Tavg-(Tavg+Tavgc)«.5 

6412  60T0  6392 

6416  END  IF 

6420! 

64241  Conpute  properties  of  uater 

6428  Kw-FNKu(Tavg) 

6432  Muua-FNHuu<  Tavg ) 

6436  Cpu-FNCpu(TavQ) 

6440  Pru-FNPr u(  Tavg ) 

6444  Rhou-FNRhou(Tovg) 

64481 

64521  Compute  properties  of  Fr eon- 1)4 

6456  L«td-Tdrop/L06((T<5)-Tsat)/(T(5)-Tdrop-Tsat>) 

6460  IF  JJ-O  THEN 

6464  Tu-Tsat+Fr«LKtd 

6468  Thetab-Ty-Tsat 

6472  Jj-1 

6476  END  IF 

6480  Tf-<Tu*Tset>».S 

6484  Rho-FNRho< Tf  > 

6488  Mu-FNMu( Tf ) 

649Z  K-FNK(Tf) 

6496  Cp-FNCp<  Tf ) 

6500  Beta-FNBeta<Tf > 

6504  Hf g-FNHf g( Tsat  > 

6508  Ni-Hu/Rho 

6512  Alpho-K/<Rho«Cp> 

6516  Pr-Ni/fllpha 

65201 

65241  Analysis 

65281  COMPUTE  HOOT 

6532  A-PI»<0o*2-Dl*2>/4 

6536  P-PI»Oo 

6540  Hdot-3.9657E-3+F«e»<3.619S5E-3-F*s»<8.82O06E-6-F«s»<  1 .Z3688E-7-F«s»4. 31897 

E-10))) 
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6544  Q-Mdo-t»Cpu*Tdrop 

65481  COMPUTE  NATURAL-CONVECTIVE  HEAT-TRANSFER  COEFFICIENT 

65521  FOR  UNENHANCED  END(S) 

6556  Hbar-190 

6560  Fe-<Hbar»P/(Kcu«A))*.5»Lu 

6564  Tenh-FNTanh<Fe) 

6568  Theta-Thetab«Tanh/Fe 

6572  Xx-(9.B1»Beta»Thetab»Do*3»Tanh/<Fe»Ni»Alpha>r .166667 

6576  Yy-< 1+< .559/Pr>-<9/16))*<8/27) 

658©  Hbarc-K/Do»<  . 6+ . 387»Xx/Yy >* 2 

6584  IF    ABS<<Hbar-Hbarc>/Hbar>>.©01    THEN 

6588  Hbar-<Hbar+Hbarc>».5 

6592  60T0  6560 

6596  END  IF 

6600 1 

66041  COMPUTE  HEAT  LOSS  RATE  THR0U6H  UNENHANCED  ENDS 

6608  QI-(Hbar»P»Kcu»Ar.S»Thetab»Tanh 

6612  Oc-0-2»01 

6616  A6-PI»D2«L 

66201  COMPUTE  ACTUAL  HEAT  FLUX 

6624  Odp-Qc/Ae 

6628  IF   Icfix-0  OR  Icfix>1    THEN  C»f-1/Cf *< 1 ./Rr > 

6632  Thetab-CBf/Cp«Hfg»<Qdp/<Hu»Hfg)«< .014/< 9.81»Rho) )* .5)*< 1/Rr)»Pr*1 .7 

6636  Ho-Qdp/Thetab 

6640  Onega-Ho/Cf 

6644  Uo-Q/<PI»Do«L»Lwtd> 

6648  Vy-Mdot/(Rhou»PI»Di'2/4) 

6652  RewRhou*Vw*Di/Huua 

6656  Twl-Tw40*Ru/(PI*Do«L) 

6660  6ana-Ku/Di*Rew* . 8*Pru* (1/3. )•(  Muua/FNMuw(  Tyi ) ) * . 14 

6664!  PRINTER  IS  1 

6668  Yw-(  l./Uo-Ry )• Onega 

6672  Xy0nega*Do/(6ana*Di> 

6676  Sx-Sx+Xu 

6680  Sy-Sy+Yu 

6684  Sxy-Sxy*Yu»Xy 

6688  Sx2-Sx2+Xu«Xy 

6692  Sy2-Sy2+Yw»Yy 

6696  NEXT  I 

6700  ASSI6N  CFile  TO  • 

6704  H-<Sx»Sy-Nrur>*Sxy)/(Sx«Sx-Nrun«SxZ> 

6708  C-(Sy-Sx»M)/Nrun 

6712  IF  Icfix-0  OR  Icfix-3  OR  Icfix-4  THEN 

6716  Clc-1/M 

6720  Cfc-1/C 

6724  END  IF 

6728  IF  Icfix-1  THEN 

6732  Cic-1/M 

6736  Cfc-Cf 

6740  END  IF 

6744  IF  Icfix-Z  THEN 

6748  Clc-Ci 

6752  Cfc-1/C 

6756  END  IF 

6760  IF  ABS((Ci-Cic)/Cic)>.001  OR  ABS(  ( Cf-Cf c >/Cf c> > .001  THEN 

6764  Ci-<Ci+Cic)«.5 

6768  Cf-<Cf*Cfc)».S 

6772  PRINTER  IS  1 

6776  PRINT  USING  "ZX.M"  Cef  -  "" ,HZ. 3DE ,ZX.""  Ci  -  "" .HZ. 3DE" «C»f .Ci 

6780  PRINTER   IS  701 
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6784  GOTO  G328 

6788  END  IF 

G79Z  PRINT 

679E  PRINTER  IS  701 

680©  PRINT  USING  "Z3X,""   Cf  Ci 

6804  PRINT  USING  "8X." "ASSUMED      "" .HZ . 3DE , 3X .HZ . 3DE" %U a .Cia 

6808  PRINT  USING  "8X , ""CALCULATED   "" ,MZ. 3DE .3X.HZ. 3DE" ?Cof ,Ci 

68 1 Z  PRINT 

6816  Sui"iZ-SyZ-Z«M»Sxy-Z»OSy+rrZ»SxZ+Z»H«OSx+Nrun»C*Z 

68Z0  PRINT   USING   "  10X  .""Sum   of    Squares   -   ""  ,Z.  3DE"  jSuinZ 

68Z4  PRINT   USING   "10X .""Coef f icient    -   "" .D.DOO" «Rr 

68Z8  SUBEND 

683Z  DEF   FNHuu(T) 

6836  A-Z47.8/(T+133.I5> 

6840  Hu-Z.4E-5»10*A 

6844  RETURN   Hu 

6848  FNEND 

685Z  DEF   FNCpw(T) 

6856  Cpu-4.ZnZ0858-T«<Z.Z68Z6E-3-T«<4.4Z361E-542.7UZ8E-7»T)> 

6860  RETURN  Cpw*1000 

6864  FNEND 

6868  DEF   FNRhou(T) 

687Z  Ro-999.5Z946*T«<  .01Z69-T»(5.482513E-3-T« 1 .Z34147E-5) > 

6876  RETURN   Ro 

6880  FNEND 

6884  DEF   FNPrw(T) 

6888  Pry-FNCpu<  T ) »FNHuw< T ) /FNKu<  T ) 

6B9Z  RETURN  Pru 

6896  FNEND 

6900  DEF   FNKw(T) 

B904  X-(T+Z73.15)/Z73.15 

6908  Kw— .9ZZ474X«(Z. 8395-Xt( 1.8007- X*( .52577- .07344«X) ) ) 

69 1Z  RETURN  Kw 

6916  FNEND 

69Z0  SUB  Plot 

B9Z4  CON  /Cply/   A< 10, 10) ,C< 10) ,B< 5) .Nop.Iprnt  .Opo.Ilog.If n.I join.Njoin 

69Z8  DIM  Bs(3) 

B93Z  INTEGER  It 

6936  PRINTER  IS  1 

6940  Idv-0 

6944  BEEP 

6948  INPUT  "LIKE  DEFAULT  VALUES  FOR  PLOT  (  1-Y.0-N)?" .Idv 

69 5Z  Opo-e 

6956  BEEP 

6960  PRINT  USING  "4X. ""Select  Option 

6964  PRINT  USING  "6X.M"0  q  versus  delta-T 

6968  PRINT  USING  "6X.""1   h  vereue  delta-T 

697Z  PRINT  USING  "6X.-'Z   h  vereue  q 

6976  INPUT  Opo 

6980  BEEP 

6984  INPUT  "SELECT  UNITS  <0-SI . 1 -ENGLISH)" .Iun 

6988  PRINTER  IS  705 

699Z  IF  IdvOI  THEN 

6996  BEEP 

7000  INPUT  "ENTER  NUHBER  OF  CYCLES  FOR  X-AXIS", Cx 

7004  BEEP 

7008  INPUT  "ENTER  NUHBER  OF  CYCLES  FOR  Y-AXIS", Cy 

70 1Z  BEEP 

7016  INPUT  "ENTER  HIN  X-VALUE  (HULTIPLE  OF  10)".X«in 

70Z0  BEEP 
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7024  INPUT  "ENTER  MN  Y-VALUE  (MULTIPLE  OF  10>".Ynin 

702  8  ELSE 

7032  IF  Opo-0  THEN 

7036  Cy-3 

7048  Cx-3 

7044  X«in-.1 

7048  Ynin-100 

7052  END  IF 

7056  IF  Opo-I  THEN 

7060  Cy-2 

7064  Cx-2 

7068  Xfiin-.l 

7072  Ywin-100 

7076  END  IF 

7080  IF  Opo-2  THEN 

7081  IF  Iun-0  THEN 
7084  Cy-3 

7088  Cx-2 

7092  Xmin-1000 

7096  Ymn-100 

7097  ELSE 

7098  Cy-3 

7099  Cx-3 

7100  XHln-100 

7101  Ynin-10 

7102  END  IF 

7104  END  IF 

7105  END  IF 
7108  BEEP 

7112  PRINT   "INiSPhlP   2300. 2200. 8300. 6B00»" 

7116  PRINT  "SC  0. 100.0. IOOiTL  2.0;" 

7120  Sfx-100/Cx 

7124  Sfy-100/Cy 

7128  BEEP 

7132  INPUT  "UANT  TO  BY-PASS  CAGE?  < 1-Y.O-N)" .Ibyp 

7136  IF  Ibyp-1  THEN  7632 

7140  PRINT  "PU  0.0  PD" 

7144  Nn-9 

7148  FOR  1-1  TO  Cx+1 

7152  Xat-Xwin»10*<I-1> 

7156  IF  I-Cx+1  THEN  Nn-1 

7160  FOR  J-1  TO  Nn 

71B4  IF  J-1  THEN  PRINT  "TL  2  0" 

7168  IF  J-2  THEN  PRINT  "TL  1  0" 

7172  Xa«Xat«J 

7176  X-L6T<Xa/Xwin>«Sfx 

7180  PRINT  "Pft'iX.-.Oi  XTj" 

7184  NEXT  J 

7188  NEXT  I 

7192  PRINT  "PA  lOO.OiPUt" 

7196  PRINT  "PU  PA  0.0  PD" 

7200  Nn-9 

7204  FOR  1-1  TO  Cy+1 

7208  Yat-Y«in»10*(I-1 ) 

7212  IF  I-Cy-H  THEN  Nn-1 

7216  FOR  J-1  TO  Nn 

7220  IF  J-1  THEN  PRINT  "TL  2  0" 

7224  IF  J-2  THEN  PRINT  "TL  1  0" 

7228  Ya-YaW 

7232  Y-LGT<Ya/Y«in>»Sfy 
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7Z36 

PRINT    "PA  0."»Y."YT" 

724® 

NEXT    J 

7244 

NEXT    I 

7248 

PRINT    "PA  0.100  TL    0  2" 

7252 

Nn-9 

72  56 

FOR   1-1    TO   Cx+1 

72  60 

Xat-Xnin»10*(I-1 ) 

7264 

IF    I-Cx+1    THEN  Nn-1 

7268 

FOR  J-1    TO  Nn 

7272 

IF    J-1    THEN  PRINT    "TL   0  Z" 

7276 

IF   J>1    THEN  PRINT    "TL   0    1" 

7280 

Xa-XaW 

7284 

X-L6T(Xa/X«in)»Sfx 

7288 

PRINT   "PA" j X,", 100;    XT" 

7292 

NEXT   J 

7296 

NEXT   I 

730® 

PRINT   "PA    100.100  PU  PA    100,0  PD" 

7304 

Nn-9 

7308 

FOR  1-1    TO  Cy+1 

7312 

Ya+.-Ynin«10*<I-1> 

7316 

IF   I-Cy+1    THEN  Nn-1 

732® 

FOR  J-1    TO  Nn 

7324 

IF   J-1    THEN  PRINT   "TL  0  2" 

73Z8 

IF   J>1    THEN  PRINT   "TL  0   1" 

7332 

Ya-Ya+,»J 

7336 

Y-L6T<Ya/Y«in)«Sfy 

7340 

PRINT   "PD  PA    100." ,Y,"YT" 

7344 

NEXT  J 

7348 

NEXT   I 

7352 

PRINT   "PA    100.100  PU" 

7356 

PRINT   "PA  0.-2    SR    1.5.2" 

7360 

Ii-L6T(Xnin) 

7364 

FOR   1-1    TO  Cx+1 

7368 

Xa-X«in«10*<I-1) 

7372 

X-L6T<Xa/X«in>«Sfx 

7376 

PRINT   "PA"iX,".0i" 

7380 

IF   Ii>-0  THEN  PRINT   "CP   -Z ,-2iLB10iPR  -Z.ZiLB" 

7384 

IF    IK0  THEN  PRINT   "CP   -Z  .-ZiLB10;PR  O.Z»LB"iI 

7388 

Ii-Ii+1 

7392 

NEXT   I 

7396 

PRINT   "PU  PA  0,0" 

7400 

Ii-L6T<Y«in) 

7404 

Y10-1O 

7408 

FOR   1-1    TO  Cy+1 

7412 

Ya-Ynin«10*(I-1) 

7416 

Y-L6T<Ya/Y«in)»Sfy 

7420 

PRINT   "PA  0,"|Y."" 

7424 

PRINT   "CP  -4,-.Z5iLB10iPR  -Z ,2iLB"iIii"" 

7428 

Ii-Ii+1 

7432 

NEXT   I 

7436 

BEEP 

7440 

INPUT   "UANT  USE   DEFAULT   LABELS   (  1-Y.O-N)?" ,Idl 

7444 

IF   IdlOl    THEN 

7448 

BEEP 

7452 

INPUT   "ENTER  X-LABEL" .Xlabeli 

7456 

BEEP 

7460 

INPUT   "ENTER  Y- LABEL " .Ylabelt 

7464 

END  IF 

7468 

IF   Opo<2   THEN 

7472 

PRINT   "SR   I.ZiPU  PA  40,-14i" 
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7476  PRINT  "LB<T;PR  -1.6,3  PD  PR  1.2,0  PUiPR  . 5 ,-4; LBuoiPR  .5,1i" 

7480  PRINT  "LB-TjPR  .5.- 1 1 LBsat iPR  .5.1;" 

7484  IF  Iun-0  THEN 

7488  PRINT  "LB)  <KP 

7492  ELSE 

7496  PRINT  "LB)  <F)" 

7500  END  IF 

7504  END  IF 

7508  IF  Opo-2  THEN 

7512  IF  Iun-0  THEN 

7516  PRINT    "SR    1.5,2»PU  PA  40.-14iLBq   (U/puSR    I.I.SjPR  0.5,1iLB2»SR    1.5.2iPR 

0.5,-ltLB)" 

7520  ELSE 

7524  PRINT   "SR    1.5.2|PU  PA    34.-14iLBq   (Btu/hnPR    .5,.5iLB.iPR    .5, -.5]" 

7528  PRINT  "LBft;PR  .5,1»SR  1,1.5(LB2)SR  1.5.2iPR  .5,-1jLB>;" 

7532  END  IF 

7536  END  IF 

7540  IF  Opo-0  THEN 

7544  IF  Iun-0  THEN 

7548  PRINT   "SR    1.5.2|PU  PA   -12,40iDI    0.1iLBq   <U/pmPR  -1,0.5jSR    1.1.5|LB2;SR    1 

.5,Z«PR    1 ,.5«LB)" 

7552  ELSE 

7556  PRINT  "SR  1.5,2;PU  PA  -  1 2 , 32 « DI  0.1iLBq  <Btu/hr»PR  -.5. .5»LB. jPR  .5,.5j" 

7560  PRINT  "LBftjSR  I.I.SiPR  -I..S»LB2|PR  1,.5iSR  I.S.ZlLB)" 

7564  END  IF 

7568  END  IF 

7572  IF  Opo>0  THEN 

7576  IF  Iun-0  THEN 

75801  PRINT  "SR  I.S.ZtPU  PA  -1Z.38iOI  0,1iLBh  <U/«iPR  -1,.5«SR  1.1.5iLBZtSR  1. 

5.2iPR  .5..5i" 

7584  PRINT  "SR  t.Z,Z.4»PU  PA  -12,37jDI  0.1iLBhsPR  1.0.5tLBoiPR  -1.0.5iLB  (U/n 

7588  PRINT  "PR  -1,.5»SR  I.I.SiLBZiSR  1.5,Z|PR  .S..5iLB.iPR  .5.0iLBK)" 

7592  ELSE 

7596  PRINT   "SR    I.S.ZtPU  PA   -IZ.ZBiOI   0.1iLBh   (Btu/hrjPR  -.S,.5|L8. iPR   .5,.Sj" 

7B00  PRINT   "LBfiiPR  -1,.5tSR   I.I.SiLBZiSR    I.S.ZiPR   .S..5iLB.iPR   ,5..5jLBF)" 

7604  END   IF 

7608  END  IF 

7612  IF    Idl-0  THEN 

7616  PRINT   "SR   I.S.ZtPU  PA  50,-16  CP"s-LEN( Xlabelt)/2i "0iLB" iXlabelti"" 

7620  PRINT   "PA  -14,50  CP  0." |-LEN< YlabelS>/2«5/6i"DI    0, 1 iLB"i Ylabelti"" 

7624  PRINT   "CP   0,0  DI" 

7628  END  IF 

7632  Ipn-0 

7636  X11-1.E+6 

7640  Xul—  l.E+6 

7644  Icn-0 

7648  Hn-0 

7652  Ijoin-1 

7656  BEEP 

7660  INPUT   "WANT   TO  PLOT  DATA  FROM  A  FILE   (  1-Y.0-N)?" ,0k 

7664  IF   Ok-1    THEN 

7668  BEEP 

7672  INPUT   "ENTER  THE   NAME   OF   THE   DATA  FILE" ,D_f lie* 

7676  ASSI6N  "File   TO  D.fileS 

7680  BEEP 

7684  BEEP 

7688  INPUT   "ENTER  THE   BE6INNING  RUN  NUMBER". Md 

7692  BEEP 

7696  INPUT  "ENTER  THE  NUMBER  OF  X-Y  PAIRS  STORED" .Npairs 
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77031 

7704! 

7708 

7712 

771B 

7720 

772  4 

7728 

7732 

773B 

7740 

7744 

7748 

7752 

7756 

7760 

7764 

7768 

7772 

7776 

7780 

7784 

7788 

7792 

7796 

7800 

7804 

7808 

7812 

7816 

7820 

7824 

7828 

7832 

7836 

7840 

7844 

7848 

7852 

7856 

7860 

7864 

7868 

7872 

7876 

7880 

7884 

7888 

7892 

7896 

7900 

7904 

7908 

7912 

7916 

7920 

7924 

7928 

7932 

7936 


"CONNECT  DATA  UITH  LINE  < 1-Y ,0-N)7" . Icl 


■6X, 
"6X, 


'"  1 
""3 
'"5 
"'6 
'"7 


Star 

Circle 

RopibuB""" 

Right-Bide-up  triangle""" 

Up-side-doun  triangle 


"snv 

"SMo1 


BEEP 

INPUT 

BEEP 

PRINTER  IS  1 

PRINT  USING  "4X. ""Select  a  oynbol 

PRINT  USING 

PRINT  USING 

PRINT  USING  "6X 

PRINT  USING  "6X 

PRINT  USING  "6X 

INPUT  Syn 

PRINTER  IS  705 

PRINT  "PU  DI" 

IF  Syw-1  THEN  PRINT  "SM» 

IF  Sy«-2  THEN  PRINT 

IF  Sy«-3  THEN  PRINT 

IF  Md>l  THEN 

FOR  1-1  TO  (Hd-1 ) 

ENTER  CFileiYa.Xa 

NEXT  I 

END  IF 

FOR  1-1  TO  Npeira 

ENTER  SFileiYa.Xa 

IF  1-1  THEN  Ql-Ya 

IF  I-Npairs  THEN  Q2-Ya 

IF  Opo-1  THEN  Ya-Ya/Xa 

IF  Opo-2  THEN 

O-Ya 

Ya-Ya/Xa 

Xa-Q 

END  IF 

IF  Xa<Xll  THEN  XI 1 

IF  Xa>Xul  THEN  Xul 

IF  Iun-1  THEN 

IF  Opo<2  THEN  Xa-Xa*1.8 

IF  Opo>0  THEN  Ya-Ya».1761 

IF  Opo-0  THEN  Ye«Ya» 

IF  Opo-2  THEN  Xa-Xa» 

END  IF 

X-L6T(Xa/Xnin)*Sfx 

Y-L6T(Ya/Y«in)«Sfy 

Kj-0 

CALL  Synb(X,Y,Syn,Icl,Kj> 

60T0  7924 

IF  Sy«>3  THEN  PRINT 

IF  Sy«<4  THEN  PRINT 

IF  Icl-0  THEN 

PRINT  "PA'.X.Y."" 

ELSE 

PRINT  "PA"IX,Y,"PD" 

END  IF 

IF  Sy«>3  THEN  PRINT  "SR  1.2.1.6 

Sy«-4  THEN  PRINT 

Sy«-S  THEN  PRINT 
IF   Sy«-6   THEN  PRINT 
IF   Syw-7  THEN  PRINT 
NEXT   I 
PRINT   "PU" 
BEEP 
INPUT   "WANT   TO  LABEL?    < 1-Y.O-N)" .Ilab 


Plus   8ign" 
Square 


■Xa 
■Xa 


,317 
.317 


•SH" 
'SR 


1.4.2.4' 


IF 
IF 


"UC2 ,4. 99, 0.-8. -4. 0.0.8.4.0«" 
"UC3.0.99.-3.-6.-3.6.3.6.3.-6I' 
"UC0.5.3.99.3.-8.-6.0.3.8i" 
"UCO.-5.3.99.-3.B.6.0.-3,-8j" 
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794®  IF  Ilab-1  THEN 

7944  PRINT  "SP0tSP2" 

7948  BEEP 

7952  IF  Klab-0  THEN 

795G  Xlab-5 

7960  Ylab-85 

7964  INPUT  "ENTER  INITIAL  X.Y  LOCATIONS" , XI ab.Ylab 

7968  Xtt-Xlab-S 

7972  Ytt-Ylab+8 

7976  PRINT  "SR  1,1.5" 

798®  PRINT  "SM|PA"tXtt,Ytt,"LB       X   Heat   File" 

7984  Ytt-Ytt-3 

7988  PRINT  "PA" .Xtt , Ytt , "LB      Oil   Flux   Name" 

7992  IF  Sy«-1  THEN  PRINT  "SM»" 

7996  IF  Syn-2  THEN  PRINT  "SM+" 

8000  IF  Syw-3  THEN  PRINT  "SHo" 

8004  Klob-1 

8008  END  IF 

8012  Kj-1 

8016  CALL  Synb(Xlab.Ylab,Syn.Icl.Kj ) 

8020  PRINT  "SR  1 . I.SjSH" 

8024  IF  Sy«<4  THEN  PRINT  "PR  2.0" 

8028  BEEP 

8032  INPUT   "ENTER  BOP". Bop 

8036  IF   Bop<10  THEN  PRINT   "PR  2 ,0«LB" iBopi "" 

8040  IF   Bop>9   THEN  PRINT   "PR   .S,0tLB" iBopi "" 

8044  Ihf-0 

8048  IF  01 >Q2  THEN  Ihf-I 

8052  IF  Ihf-0  THEN  PRINT  "PR  4.0iLBInc" 

8056  IF  Ihf-1  THEN  PRINT  "PR  4.0jLBDec" 

8060  PRINT  "PR  2,0jLB"»D_file»j"" 

8064  PRINT   "SP0;SP1|SR    1.5.2" 

8068  Ylab-Ylab-5 

8072  END  IF 

8076  BEEP 

8080  ASSI6N  IF  lie  TO  • 

8084  Xll-Xll/1.2 

8088  Xul-Xul»1.2 

8092!  60T0  8040 

8096  END  IF 

8100  PRINT  "PU  SM" 

8104  BEEP 

8108  INPUT  "WANT  TO  PLOT  A  POLYNOMIAL  (  1-Y.0-N)?" ,6o_on 

8112  IF  6o_on-1  THEN 

8116  BEEP 

8120  PRINTER  IS  1 

8124  PRINT  USING  "4X, ""Select  line  type 

8128  PRINT  USIN6  "GX.""0     Solid  line""" 

8132  PRINT  USIN6  "6X.""1     Oaehed""" 

8136  PRINT  USIN6  "BX.""Z , , ,5  Longer  line  -  daeh""" 

8140  INPUT  Ipn 

8144  PRINTER  IS  705 

8148  BEEP 

8152  INPUT  "SELECT  (0-LIN. 1-L06(  X.Y) )" .Ilog 

8156  Iprnt-1 

8160  CALL  Poly 

8164  IF  Ifn-0  THEN 

8168  BEEP 

8172  INPUT   "ENTER  NUMBER  OF   FILES  TO  JOIN?".Njoln 

8176  END  IF 
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8180  Ijoin-0 

8184  IF    Ifn<Njovn   THEN   Ijoln-1 

8188  IF    Ifn-0  OR   Ijoin-1    THEN 

8132  FOR  Ij-0  TO   3 

819B  Be(Ij )-Bs<Ij )+B(Ij ) 

8200  NEXT   Ij 

B204  Ifn-Ifn+1 

8208  END  IF 

8212  IF  Njoin-Ifn  THEN 

8216  FOR  Ij-0  TO  3 

8220  B<Ij)-Bs<Ij )/Njoin 

8224  Bs(Ij>-« 

8228  NEXT  Ij 

8232  Ifn-8 

823B  ELSE 

8240  GOTO  7B56 

8244  END  IF 

8248  BEEP 

8252  INPUT  "ENTER  Y  LOUER  AND  UPPER  LIMITS" ,Y11 ,Yul 

8256  FOR  Xx-0  TO  Cx  STEP  Cx/200 

8260  Xo-Xnin«1®-Xx 

8264  IF  Xa<Xll  OR  Xa>Xul  THEN  8372 

B268  Icn-Icn+1 

827Z  Pu-0 

8276  IF  Ipn-1  THEN  Idf-Icn  HOD  Z 

8280  IF  Ipn-2  THEN  Idf-Icn  HOD  4 

8284  IF  Ipn-3  THEN  Idf-Icn  HOD  8 

8288  IF  Ipn-4  THEN  Idf-Icn  HOD  16 

6292  IF  Ipn-5  THEN  Idf-Icn  HOD  28 

8296  IF  Idf-1  THEN  Pu-1 

8300  IF  Opo-0  THEN  Ya-FNPoly( Xa) 

8304  IF  Opo-2  AND  Ilog-0  THEN  Ya-Xa/FNPoly(  Xa) 

8308  IF  Opo-2  AND  Ilog-1  THEN  Ya-FNPoly(  Xa) 

8312  IF  Opo-1  THEN  Ya-FNPoly< Xa) 

8316  IF  Ya<Y«in  THEN  8372 

8320  IF  Ya<Yll  OR  Ya>Yul  THEN  8372 

8324  IF  Iun-1  THEN 

8328  IF  Opo<2  THEN  Xa-Xa* 1 . 8 

8332  IF  Opo>0  THEN  Ya-Ya«.1761 

8336  IF  Opo-0  THEN  Ya-Ya*.317 

8340  IF  Opo-2  THEN  Xa-Xa*. 317 

8344  END  IF 

8348  Y-L6T<Ya/Y*in)*Sfy 

8352  X-L6T<Xa/X*in)«Sfx 

8356  IF  Y<0  THEN  Y-0 

83B0  IF  Y>100  THEN  GOTO  8372 

8364  IF  Pu-0  THEN  PRINT  "PA" ,X,Y,MPD" 

8368  IF  Pu-1  THEN  PRINT  "PA" tX.Y,"PU" 

8372  NEXT  Xx 

8376  PRINT  "PU" 

8380  60T0  7656 

8384  END  IF 

8388  BEEP 

8392  INPUT  "UANT  TO  PLOT  REILLY'S  DATA?  (  1-Y.0-N)" .Irly 

8396  IF  Opo-0  OR  Opo-1  THEN 

8400  Xll-3 

8404  Xul-20 

8408  END  IF 

8412  IF  Opo-2  THEN 

6416  XI 1-10000 
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84Z0  Xul- 100000 

84Z4  ENO  IF 

84Z8  IF  Irly-I  THEN 

843Z  Y11-Z0 

8436  Yul-70 

8440  BEEP 

8444  INPUT  "ENTER  LOWER  AND  UPPER  Y-LIMTS  FOR  PL0TTIN6"  .Yll .  Yul 

8448  FOR  Xx-0  TO  Cx  STEP  Cx/Z00 

B45Z  Xa-Xnln*10'Xx 

8456  IF  Xa<Xll  OR  Xa>Xul  THEN  8508 

84E0  Xl-LOG(Xa) 

8464  IF  Opo-0  THEN  Yl  — Z  .S40Z837E-1 -fXl^  4.97Z01  51 -X1»Z  .51  34787E- 1  > 

8468  IF   Opo-1    THEN  Yl— Z.540Z837E- 1+X1K 3. 97Z01 S1-X1«Z .51 34787E- 1 ) 

847Z  IF  Opo-Z  THEN  Yl  —  3.7073801E- 1 +X1»(  B.7ZS9 190E- 1 -Xl»6.88Z684ZE-3> 

8476  Ya-EXP(Yl) 

8480  Y-L6T< Ya/Y«in)«Sfy 

8484  X-L6T<Xa/Xmn)«Sfx 

8488  Ipu-0 

B49Z  IF  Y<Y11  THEN  Ipu-1 

8496  IF  Y>Yul  THEN  60T0  8508 

8500  IF  Ipu-0  THEN  PRINT  "PA" ,X.Y,"PD" 

8504  IF  Ipu-1  THEN  PRINT  "PA" ,X.Y,"PU" 

6508  NEXT  Xx 

851 Z  PRINT  "PU" 

8516  END  IF 

85Z0  BEEP 

85Z4  INPUT  "WANT  TO  PLOT  ROHSENOU  CORRELATION?  (  1-Y.O-N)" .Iroha 

85Z8  IF  Iroh»-t  THEN 

853Z  Yll-15 

8536  Yul -80 

8540  BEEP 

8544  INPUT  "ENTER  T»at  ( Dag  C)",Taat 

8548  Cef-.0040 

855Z  BEEP 

8556  INPUT  "ENTER  Cef  ( DEF-0.004 )" .Csf 

8560  Tf-Teat+Z 

8564  FOR  Xx-0  TO  Cx  STEP  Cx/Z00 

8568  Xa-Xnin*10~Xx 

857Z  IF  Xa<Xll  OR  Xa>Xul  THEN  8684 

8576  Xl-L06<Xa) 

8580  IF   Opo<Z    THEN  Tf-Tsat+Xa/Z 

8584  Rho-FNRho<  Tf ) 

8588  K-FNK(Tf) 

859Z  Mu-FNMu(Tf) 

8596  Cp-FNCp( Tf ) 

8600  Hf Q-FNHf g<  Teal ) 

6604  Ni-Mu/Rho 

6608  Pr-Cp«Hu/K 

861 Z  OneQa-C8f«Hfg/Cp»<(.014/<9.81»Rho))*.5/<Hu«HfO))*(1./3)«Pr*1.7 

8616  IF  Opo-0  THEN  Ya-( Xa/0naga)*3 

B6Z0  IF  0po-1  THEN  Ya-( Xa/0«ega)*3/Xa 

86Z4  IF  Opo-Z  THEN  Yb-Xb*(  2  .  /3)/0P.ega 

86Z6  IF  Opo-Z  THEN 

863Z  Tf  c- Tsat  +  Xa/ Ya« . 5 

8636  IF  ABS(Tf-Tfc)>.01  THEN 

8640  Tf-(Tf+Tfc)».S 

8644  GOTO  8594 

8648  END  IF 

865Z  END  IF 

6656  Y-L6T<Ya/Y«in)«Sfy 
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8660  A-L6T(Xo/X«in)»Sfx 

Jb64  Ipu-0 

8668  IF  Y<Y11  THEN  Ipu-1 

867Z  IF  Y>Yul  THEN  8684 

8676  IF  Ipu-0  THEN  PRINT  "PA" .X , Y ."PD" 

8680  IF  Ipu-1  THEN  PRINT  "PA" ,X, Y."PU" 

8684  NEXT  Xx 

8688  PRINT  "PIT 

869 Z  END  IF 

8696  BEEP 

8700  INPUT  "WANT  TO  QUIT  < 1 -Y ,0-N)" , Iqt 

8704  IF  Iq-t-1  THEN  87 1 Z 

8708  GOTO  7636 

871 Z  PRINT  "PU  PA  0,0  SPO" 

8716  SUBEND 

B7Z0  SUB  SyRb<X,Y,Syn,Icl,Kj  ) 

87Z4  IF  Syw>3  THEN  PRINT  "SM" 

87Z8  IF  Sy«<4  THEN  PRINT  "SR  1.4.2.4" 

8732  Yad-0 

8736  IF  Kj-1  THEN  Yad-.8 

8740  IF  Icl-0  THEN 

8744  PRINT  "PA"  .X.  Y-»Yad."" 

8748  ELSE 

8752  PRINT  "PA" ,X , Y+Yad."PO" 

8756  END  IF 

8760  IF  Sy«>3  THEN  PRINT  "SR  1.2.1.6" 

8764  IF  Sy«-4  THEN  PRINT  "UC2 .4 .99 . 0.-8. -4. 0.0. 8. 4. Oi " 

8768  IF  Sy*-5  THEN  PRINT  "UC3.0.99 ,-3, -6, -3,6,3,6,3,-61" 

8772  IF  Syw-B  THEN  PRINT  "UCO. 5. 3,99 .3.-8.-6.0.3.81" 

877B  IF  Syw-7  THEN  PRINT  "UCO. -5. 3. 99 .-3.8,6.0.-3.-81" 

8780  IF  Kj-1  THEN  PRINT  "SHiPR  O.-.B" 

8784  SUBEND 

8788  SUB  Purg 

8792  BEEP 

8796  INPUT  "ENTER  FILE  NAME  TO  BE  DELETED" ,File» 

8800  PURGE  File* 

8804  60T0  8792 

8808  SUBEND 

8812  SUB  Tdcn 

8816  COM  /Cc/  C(7).Ical 

8820  01 H  Enf(l) 

8824  DATA  0. 1008B091 .25727.94369 .-767345. 8295. 78025595. 81 

8828  DATA  -9247486589 .6. 9768BE-H  1 .-2. G6192E+ 13.3. 94078E+ 14 

8832  READ  CO) 

8836  BEEP 

8840  INPUT  "GIVE  A  NAME  FOR  FILE  TO  BE  CREATED" .File* 

8844  BEEP 

8848  INPUT  "SELECT  TUBE  (O-UH. 1-HF ,2-UT)" .Iti 

8852  BEEP 

8856  INPUT  "SELECT  THERMOCOUPLE  TYPE  (O-NEU, 1 -OLD)" .leal 

8860  IF  I«<2  THEN  Dl-,0127 

8864  CREATE  BDAT  File*. 4 

8868  ASSI6N  CFlIe  TO  File* 

8872  OUTPUT  ■Fileiltt 

8876  J-0 

8880  BEEP 

8884  OUTPUT  709 i "TIMEDATE" 

8888  ENTER  709 iDi 

8900  PRINTER  IS  1 

8904  PRINT 


205 


8908  PRINT  "  Month,  date  and  tine:  "  \  DATES<  Dt ) ,TIME«<  Dt ) 

8912  PRINT 

891E  PRINT  USING  "10X,""  Fne      Tin      Tev      Vu*2      Tdrop 

892®  IF  K-0  THEN 

8924  PRINTER  IS  701 

8926  PRINT 

8932  PRINT  "  Month,  date  and  time:  " ;DATE»< Dt ) .TIMES*  Dt ) 

8936  IF  Itt-0  THEN  PRINT  USIN6  " 10X.""Tube  Type:  Uieland  Smooth 

894®  IF  Itt-1  THEN  PRINT  USING  " 10X .""Tube  Type:  High  Flux 

8944  IF  Itt-2  THEN  PRINT  USING  " 10X .""Tube  Type:  Turbo-B 

8948  PRINT 

8952  PRINT  USING  "10X.""  Fns      Tin      Tev      Vy"2      Tdrop""" 

8956  PRINTER  IS  1 

896®  K-1 

8964  END  IF 

8968  BEEP 

8972  INPUT  "ENTER  FLOUMETER  READING" , Fne 

8976  OUTPUT   709t"AR  AF0  AL4   VR1" 

8980  FOR  L-0  TO  4 

8984  OUTPUT  709 i" AS  SA" 

8988  IF  L>0  AND  L<4  THEN  9020 

8992  S-0 

8996  FOR  1-0  TO  9 

9000  ENTER  709 i E 

9004  S-S+E 

9008  NEXT  I 

3012  IF  1-0  THEN  E«f< 0)-ABS( S/ 10) 

9016  IF  1-4  THEN  E«f ( 1 )-ABS< S/ 10> 

9020  NEXT  L 

9024  OUTPUT   709 i "AR  AF20  AL20  VR1" 

9028  OUTPUT  709 t "AS   SA" 

9032  Etp-0 

9036  FOR   1-0  TO   9 

9040  ENTER  709 .Et 

9044  Etp-Etp4Et 

9048  NEXT    I 

9052  Etp-Etp/10 

9056  Tin-FNTv»v(E*f< 1 >> 

9060  Tev-FNTvev(E*f<0>) 

9064  6rad-37.98S3+.104388»Tin 

9068  Mdot-3.96S7E-3+F«e»<3.6195SE-3-F«e»<8.8ZO06E-6-F«e«< 1 .Z36B8E-7-F*a«4. 31897 

E-10))) 

9072  Vu-Hdot/< 1000»PI»Di*Z>»4 

9076  Tdrop-Etp»1.E4B/< 10»Grad) 

9080  PRINT  USING   "10X, 3<DD.DO,4X> . 1X.Z.00.4X.MZ.4D" iF*e.Tin.Tev.Vy*Z .Tdrop 

9084  BEEP 

9088  INPUT  "UANT  TO  ACCEPT  THIS  DATA  SET?  < 1-Y.0-N)" .Ok 

9092  J-J-H 

9096  IF  Ok-0  THEN 

9100  J-J-1 

9104  60T0  8968 

9108  ELSE 

9112  OUTPUT  •FilejF«e.E«f<»),Etp 

9116  PRINTER  IS  701 

9120  PRINT  USING  " 10X . 3< DO.DO.4X) , 1X.Z.D0.4X.MZ .*D" l Fne .Tin, Tev , Vu*Z .Tdrop 

9124  PRINTER  IS  1 

91Z8  BEEP 

913Z  INPUT   "WILL   THERE   BE   ANOTHER  DATA  SET?   (  1-Y.0-N)" .Go_on 

9136  IF  6o_on-1    THEN  8968 

9140  END  IF 
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9144  ASSIGN  BFile  TO  • 

9148  PRINTER  IS  701 

9152  PRINT 

9156  PRINT  USING  " 10X .""NOTE:  "".ZZ,""  data  sets  are  stored  In  file  "",15A"|J.F 

lie* 

9  IB©  PRINTER  IS  1 

9164  SUBENO 

9168  SUB  Uoprt 

917Z  PRINTER  IS  1 

9176  BEEP 

91 BC  INPUT   "Enter   Uo   File   Na«e",File* 

91 84  BEEP 

9188  INPUT  "Number  of  Data  Rune" .Nrun 

919Z  INPUT  "Do  You  want  a  Ploi  File?< 1-Y.©-N)" , Iplot 

9196  BEEP 

320©  IF    Iplot-1    THEN 

9204  INPUT   "6ive  Plot    File   Name" ,P_f ile« 

9208  CREATE   BOAT  P_file«,4 

9212  ASSI6N  ©Plot    TO  P_file« 

9Z16  END  IF 

9220  PRINTER  IS  701 

9224  PRINT 

9228  PRINT 

9232  PRINT  USING  "10X.""    Uater  Vel        Uo 

9236  ASSI6N  CFile  TO  FlleS 

9240  IF  Iplot-1  THEN  ASSIGN  OFilel  TO  P_f ile» 

9244  FOR  I- I  TO  Nrun 

9248  ENTER  BFileiVw.Uo 

9252  IF  Iplot-1  THEN  OUTPUT  OFileliVy.Uo 

9256  PRINT  USING  " 15X,0.D0.6X.H2.3DE"j Vu.Uo 

9260  NEXT  I 

9264  ASSIGN  (File  TO  • 

9268  ASSI6N  OFllel  TO  • 

9272  PRINT  USING  "1 OX, ""NOTE:  "".ZZ,""  data  sets  are  stored  in  file  "".15A"iNru 

n, Filet 

9276  IF  Iplot-1  THEN 

9280  PRINT  USING  " 10X."MNOTE:  "".ZZ.""  X-Y  Pairs  are  stored  in  file  "" . »5A" ;Nru 

n.P.filet 

9284  END  IF 

9288  PRINTER  IS  1 

9292  SUBEND 

3296  SUB  Select 

330©  COM  /Idp/  Idp 

9304  BEEP 

9308  PRINTER  IS    1 

9312  PRINT  USING   "4X, ""Select   option: 

3316  PRINT  USING   "6X.""   0  Taking  data  or   re-processing  previous   data 

3320  PRINT  USING   "6X.""    1    Plotting  data  on  Log-Log   

3324  PRINT  USING   "6X.""    2   Plotting  data  on  Linear 

3328  PRINT   USING   "6X,""    3  Hake  cross-plot    coefft    file""" 

3332  PRINT  USIN6  "6X.""   4  Re-circulate  yater 

3336  PRINT  USING   "6X,""   5  Purge 

3340  PRINT  USING  "6X.""   6   T-Drop  correction 

3344  PRINT  USING   "BX,""   7  Print   Uo  File 

3345  PRINT  USING   "BX,""    8  Modify   X-Y   file 

3346  PRINT  USING   "6X.""    9   Move 

3347  PRINT  USING  "6X.""10  Co*b/Fixup 

3348  INPUT   Idp 

3352  IF  Idp-0  THEN  CALL  Main 

3356  IF  Idp-1  THEN  CALL  Plot 
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9360  IF  Idp-2  THEN  CALL  PI  in 

9364  IF  I dp- 3  THEN  CALL  Coef 

93B8  IF  Idp-4  THEN  CALL  Main 

937Z  IF  Idp-5  THEN  CALL  Purg 

9376  IF  Idp-B  THEN  CALL  Tdcn 

938®  IF  Idp-7  THEN  CALL  Uoprt 

936)  IF  Idp-8  THEN  CALL  Xywod 

9382  IF    Idp-9    THEN  CALL    Move 

9383  IF    Idp-10  THEN  CALL   Conb 
9385  SUBENO 

9394  SUB   Xynod 

9434  PRINTER   IS    1 

9444  BEEP 

9454  INPUT  "ENTER  FILE  NAME" .File* 

9464  ASSI6N  «File1  TO  File* 

9474  BEEP 

9484  INPUT  "ENTER  NUMBER  OF  X-Y  PAIRS", Np 

9494  BEEP 

9504  INPUT   "ENTER  NEU  FILE   NAHE'.FileZ* 

9514  CREATE   BOAT   FileZS.5 

95Z4  ASSI6N  IFileZ   TO  FileZ* 

9534  BEEP 

9544  INPUT  "ENTER  NUMBER  OF   X-Y  PAIRS  TO  BE   DELETED" .Ndel 

9554  IF   Ndel-0  THEN  9594 

9564  FOR  1-1    TO  Ndel 

9574  BEEP 

9584  INPUT  "ENTER  DATA  SET  NUMBER  TO  BE  DELETED" ,Nd< I > 

9594  NEXT  I 

9604  FOR  J-1  TO  Np 

3614  ENTER  OFileljX.Y 

96Z4  FOR  1-1  TO  Ndel 

9634  IF  Nd< I )-J  THEN  9674 

9644  NEXT  I 

9654  OUTPUT  GFileZiX.Y 

9664  PRINT  J.X.Y 

9674  NEXT  J 

9675  PRINTER  IS  701 
9684  ASSI6N  SFllel  TO  • 
9694  ASSI6N  BFileZ  TO  • 
9704  SUBEND 

9714  SUB  Hove 

9718!  FILE  NAME:  HOVE 

97ZZI 

97Z6  DIM  Bop(66),A(B6),B(6G>.C(66),D<66),E(66),F(66),6<66).H(66),J<66>.X(66).L( 

66>.H<66> 

9730  DIM  Told(66).N(66),Vr(66),Ir(6G) 

9734  BEEP 

9738  INPUT   "OLD  FILE   TO  MOVE" ,DZ_fi let 

9742  ASSI6N  OFileZ   TO  OZ.filet 

9746  ENTER  •FileZiNrun.Dt .Ldicl .LdicZ ,Itt 

9750  FOR  I- I    TO  Nrun 

9754  ENTER  OFileZiBop< I ) ,Told( I ) 

9758  ENTER  «F-ileZiA(  I )  ,B<  I  >  ,C<  I )  ,D<  I  >  ,E<  I  >  ,F<  I )  ,6<  I  >  ,H<  I ) .  J<  I  >  ,K(  I )  ,L(  I )  ,M<  I )  ,N 

<I> 

976Z  ENTER  •FileZ»Vr(I).Ir<I) 

9766  NEXT   I 

9770  ASSI6N  «FileZ   TO   • 

9774  BEEP 

9778  INPUT  "SHIFT  DISK  AND  HIT  CONTINUE" ,0k 

978Z  BEEP 
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9786  INPUT   "INPUT  BDAT   SIZE", Size 

9790  CREATE   BOAT   DZ_f ile«,Slze 

9794  ASSIGN  «File1    TO  DZ_file« 

9798  OUTPUT  OFilel jNrun.Dt  .Ldtcl .LdtcZ ,Itt 

9B0Z  FOR  1-1    TO  Nrun 

9806  OUTPUT  OFilel ;Bop< I >,Told< I  ) 

9810  OUTPUT  •Flle1iA(I),B(I).C(I),0(I).E(I)lF(I),G(I).H(I),J(I),K(I),L(I),n<I>. 

N(I) 

9BU  OUTPUT  BFllel|Vr(I).Ir(I) 

9818  NEXT  I 

98ZZ  ASSIGN  CFlIel    TO   • 

9BZG!  RENAME    "TEST"   TO   DZ_file» 

3830  BEEP  ZO00..Z 

9834  BEEP  4000.. Z 

9838  BEEP   4000,. Z 

984Z  PRINT   "DATA  FILE  MOVED" 

9846  SUBEND 

9874  SUB  Corxb 

98781  FILE  NAME:  COMB 

988Z! 

9886  DIM  Enf(IZ) 

9890  BEEP 

9894  INPUT  "OLD  FILE  TO  FIXUP" ,DZ_f ilet 

9898  ASSI6N  CFileZ  TO  DZ.flle* 

990Z  D1_filet-"TEST" 

9906  CREATE  BOAT  DWllet.ZO 

9910  ASSI6N  BFilel    TO  Ol.fllet 

9914  ENTER  OFileZiNrun.Dt  .Ldicl .LdtcZ ,I*i 

9915  Nrunn-ZO 

9918  IF  K-0  THEN  OUTPUT  OFilel iNrunn.Dt .Ldtcl .LdtcZ .Itt 

99ZZ  FOR  1-1  TO  Nrun 

99ZG  ENTER  •FileZ;Bop.Told.E*f < • ) .Vr.Ir 

9930  OUTPUT  SFilel iBop, Told. Enf( •) .Vr.Ir 

9934  NEXT  I 

9938  ASSIGN  BFlleZ  TO  • 

994ZI  RENAME  "TEST"  TO  DZ_flle« 

9946  BEEP  ZOO0..Z 

9950  BEEP  4000.. Z 

9954  BEEP  4000.. Z 

9958  BEEP 

996Z  INPUT  "WANT  TO  ADO  ANOTHER  FILE  ( 1-Y.O-N)?" .Oka 

99G6  IF  Ok o-l  THEN 

9970  K-1 

9974  BEEP 

9978  INPUT  "GIVE  NEU  FILE  NAHE".Nfile$ 

998Z  ASSI6N  IFileZ  TO  Nfilet 

9986  60T0  9914 

9990  END  IF 

9994  ASSIGN  OFileZ  TO  • 

9998  SUBEND 
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